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long charges which can be 
loaded in even badly 
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Delay Units 
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The MFR-1, Universal Magnetic Tape 
Recorder, is designed as a portable drum 
type recorder without sacrificing any de- 
sirable features of truck mounted equip- 
ment. A stainless steel belt drive provides 
exacting speed control that exhibits itself 
in a 54 db signal to noise ratio—without 
utilizing noise cancelling. 


The bank of recording heads is installed 
on precision ways that afford easy access 
in conjunction with rigid and absolute 
geometric positioning. It is not necessary 
to disassemble the recorder to clean the 
heads. Each magnetic head is individually 
mounted with its attendant pressure 
spring and is easily replaced or aligned. 
The MFR-1 may be operated as a direct 
recorder, but in order to realize full ad- 
vantage of the exacting design features, 
the recorder should be operated in con- 
junction with FM accessories, 

The auxiliary power-pac, weighing less 
than 30 pounds and having a 2 cubic- 
foot volume, contains the time standard, 
phase comparator and motor drive oscil- 
lator with its attendant amplifier. Only 12 
volts at 4 amps is required to operate the 
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For information write Electro-Tech 
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rollers allow smooth, easy movement up 
or down when drawing lines. A spiral 
index window automatically indicates dis- 
tances between horizontal and angular 
lines as close as 1/16”. Lines can be 
spaced the same or varied distances 
apart. Vertical lines are made by placing 
ball point or pencil in any hole and roll- 
ing ruler up or down. Circles and arcs 


(Continued on page 31) 
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This is a reprint of the supplement published in 
GEOPHYSICS, v. 22, n. 1 (January 1957). 


GEOPHYSICAL CASE HISTORIES, 
VOLUME 1 (1948) 


$10.00 
Members of SEG, AAPG & AIME $ 6.00 


Foreign postage, per copy $ 0.50 


Edited by L. L. Nettleton. 680 pages. Fully 
illustrated. 7 x 10. Cloth. 


This volume, first of a series to be published at five 
to seven year intervals, is a collection of 60 papers 
by 61 authors on geophysical observations made 
under a wide variety of field circumstances. The 
purpose of the series is to provide material by which 
geophysical surveys can be judged from later devel- 
opment and thus aid in the interpretation and eval- 
uation of other geophysical work, Contents: 3 general 
and historical papers; 21 salt dome case histories 
(Texas, Louisiana and Mississippi); 17 mid-continent 
case histories (Arkansas, Illinois, Oklahoma and 
Texas); 4 Rocky Mountain case histories: 11 Califor- 
nia case histories; 4 foreign case histories. Features 
a complete title index cross-referenced to a classified 
index of all maps and figures. 


GEOPHYSICAL CASE HISTORIES, 
VOLUME II (1956) $10.00 
Members of SEG, AAPG & AIME $ 6.00 
Foreign postage, per copy $ 0.50 


Edited by Paul L. Lyons, Sinclair Oil & Gas 
Co. 690 pages. Fully illustrated. 7 x 10. 
Cloth. 


The second volume of the series contains 53 papers 
by 75 authors. In this volume the case histories are 
grouped by type of trap rather than by geologic 
province. It contains 17 foreign papers—13 more 
than Volume I, and more than half the material is 
original, being published for the first time. The 
second volume also contains a complete index cross- 
referenced to a classified index of all maps and 
figures. Contents: 6 general and historical, 6 salt 
dome, 6 reefs, 23 anticlines, 3 stratigraphic traps, 
3 mining and 5 new uses. 


CUMULATIVE INDEX, 1931-1953 (1955) $5.00 
Members of SEG $3.00 
Foreign postage, per copy $0.50 


Edited by Kenneth L. Cook, Department of 
Geophysics, University of Utah. 322 pages. 
7 x 10. Cloth. 


All publications of the Society (GEOPHYsICSs, EARLY 
GEOPHYSICAL PAPERS and GEOPHYSICAL CASE HISTORIES, 
VOLUME 1) during the period are indexed by subject, 
title and authors. Geophysical patents abstracted in 
GeopHysics from 1939 through 1953 are indexed 
separately by subject, patent number and inventor. 
Coverage of all literature reviews substantially ex- 
tends the usefulness of this index, since most of the 
significant literature of exploration geophysics since 
1936 not published by the Society has been reviewed 
in GEoPHysics. Whether or not the reader has a 
complete file of sec publications, he will find this 
index most useful. 
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SIXTH SUPPLEMENT $1.00 
Members of SEG $0.50 


This is a reprint of the supplement published in the 
1957 GEOPHYSICS YEARBOOK (December 1957). 


SEVENTH SUPPLEMENT $1.00 
Members of SEG $0.50 


This is a reprint of the supplement published in 
GEoPHysics, v. 23, n. 5 (December, 1958). 


LESSONS IN SEISMIC COMPUTING 


(1959) $5.50 
Members of SEG $4.25 
Foreign postage, per copy $0.50 


Written by M. M. Slotnick, edited by R. A. 
Geyer. 350 pages, 7 x 934. Fully illustrated. 
Cloth. 


An elementary text and problem book. Its 350 pages 
contain 44 lessons in seismology arranged for selec- 
tion or combination to cover the normal 36-week 
course, or for condensation into an 18-week course. 
Written by Dr. Slotnick for the Geophysics Depart- 
ment of Humble Oil & Refining Co., this series of 
lessons presents the basic physical laws relevant to 
seismic interpretation. The lessons begin without 
assuming more than secondary school mathematics. 
An elementary knowledge of the Calculus is desirable, 
though not required, for the last half of the book. In 
its scope, detail and clarity of style, this work pro- 
vides au authoritative reference for individual study, 
as weil as a textbook for formal courses. 


YEARBOOK of the Society of Explora- 
tion Geophysicists (Available Jan. 


1961) $3.00 
Members of SEG $2.50 
Foreign postage $0.20 


Published annually in January 


World-wide professional directory of exploration geo- 
physicists. Alphabetical listing of SEG membershi 
with the grade of membership indicated. Geographic 
listing by companies. Annual report of the SEG. New 
committees for the current year. Directory of SEG 
local sections. Listing of affiliated societies. Pictures 
of SEG officers. Report of the annual international 
meeting of SEG and abstracts of the papers. Consti- 
tution and by-laws, as amended. Calendar of future 
meetings. Code of ethics 


GEOPHYSICS (Bimonthly) 
Subscription, per year $10.00 
Foreign postage, per year $ .50 
Members of SEG No charge 


The official journal of the SEG is published in Feb- 
ruary, April, June, August, October and December 
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SEG PUBLICATIONS 


each year beginning 1960. (The first issue was pub- 
lished in January, 1936 as Volume I, No. 1. Publica- 
tion was quarterly through 1957. There were five 
issues in both 1958 and 1959.) Volumes run on a 
calendar year basis. An average of eight technical 
peas appear in each issue, all contributed by mem- 

$ and others engaged in exploration geophysics. 


GEOPHYSICS—Clothbound Volumes 


SEG Foreign 

Reprints Members Others Postage 
Vol. I (1936) $ 6.00 $10.00 $0.50 
Vol. 2 (1937) 6.00 10.00 0.50 
Vol. 3 (1938 6.00 10.00 0.50 
Vol. 4 (1939 6.00 10.00 0.50 
Vol 1940) 6.00 10.00 0.50 
Vol. 6 (1941) 6.00 10.00 0.50 
Vol. 7 (1942 8.00 12.00 0.50 
Vol. 8 (1943 8.00 12.00 0.50 
Vol. 9 (1944 8.00 12.00 0.50 
Vol. 10 ee 8.00 12.00 0.50 
Vol. II (1946 8.00 12.00 0.50 
Vol. 12 (1947) 8.00 12.00 0.50 
Vol. 13 (1948 10.00 14.00 0.50 
Vol. 14 (1949 10.00 14.00 0.50 
Vol. 15 (1950 10.00 14.00 0.56 
Vol. 16 (1951 10.00 14.00 0.56 
be 17 (1952 10.00 14.00 0.50 

18 (1953 10.00 14.00 0.50 
Orginal 

Vol. 22 (1957 6.00 10.00 0.50 
Vol. 23 (1958) 6.00 10.00 0.50 
Vol. 24 (1959) 6.00 10.00 0.50 
Vol. 25 (1960) 6.00 10.00 0.50 
Available issues, each $3.00 

Members of SEG $2.50 
Foreign postage $0.20 
Volume No. Issue Year 
19 1 January 1954 
19 2 April 1954 
19 3 Jul 1954 
19 4 October 1954 
20 | January 1955 
20 2 April 1955 
20 3 July 1955 
20 4 October 1955 
21 1 January 1956 
21 2 April 1956 
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23 i January 1958 
23 2 April 1958 
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23 4 October 1958 
23 5 December 1958 
24 I February 1959 
24 2 April 1959 
24 3 July 1959 
24 a October 1959 
24 5 December 1959 
25 I February 1960 
25 2 April 1960 
25 3 June 1960 
25 4 August 1960 
25 5 October 1960 
25 6 December 1960 
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GEOPHYSICS 


A CASE HISTORY OF THE BRONTE (ELLENBURGER) AND 
RAWLINGS FIELDS, COKE COUNTY, TEXAS* 


ROLLIN E. PHIPPSt 


ABSTRACT 


The North Bronte Area interested Lion Oil Company on the basis of surface, isopachous and 
magnetic studies. A large block of acreage was acquired and explored by the seismograph. A well 
was drilled on the shooting and proved the shooting to be reliable. A velocity survey in this well, 
however, indicated that the map of the reflective event identified as Ellenburger was actually Caddo. 
Additional shooting was done before the Lion No. 1 Edwards, the discovery well of the Bronte 
(Ellenburger) Field, was drilled. The discovery well of the Rawlings Field was drilled by Humble on 
seismic information. There are eight producing zones in these fields. Shooting costs and permit fees 
have remained amazingly constant through the years. 


INTRODUCTION 


The Bronte (Ellenburger) and Rawlings Fields are located in the eastern 
portion of Coke County, Texas, on the Eastern Shelf area of the Midland Basin. 
At present the Bronte Field represents the southernmost area of Cambrian pro- 
duction of the Cambrian Trend of West-Central Texas. Cambrian production 
was first established on this trend in the Bronte and Rawlings Fields. The Concho 
Arch is to the south and the axis of the Bend Arch is located 60 miles to the east. 
(See Figure 1.) 

GEOPHYSICAL PROCEDURES 


In 1946, Lion was offered a block of more than 11,000 acres. A surface struc- 
ture map, made by L. G. Welch in 1928, was submitted along with the acreage. 
A. E. Barnes was Lion’s District Geologist in Abilene at that time. A surface 
anomaly was indicated in the area and after additional study of the regional sub- 
surface structure of the Permian with special emphasis on Permian Isopachous 
intervals, Lion took the acreage. Mr. Barnes was convinced that the Permian 
Isopachous intervals were important and should indicate thinning or draping 
over, or the presence of, pre-Permian structures. A magnetic survey seemed to 
indicate a positive area substantiating the surface and subsurface feature. An 
Ellenberger subsurface study indicated nothing more than a structural ‘‘nosing” 
to the west across the block. 

* Presented at the 13th Annual Midwestern SEG Meeting, Shreveport, La., April 5, 1960. Manu- 

script received by the Editor May 26, 1960. 

t Lion Oil Company, Midland, Texas. 
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Fic. 1. Geologic Provinces West Texas & Southeast New Mexico. 


The initial seismic work was done from January to June, 1950. A single 
velocity survey was available in the Union No. 1 McCutchen, approximately 4 
miles southwest of the surveyed area. 

Figure 2 shows a portion of the seismic program conducted at this time. Ac- 
tually, the Bronte (Ellenburger) and Rawlings Fields are contiguous. The dis- 
covery well of the Bronte (Ellenburger) Field, the Lion No. 1 Edwards in the 
southeast quarter of Section 365, was completed February 11, 1953. The 
Rawlings Field discovery, the Humble No. 1 Rawlings in the southeast quarter 
of Section 397, was spudded in March 1952. The Bronte (Ellenburger) produces 
from the Ellenburger (Cambro-Ordovician), while the Rawlings Field produces 
from one Cambrian and six Pennsylvanian horizons. The Fort Chadbourne Field 
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BRONTE AND RAWLINGS FIELDS 


Fic. 2. Coke County, Texas, North Bronte Area. 
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is 43 miles northeast. It is a westward-dipping monoclinal structure and produces 
from the Lower Pennsylvanian (Odom and Caddo). The West Fort Chadbourne 
Field is 53 miles north and is Pennsylvanian production associated with the Fort 
Chadbourne fault system. The Bronte Field is 3 miles south and is another multi- 
pay structural field. 

The topography in most of the area was characterized by steep-sided buttes 
capped by Edwards Limestone. The buttes are shown on Figure 2, and only one 
shotpoint is located on a butte. Record quality was poor on the buttes, as shown 
on Figure 3. These records were shot less than 1,000 ft apart. The record at the top 
was shot on top of a 200-ft butte, and the lower one in the valley below. You will 
notice that the deeper energy is much better on the lower record. Needless to say, 
the crew stayed in the valleys. Tomlinson Geophysical Company shot 570 holes 
at an average depth of 78 ft, using 30 pounds of powder per shot, during the 
initial program. Five geophones per group were used, and slight ‘‘mixing’’ was 
deemed necessary. Permits cost $10.00 per hole. By computing the records to a 
horizontal datum plane, a general tilting of the seismic data to the southeast was 
noted. Rather than apply a velocity gradient, since so little velocity information 
was available, the records were recomputed and referred to the shot depth. This 
was effective in correcting the mis-ties. 

Good reflective events were tentatively identified as Vale shale, Saddle Creek 
lime, Palo Pinto lime, Lower Capps lime and Ellenburger dolomite. (See Correla- 
tion Chart, Figure 4.) The event identified as Palo Pinto (Lower Canyon) offered 
a very persistent continuous reflection in the southern portion of the area but 
deteriorated rapidly to the north. This deterioration is probably associated with 
a lime and reef buildup. The most persistent event on the records was the Lower 
Capps. Several leads were developed by the shooting. Probably due to a mis- 
correlation, an Ellenburger map, prepared by the geological staff while the shoot- 
ing was in progress, showed the presently producing area as an anomaly onthe 
down side of a fault. For this reason a rather small localized feature to the west 
was chosen as the best prospect. The Lion No. 1 Waldrop was drilled on this 
feature in October, 1950. It was located approximately 4 miles west of the Bronte 
(Ellenburger) Field, and northeast of the McCutchen (Swastika sand) Field. 
(See Figure 2.) The well checked the shooting at the Lower Capps (Upper Strawn) 
horizon within 75 ft, but the prospective Swastika sand of the Lower Cisco (Upper 
Pennsylvanian) was not present. A velocity survey was run in the Waldrop, and 
it was discovered that the event previously identified as the Ellenburger was 
actually Caddo. 

The results of the initial shooting in the field area are shown on Figure 5. 
This map was contoured on a reflective event identified as the Lower Capps, 
Lower Pennsylvanian in age. Actually the time horizon mapped is probably 
slightly above the Capps within the Upper Strawn. 

The importance of the isopach interval was recognized, and the crew was 


(Text continued on page 1174) 
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Fic. 4. Stratigraphic relationship of various units of the Permian and Pennsylvanian in West Texas. 
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Fic. 5. North Bronte, Rawlings Field, Coke County, Texas. Lower Capps. 
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requested to present the Saddle Creek to Lower Capps map shown on Figure 6, 
The discovery well on the Bronte (Ellenburger) Field is shown in the southeast 
quarter of Section 365. The separation of the two thin areas was proved by a dry 
hole in the northeast quarter of Section 390, plugged and abandoned July 7, 1953. 
A Lower Capps to Ellenburger Isopach was also prepared (Figure 7). The thin 
area is still persistent to the south; however, the thin area in Section 365 on Figure 
6 was not apparent on this interval. The inferred thickening to the east in Section 
397 is more abrupt than shown here since it is probably associated with faulting. 
The Ellenburger is eroded above the top of the Rawlings structure; but since 
the reflective event mapped was above the Ellenburger, this is a valid indication 
of Lower Pennsylvanian structure. 

The same crew was sent back for some additional seismic lines in May, 1951. 
The results of the additional lines across the now productive area are shown on 
Figure 8. The new line across Section 397 showed approximately 100 ft of addi- 
tional relief. The steep east dip, or down-to-the east faulting, on the east side of 
the field was indicated. Actually, the seismic interpretation would preclude any- 
thing but faulting. An offset well in the southwest quarter of Section 396, drilled 
by Humble in October of 1952, was 900 ft low to the discovery well. This was no 
surprise to Lion personnel, due to our seismic interpretation. Humble has credited 
the seismograph with the discovery of the Rawlings Field, the Bronte Field to 
the south, and the Fort Chadbourne to the north-east. As previously mentioned, 
the Fort Chadbourne Field is a monoclinal structure and was found due to a 
velocity anomaly which indicated a pseudo high closure. 

Figure 9 is a record section (A-A’) shot by Humble in April, 1948. These 
records were shot with 1,000-foot to 1,800-foot spreads, four geophones per group. 
Forty-pound charges were shot at approximately 120 ft. An uphole survey indi- 
cated a 9,500 ft/sec computing velocity, and the records were corrected to a 
1,600-ft datum plane. This cross-section runs east-west approximately bisecting 
the Rawlings Field. (See Figure 8.) 

Figure 10 is a record section (B-B’) which closely duplicates the location of 
the previous section. These records were shot for Lion by Tomlinson Geophysical 
Service in 1950. The method of shooting was described previously. These record 
sections are shown merely to indicate the record quality in the area. The present 
Humble and Lion supervisors would have preferred less mixing or compositing, 
but these records found oil. The shooting cost averaged $121.00 per location, or 
$13,800.00 per month, including explosives, permits, drilling and all other mis- 
cellaneous charges. 

Lion acquired gravity coverage in 1955 because of the “reef play” in the 
western portion of Coke County. Figure 11 shows the area as presently inter- 
preted from Exploration Survey’s derivative data. The lines of traverse of the 
gravity crew and the derivative data points are shown on this map. It will be 
noted that even with the knowledge of the anomaly, the gravity indication is very 

(Text continued on page 1178) 
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Fic, 6. Saddle Creek—Lower Capps Isopach. 
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Fic. 7. Lower Capps—Ellenburger Isopach. 
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Fic. 8. North Bronte, Rawlings Field, Coke County, Texas. Lower Capps. 
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subtle. This certainly demonstrates the need for detailed gravity control in an 
area such as this. Notice the zero gravity line. This is often indicative of faulting 
(all faults must be associated with a zero derivative contour, according to theory). 
Some substantiation is added to the seismic fault. 

Figure 12 is a revised composite of the seismic and subsurface interpretation 
on the Gardner sand (Strawn) after S. D. Ayers. This horizon is approximately 
300 ft below the lower Capps horizon, previously shown. Note the steep dip on 


< ib 


1179 


AND RAWLINGS FIELDS 


BRONTE 


ep 


piste 


SQYSS 


“Fe 
SE 


10. 


Fic. 


the present producing area of the Bronte 


igure 13 shows 


Fi 


ld. 


fie 


the east side of the 


demonstrates the narrow structural 


is map 


ields. Thi 
of 1953 all production 


(Ellenburger) and Rawlings I 


feature of the arez 


his area, previously ident 


s the ‘‘Rawlings”’ Fie 


int 


In May 
offic 


a. 


Ids, 


classified a 


ally be 


. 


” 


North Bronte 


fied as 


) 


(Text continued on page 1183 


|. 
7 
= . 
4 
ie 
. 


ROLLIN E. PHIPPS 


-1.3 
° 


° 


| 


Fic. 11. North Bronte, Rawlings Field, Coke County, Texas. Gravity. 
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An interesting and refreshing sidelight to the history of this seismic program is 
a comparison with present conditions in the area. A contract crew which shot in 
the vicinity last year reports that permits are still $10.00 to $15.00 per hole and 
the average cost per shot hole is $119.00. Midwest Oil Corporation reports that a 
program which was conducted in 1958, in an area just north of the subject area, 
had identical shooting conditions. The best record was produced with a 20-pound 
charge at about 100 ft, 1,320-ft spreads with six geophones per group. However, 
they had to pay up to $25.00 for permits. Midwest Oil Corporation conducted a 
seismic program west of the Rawlings Field in 1954. Permits were $25.00 per hole, 
and they shot 10-pound charges at 60-ft hole depths. Two-hundred feet of high 
velocity limestone build-up caused the drilling of the Midwest No. 1 Simpson, 
shown on Figure 7, 1 mile southwest of the Rawlings Field. The consistency of the 
shooting cost is probably due to a rapid exploration and “‘cooling”’ off of the 
area. In 1948, 100 crew weeks were recorded in Coke County, while 1950 showed 
only 65. By 1952 only 22 crew weeks were expended, and this figure dwindled 
to 18 in 1958. 

CONCLUSION 


These fields were found independently by seismograph crews working for 
Humble and Lion. They have now produced over 1,500,000 barrels of oil from 
19 wells. The initial shooting picture had to be detailed as much as topography 
would allow. Gravity could have indicated the lead; however, more detailed work 
would have been necessary to delineate the prospect area. The area has been a 


victim of rapid exploration and equally rapid ‘“‘cooling off”? period. 
I qually rap g 
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NUMERICAL RESISTIVITY INTERPRETATION: 
GENERAL INHOMOGENEITY* 


K. VOZOFFT 


ABSTRACT 


A linear approximation is developed for the equation of conduction in a medium where resistivity 
is an arbitrary function of x, y, and z. This is applied by assuming the earth to be subdivided into 
small, homogeneous blocks of arbitrary resistivity. Under this approximation, the apparent resistivity 
is just the sum of the effects of the individual blocks. The equations are linear, and surface apparent 
resistivity data can be inverted to yield block resistivities. 

The quality of the approximation has been checked by comparison with model measurements in 
two situations: remote current source (telluric method), and local current source (resistivity method). 
It was found that the results are satisfactory provided that the proper type of expression is used for 
the effect of the resistivity contrast of each block. 


The ultimate objective of this work is to enable a geophysicist to obtain a 
complete interpretation of a dc resistivity survey with a minimum of subjective 
adjusting. Ultimately, these techniques will probably develop to a stage where 
measurements will come in from the field, be put directly into a computer, and 
an interpretation will come out that uses not only all of the measurements to 
their full accuracy, but also any additional information which may be available, 
such as that from physical property measurements and geological observations. 
To some extent this is already being done with seismic data. The one-dimensional 
and two-dimensional resistivity problems have already been carried to the com- 
puter stage (Alfano, 1959; Ness, 1959; Vozoff, 1958). The following discussion 
applies to the three-dimensional problem. 


INTRODUCTION 


In order to do this problem, we replace the real, nonuniform earth with a 
model earth made up of small, uniform blocks. Then we can approximate any 
distribution of conductivity of interest by making the blocks small enough and 
by assigning to each of them an appropriate conductivity. We first set up the 
equation of conduction for this model earth, which allows us to predict the sur- 
face apparent resistivities for any assumed distribution of conductivity (forward 
problem). This was first done by Stevenson (1934) and the next section follows 
Stevenson quite closely. If we take the first two terms of the series, we have an 
approximation which is linear in that the effect of each block is added to give 
the net effect (no interaction). We are able to write a set of linear equations 
relating block conductivities with surface apparent resistivities. It is then a 


* Revised from a paper with the same title presented at the 29th Annual Meeting of the Society 
of Exploration Geophysicists, Los Angeles, California, November 11, 1959. Manuscript received 
by the Editor March 21, 1960. 

+ Department of Physics, University of Alberta, Edmonton, Alberta. 
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straightforward problem of working backward by matrix manipulation to find 
conductivities from apparent resistivity measurements (inverse problem). 

The last part of the work concerns some model measurements made to check 
the quality of the approximation, and an example of an actual solution is shown. 
EQUATION OF CONDUCTION 

Considering a conduction medium and direct current, at a point P where the 
conductivity is op and the vector current density is Jp, we have 
Ep = — Vor = Jr/op 
which defines the field at that point. Taking the divergence, we obtain 
V-Je = — Vor:Vor — opV’op. 
Rearranging, we obtain 
= — (V-Je + Vop)/or 


which is similar to Poisson’s equation. We solve for ¢p by integrating sources 
over all space, 


V-Jrt+ Vor: Vor 
TRIRP 


dR = dridrodr; = (r1 pi)? + p2)? + (rz ps)’. 


dR, 


For a source of current / located at point Q, the V- J term is zero everywhere ex- 
cept at Q, so the first term nae simply, and we have, assuming that op=ae, 


Vor: 
dee pr ORT RP 


rpg’ = (pi — qi)? + (p2 — G2)? + (ps — Qs)”. 


The first term is the well-known formula for potential in a uniform material of 
conductivity op and the second term is a correction because of the fact that the 
material is not uniform. This correction term is the potential due to charge dis- 
tributions which are induced wherever conductivity is changing. The size of the 
charge distribution is proportional to the gradient of the conductivity and to the 
strength of the electric field at that point. It is maintained by the current. If @ is 
constant, Vo=0 and the second term drops out. In the case of a half space with 
both P and Q on the surface, both terms are doubled and we have equation (1). 
We can call the first term the undisturbed potential, and the second, the cor- 
rection potential. 


Vor: Vor: Vor (1) 
or = dR 
TRIRP 


op “+ op’. 
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The assumption that op=og requires that there be a thin surface layer of uniform 
conductivity. 
THE APPROXIMATION 


Equation (1) would be a fairly routine integral equation for ¢ in terms of o 
were it not for the V@ under the integral. As it is, an exact solution in a finite 
number of terms cannot be written down. 

What is commonly done in a case like this is to approximate ¢r under the 
integral by the first of the two terms, i.e. for purposes of finding Vor we write 

I 


and take the gradient 
Vor 


Pr 


rer = (ri — + (r2 — + (rs — 


Substituting into the integral for dp’, we obtain 


I rar 
p rr 


Sr = In Or. 


The complete approximation for potential becomes equation (2). 


I Ter 
op = - dR}. (2) 
oe 


The process of successive approximation could continue ad infinitum, but it 
is of interest to see how good this approximation is. 

If we go back to the model made up of uniform cubes, the integrand in equa- 
tion (2) will be zero except at the interfaces between cubes, since within any 
cube a is constant. Assume that the cubes are to be oriented along our n, fe, 73 
axes, and that we are making measurements along, say, the r2 direction. For a 
surface A; (see Figure 1) perpendicular to the 7; axis, whose coordinates are 
(a1, 72, 13), 


~ 


rer = (a — + (re — g2)f + (ra — 
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gd Q 


SURFACE 


Fic. 1. Coordinate System and Physical Situation. 


The second term in equation (2), dp’, becomes 


p 


dS 
(a—gi) —— 
dr, 


which integrates to give us a function of conductivity multiplied by a function 
which depends only on geometry. 


—JAS 
op = f f 
p Ay 
drodrs 


+/AS 
= gi(P), 


p 


(3) 


where AS =In (02/01), the change in In o at the surface, and g:(P) represents the 
geometric contribution of that one surface. For a complete cube there will be one 
term for each surface, 


g(P)  TASG(P) 
p = 


p p 


The total potential is given by equation (4) where G(P) is the sum of the effects 
of the six surfaces of one cube. 
TAS 
op = gp? + G(P). (4) 


4a p 
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When there are » cubes, we can add up their effects linearly because changes in 
conductivity enter as the logarithm 


: > AS;G;(P). 


i=l 


op = op? + 
Since each of P and Q can be anywhere on the surface 


n 
$(P;, Or) ~ Or) + —— AS.Gi(P;, Qe), 


p 
or, simplifying notation, 


I n 
O° + AS jx (5) 
4r’op 

which is the most general form of the expression. This is a set of linear equations 
in m unknowns, the values of S;. Since we can take readings for as many locations 
P; and Q; as we choose, we can make the total number of readings equal to m and 
have a completely determined set of equations, or we can make the total number 
of readings greater than m and have an overdetermined set. In the latter case 
we can use least squares or some other form of fitting. 


QUALITY OF THE APPROXIMATION 


Consider again equation (4), which is the approximate potential at P due to 
a source at Q when there is one cube (0;) in an otherwise uniform half space(ep). 


4 I 
(=) G(P). (4) 


PT PQ op 


Compare this with the expression for a single sphere (o;) in a uniform full space 


(cp), 


nr,?"*! P,,(cos 8) 
sat [nor + (n + 


(6) 


where 7; is the radius of the sphere, and @ is the angle between reg and rep sub- 
tended at the center of the sphere. As o; becomes infinite or zero, equation (6) 
tends to some finite value, whereas equation (4) tends to plus or minus infinity. 
The behavior of (4) is contrary to experience, which indicates that once the ratio 
o1/op exceeds some moderate value, increasing it still further makes no detectable 
difference in the measured apparent resistivities. Likewise, a decrease in the ratio 
beyond moderate limits causes no detectable change in measured results. 

Assuming that the term In (0:/ap) of equation (4) could be replaced by a 
contrast factor of the form 
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Fic. 2. Variation of C with Conductivity Contrast. 


0, — Op 


—=C 7 


a series of model measurements were made to determine empirically the constant 
c. These measurements were made on single blocks of several conductivities in 
water at two different conductivities. This was done in two ways: with point 
sources of current, and with uniform current to simulate dc telluric fields. The 
results are summarized in Figure 2, and indicate that 


cm 3.6. 


Replacing equation (4) by equation (8) and using 


I €P 


pr pag 4a p + 20 p 


G(P) as calculated by the method of equation (3), the effects of several cubes 
were predicted and compared with measurements. Figures 3 and 4 show two 
typical examples. The data are shown as percentage difference in apparent re- 
sistivity, from that of a uniform medium due to the presence of the cube. Figure 
3 is the result in a uniform field. The inset at the right shows the model arrange- 
ment to the same horizontal scale as the curves. The traverse line passes over the 
center of a unit cube whose top surface is at a depth of one-half unit. The measure- 
ments were made on a graphite cube in tap water. 

Figure 4 is an example of a point-source measurement over the same cube 
at a depth of one unit. The presentation is the standard one for this type of 
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Fic. 3. Telluric Measurement. Percentage difference in field due to a cube. 
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Fic. 4. Apparent Resistivity Measurement. Percentage difference in apparent resistivity due 
to a cube. Pole-dipole configuration. 
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measurement (pole-dipole) and is discussed by Marshall and Madden (1959, 
p. 793). 

The solution might be improved by a method of successive approximations. 
That is, once a solution has been found, calculate a new potential and use this in 
equation (1). This has not yet been attempted. 

Figure 5 is an example of an analysis of model data, using this approach. It 
is taken from Ness (1959) with the kind permission of the author. In his thesis, 
Dr. Ness presents several examples of computer analysis of both model and field 
data, which indicate that the method is successful. The model (dark outline) was 
approximated by the group of eight smaller blocks shown. Blocks and model 
were six units long perpendicular to traverse, and the traverse line ran three units 

1 2 3 7 8 


Fic. 5. Result of Analysis of Model Measurements. Model (heavy outline) of conductivity 
3.8, can be replaced by eight smaller blocks of indicated conductivities. (After Ness, 1959). 


from the ends. The resistivities calculated from the approximation were fitted 
(least squares) to those measured on the model. The ratio o;/op was 3.8, for which 
equation (7) gives a contrast factor of 1.7, as compared to the calculated value 
of 1.6. 

LINEARITY OF THE APPROXIMATION 


We were able to put equations (4) and (5) into a form linear in the o; only 
because the effect at a boundary was, apparently, proportional to the logarithm 
of the ratio of the two conductivities separated by that boundary. We found 
later that the effect is more nearly proportional to an algebraic function of the 
sums and differences of the two conductivities, so our linearization is not really 
valid. 

We would expect that the approximation would be quite good for small devia- 
tions (o:=0p+e; e—0), and that it would deteriorate rapidly as | «/op| increases. 
Figure 6 shows the result of a model measurement in the worst possible situation: 
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Fic. 6. Telluric Measurement. Percentage difference due to a pair of cubes. 


a bakelite cube and a graphite cube in direct contact. From this it is apparent 
that, due to interaction, the approximation is not as good for a pair of cubes as 
it is for an isolated cube. It is improbable that the approximation can be, or 
should be, modified to compensate for this, and we must recognize that the 
solutions will be semiquantitative at best. 


DIMENSIONAL REQUIREMENTS ON DATA 


Regarding the data collection: it is not sufficient that the number of measure- 
ments equal the number of unknowns. As was indicated by Stevenson (1934), 
there must be a certain equivalence between the number of dimensions desired of 
the solution and the number of dimensions of field measurements taken.' That is, 
for the half-dimensional problem of horizontal layering, e=o(s); — <z <0, all 
we require for interpretation is a half-dimension of information, say an apparent 
resistivity curve, taken by leaving the source fixed in one location and moving 
the observation point along a curve away from the source to infinity. Vertical 
layering is a one-dimensional problem o=o(y); —~<y<™, and we need one 
full dimension of information for a complete solution, provided only that the 


1 Backus (1960), has proven this principle rigorously for the general three-dimensional problem 
under the restriction that o has only small deviations from its average over the half space. In his 
development, Backus uses the concept of a “generalized source,” i.e. a voltage applied to the entire 
surface, whose amplitude varies sinusoidally with r; and r2. The point source assumed throughout the 
present work is a special case of this generalized source, derived through a two-dimensional Fourier 
integral. 
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observer does not traverse parallel to strike. The case where anomalies are due 
to horizontal cylinders of various conductivities is a two-dimensional problem, 
and one way to solve it is to collect apparent resistivity curves for a very large 
number of source locations on a traverse line perpendicular to strike. This gives 
us a set of values which can be plotted and contoured like a map, a function of 
two variables. For the three-dimensional problem which we have been consider- 
ing, three dimensions of information are necessary. These can be obtained in 
principle by measuring potentials everywhere on the surface (two dimensions) 
for sources everywhere along a line (one dimension). In this case there is still one 
dimension that has not been used—the source is free to move everywhere on the 
surface, although we have restricted its location to one line. This extra information 
would be required if the cubes were oriented in some unknown directions other 
than the x, y, axes we used for our measurements. In that case we would need 
two additional numbers—the direction cosines of the actual cube axes—which 
could in principle be obtained by shifting the source from its original traverse. 
The ‘‘excess”’ of information could also be useful in the one-dimensional and two- 
dimensional problems with transverse anisotropy, provided that the axis of 
symmetry is horizontal. 


CONCLUSIONS 


One should not expect the solutions to be better than semi-quantitative, in 
view of the approximations involved. With this restriction, the method has been 
shown to be practical, and to yield satisfactory results. 

Several modifications might be tried. One possibility is a method of successive 
approximations to the potential. This should be tried, keeping in mind that a 
practical method will be a compromise between economics and quality. A second 
possible modification would be an alternative method of fitting, although it is not 
certain whether this is worthwhile with the present approximation. Work is con- 
tinuing on both possibilities. 

Alfano, in a recent paper (1959), presents a rigorous discussion of the problem 
of an earth made up of horizontal cylinders, and shows that the boundary-value 
problem is soluable exactly for certain special resistivity distributions. He then 
attacks the same problem, allowing an arbitrary distribution of cylinder resistivi- 
ties, and derives a solution which is very similar to that presented here. For 
geological situations in which there is a distinct strike (that is, where apparent 
resistivity changes are much more rapid in one direction than in another direction 
perpendicular to the first), the two dimensional approach is obviously to be pre- 
ferred, since it requires only the square root of the number of measurements and 
calculations required in the three-dimensional approach. 

The present work started as the author’s Ph.D. thesis, submitted to MIT in 
1956. It was supervised by Professor T. R. Madden, and was supported by a 
Morse Committee Fellowship for Computer Applications. Investigation of the 
problems was continued by Philip Hallof and Norman Ness at MIT, and more 
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recently by the author at the University of Alberta. Mr. S. A. Bukhari, a graduate 
student at Alberta, did many of the calculations for checking the model results. 
The work at Alberta has been supported in part by the Geological Survey of 
Canada. The writer is grateful for the assistance and support of each of these in- 
dividuals and groups. He is also indebted to Professor George Backus of MIT 
for a very original communication on the general problem. 
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THE REFLECTION OF RAYLEIGH WAVES BY A HIGH 
IMPEDANCE OBSTACLE ON A HALF-SPACE* 


R. W. FREDRICKSf anp L. KNOPOFFT 


ABSTRACT 


The reflection of a time-harmonic Rayleigh wave by a high impedance obstacle in shearless con- 
tact with an elastic half-space of lower impedance is examined theoretically. The potentials are found 
by a function—theoretic solution to dual integral equations. From these potentials, a “reflection 
coefficient” is defined for the surface vertical displacement in the Rayleigh wave. Results show that 
the reflected wave is 7/2 radians out of phase with the incident wave for arbitrary Poisson’s ratio. 
The modulus of the “reflection coefficient”’ depends upon Poisson’s ratio, and is evaluated as rg =0.265 
for ¢=0.25. 


INTRODUCTION 

In this paper the reflection of a time-harmonic two-dimensional Rayleigh 
wave by a high impedance obstacle in shearless contact with a half-space of 
lower impedance is examined. A theoretical treatment is made, and a “reflection 
coefficient” calculated. The mathematical boundary value problem is solvable by 
first expressing it in terms of dual integral equations and then applying the theory 
of functions to solve these integral equations. 

The results show that if the incident Rayleigh wave hasa variation cos w(t—«/y) 
in vertical particle displacement on the surface, then the vertical particle displace- 
ment in the reflected wave will be A sin w(t+x/y). This is a phase shift of r/2 
radians with respect to the incident wave. 


BOUNDARY CONDITIONS 


Consider the geometry of Figure 1, in which a rigid quarter-space is coupled 
to an elastic half-space by means of a thin film of inviscid fluid. Let u(x, y)e~' 
be the elastic displacement vector, and 7;;(x, y)e~“* be the components of the 
stress tensor, where u(x, y) is a solution to 


BV XV Kut y) = 0 (1) 


Tij ‘pB? = (a?/B? = 2)6;;V-u + Ou; + Ou; / Ox; (2) 
u = (uz, ty, 0). (3) 


Here @ and 6 are the P and S wave velocities in the homogeneous, isotropic, 
perfectly elastic half-space of Figure 1. 


* Publication No. 180, Institute of Geophysics, University of California. Manuscript received 
by the Editor April 23, 1960. 
Tt Institute of Geophysics, University of California, Los Angeles. 
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P (x,y) 
P (r,8@) 


elastic half-space 
g 


rigid, incompressible 
Rayleigh wave quarter - space 


Fic. 1. Schematic diagram of the scattering geometry. 


The boundary conditions are accordingly, 
4,0) 
Tal <2 0 
u,(x, 0) 
THE INTEGRAL EQUATIONS 
If the substitutions 
uz = + dp/dy, uy, = d¢/dy — dp/dx 


(5) 


are introduced into equation (1), then the mathematical problem reduces to that 


of integrating the two scalar Helmholtz equations 
(Vv? + ky") (x, y) = 0, w/a 
(V2? + ko*)¥(x, y) = 0, = w/8 


subject to the conditions of equation (4). 
The solutions to (6) and (7) are written in the form 


1 
(x, vy) = do(x, vy) + —{ 


1 
V(x, ») = vols») + f 


where L is the path shown in Figure 2, and 
m = (k? — nz = (k? — 
Re (m1, m2) 2 0. 
The potential @9 and Wo specify the incident wave. We choose 
oo = a — ky*)!/? 


Yo = b (ks? — 


(6) 
(7) 
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k- plane 


Fic. 2. The path of integration L in the k-plane. 


where k;=w/7 is the positive real root of the Rayleigh equation (y is the Rayleigh 
wave velocity) for positive w, 


R(k) = (2k? — ko”)? — 4h*qine = 0 
and the constant Co is given by 


_ 


These Rayleigh wave potentials satisfy the conditions of zero normal and shear 
stress on y=0O for all x. 

By a combination of equations (8)—(12), and (1)—(3), one obtains the set of 
integral equations, 


1 
{ i2kn1A(k) + (ko? — 2k?) B(k)}e*dk =0, <x (13) 


1 
| (ko? — 2k*) A(k) — i2kn2B(k)} e**dk = 0, <x<0, (14) 
WL 


1 he? 
{1A (k) + ikB(k)} = — i— citer, O<x<0. (15) 


Since A(k) and B(k) are not identically zero on the entire path L, then from 
equation (13) we require that 


(16) 


2 
1197 
“ky -k2 ~k, 
| 
‘ 
ke? — 2k;® 2ksa 
= 
! 
Pa 
B(k) = i——— A(R). 
2k? — k 2? 
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The combination of equations (14)—(16) gives the dual integral equations, 
) ‘ 
= 0, <x<0, (17) 
— 2k? 
10) 
— e*dk + 1 — (18) 
2ri — 2k? 
SOLUTIONS TO THE DUAL INTEGRAL EQUATIONS 


To solve equations (17) and (18) for A(R), we use the integral representation 


——dk=ctv, 2£>0, (19) 
k — ks 


where V,(&) is any algebraic function having all of its singularities in the lower 
half-plane Im (&) <0 (or on the negative real axis in the degenerate case of no dis- 
sipation in the elastic medium). The integral in equation (19) thus supplies an 
analytic continuation of the left-hand side of equation (18) for negative values 
of x. Combination of equations (18) and (19) leads to 


A(k) + i = 0, (20) 
anid — 28? — 

<x < Ow, 


Since A(k) and N,(k&) are not zero everywhere on JL, it follows that 


(2k? — ko?) Ni(k) 


(k) = 


Let us combine equation (21) with equation (17) to obtain 
1 Rik) N, (k) 
— e**dk = 0, —-x <xr<0. (22) 
We define a new function K(k) by 
— ke’ *)(R — hs’) 


K(k) = 
R(k) 


It is clear from (23) and from the definition of R(&) that 
K(k) 1, |k|— 


It is shown in the Appendix that K(&) can be factored into a product K,(k)- K(k) 
where K(k) has all of its singularities on the negative real axis, while K_(&) has 
all of its singularities on the positive real axis. From the Appendix 
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= exp — f 


Ki(k)(Rk + 
k+ 


M.(k) = 


it is clear that a general solution to equation (22) is 
Ni(k) = My(k)P(R), 
where P(k) is a polynomial in & such that P(ks) does not vanish. 

The polynomial P(&) is fixed by the requirement that the vertical displace- 
ment u,(x, y) be continuous everywhere on y=0. The only polynomial which 
satisfies this continuity condition is a constant. Hence 

(2k? — ko?) K(k) 
2k(k + 


B(k) = — - (26) 
(k + — k3’)K 


(25) 


The appropriate solutions to the differential equations and boundary conditions 
are 


(27) 


1 i(2k? — ky *) K(k) 


(x, y) = do(x, y) + — 
(x,y) = y) + (k — (kh? — 2y 


2k(k + 
te (k)e*=- nwdk. (28) 


THE REFLECTED RAYLEIGH WAVE 


Consider the evaluation of the integrals in equations (27) and (28) when «<0. 
Jordan’s lemma is satisfied by the integrands on infinite arcs in Im (&) <0. There- 
fore the integral over the path Z can be deformed into the integral around the circle 
I drawn about t!ie simple pole at k= —&;, plus the branch line integral Ty, taken 
about the L-shaped cuts (Figure 3). By the residue theorem the integrals over T 
are readily evaluated and we obtain 


+. (branch line integral over (29 
¥=Yot a) 4. (b h line int | (30) 
= —e ranch line integrai over 


The branch line integrals are not written explicitly since we shall not be concerned 
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;  k-plone 


Fic. 3. The deformed path of integration for x<0, w>0. 


with them, but shall concentrate on the two integrated terms in equations (29) 
and (30). It should be remarked, however, that the branch line integrals in equa- 
tions (29) and (30) represent scattered P, S, and P-S conversion body waves due 
to the corner at (0, 0). These integrals will be interpreted in a forthcoming treat- 
ment of the general diffraction of Rayleigh pulses by structures such as we are 
considering here. 

The vertical component of particle displacement u,(x, 0) is a quantity readily 
measured. From equation (5) the part of u,(x, 0) contributed by the incident at 
reflected Rayleigh waves alone, say u,(x, 0) is 


uy®)(x,0) = — (aCo + tk3)e** 
— K_(k;) 
+ ky 


— ik 3x 


i(aCo ik3) ( 


Let us define a “reflection coefficient” for u,“*(x, 0) as the ratio of the coef- 
ficient of that of Then 


K_(k3) 
ky + hy K 


where A is a real quantity. This solution is the result obtained for w>0. For w<0 
a similar analysis yields 

re = 
For a Poisson’s ratio of 0.25(k2?/ki?=3) we find 

A = 0.265. 

Consider the inverse Fourier transform of the quantity 
i w>0 
w <0. 
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Fic. 4. Paths of integration used in the decomposition of K(). 


This function is —2/zt (Erdelyi, et al, 2.2-3). Hence in order to specify the wave 
shape of the reflected wave, in the case when the incident wave is an impulse, 
let the time dependence of the incident wave be 
1G x/y). 
Then the reflected wave has a time dependence and amplitude given by 
2 a 


x 


An elementary application of the convolution theorem yields the wave shape 
when the shape of the incident pulse is not that of a delta function. 


APPENDIX 


Any function of a complex variable analytic in a horizontal strip can be de- 
composed into a product of two functions each of which is analytic in an upper 
and lower half-plane; the region of overlap is the original strip of analyticity. 
By Cauchy’s theorem, we may write 

1 log K(z) 
log K(k) = log K,(k) + log K_(k) = — ——— ds 

p 2 — 
where K(k) is the function defined by equation (23), and P is the path of Figure 4. 
The integral around P is equivalent to an integral taken along the real axis and 
closed by an arc at infinity in the lower half-plane. This integral has a contribu- 

tion only from the branch line integral B_. Thus we may write 
log K,(k) = dz. 
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A similar construction in the upper half-plane leads to 


1 
log Ki(k) = —{ 


where we use the Principal Branch of the logarithmic function. We obtain, after 
a calculation involving the definition of K(s), 
dw 


1 ke 4w?(w? k 2 1/2( Ro? w?)i/2 
ky (2w? — wtk 


ke 4w?(w? — — 
K.(k) = exp tan7! 
(2w* — 


as the required decomposition. 
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TRANSMISSION AND REFLECTION OF RAYLEIGH 
WAVES BY WEDGES* 


L. KNOPOFFf anp A. F. GANGIT 


ABSTRACT 

Experimental observations have been made of the transmission and reflection of Rayleigh waves 
by wedges. Results are reported for Rayleigh waves in aluminum wedges. It is observed that the 
wave shapes of the transmitted and reflected waves differ from that of the incident wave and depend 
on the angle of the wedge. The change of shape is attributed to an interference between part of the 


incident wave-form and the radiation from a line source placed at the vertex. A procedure is given for 
the calculation of the partition between the two terms. 


INTRODUCTION 


The problem of the reflection and transmission of Rayleigh waves by wedges 
is important. It will be shown in a later paper that this process is significant in 
the production of non-coherent seismic ‘‘noise.” The problem is difficult to solve 
mathematically because body P and S waves are also generated by a Rayleigh 
wave incident on the corner of a wedge. However the problem can be modeled 
quite simply on seismic laboratory models. In this paper we present results of 
model experiments and an interpretation of the data. 

Consider a homogeneous, isotropic, elastic medium filling the region 0<6<4, 
where (r, 6, z) isa cylindrical coordinate system centered 
on the corner of the wedge. The remainder is vacuum. The two wedge faces are 
given by the coordinate planes @=0 and @=4p (Figure 1). 

Physically one would expect that a Rayleigh wave incident on a corner would 
excite the corner as a virtual source. By Huygens’ principle, Rayleigh waves 
should be excited on both faces of the wedge and P and S body waves radiated 
into the elastic medium. In this paper only the transmitted and reflected Rayleigh 
waves will be considered. 

The simplest theory states that simple reflection and transmission coefficients 
suffice to describe the Rayleigh waves on the two faces of the wedge. These co- 
efficients are pure numbers, depending only on the angle of the wedge and the 
properties of the elastic medium. The radiated body waves, of course, have more 
complicated distribution functions; their amplitudes depend not only on the 
wedge angle and elastic properties but also upon the angular coordinate 6. A de- 
scription of reflection and transmission of Rayleigh waves by such a process is 
given by deBremaecker (1958). 


* Publication No. 179, Institute of Geophysics, University of California. Received by the Editor 
April 7, 1960. 

t Institute of Geophysics, University of California. 

t Space Electronics Corporation, Glendale, California. 
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Fic. 1. Schematic diagram of the wedge experiment. 


Experimental measurements 


Measurements of the incident reflected and transmitted Rayleigh waves were 
made in the laboratory using standard model techniques. Aluminum sheets were 
cut into wedge shapes (Figure 1). Measurements were taken for wedge angles 45 
between 45° and 135° in 5° steps and in 10° steps for wedge angles between 230° 
and 300°. 

Motions in thin plates have been found to be a good approximation to two- 
dimensional elastic wave motions where the “wave lengths” of the pulses are 
large compared to the thickness of the plates. By the wave length of a pulse we 
mean the wave length associated with the frequency of the maximum of the 
Fourier spectrum of that pulse. 

A barium titanate source transducer was glued onto the edge of the sheet along 
one of the wedge surfaces and excited to give a short stress pulse in the plate. The 
observed incident pulse had three extrema; the time interval between the first 
and last extrema was about 3.4 usec. The properties of the aluminum are: plate 
P wave velocity 5.76 mm/ysec, S wave velocity 3.26 mm/usec and a plate 
Poisson’s ratio of 0.266. The seismic wave length for P waves is thus about 2 cm. 

Rayleigh waves in two-dimensional structures should propagate unattenuated 
in perfectly elastic materials. Only because of non-elastic effects should there 
be any attenuation of the surface wave amplitude with distance. The reflected 
and/or transmitted Rayleigh waves were observed at a number of stations and 
the attenuation in aluminum was found to be negligible for the accuracy of the 
experiment. The change in amplitude due to attenuation was less than the scatter 
due to the irregular placement of the detector and interference from reflected 
body waves. Measurements were taken at a number of stations for each angle and 
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the results averaged to avoid experimental inaccuracies resulting from detector 
placement error. 

Typical observations of the normal component of motion of Rayleigh waves 
on aluminum plates cut in wedge shapes are shown in Figure 2. The details of 
the measurements are given later. Comparison of the reflected and transmitted 
Rayleigh waves with the incident Rayleigh wave shows that there is a change in 
the wave-form of these waves. This is especially clear in the case of the trans- 
mitted wave for the 80° wedge. It is apparent that a more complicated interpre- 
tation of the data is necessary than that offered above; the simple theory is not 
sufficient to fit the observations. 


Inter pretation 


The basis for the interpretation comes by noting that, in the case of the 
incidence of a simple pulse upon the degenerate wedge @=2z7, i.e. a half-plane 
discontinuity, the scattered wave does not have the same time dependence. 
In particular let an incident plane pulse be represented by K6(¢—x/c) where K 
is some amplitude and 6(¢—.x/c) is the delta function. Then the wave scattered 
by a two-dimensional half-plane has the form K’(t—r/c)~?H(t—r/c) where 
H(t—r/c) is the unit step function, r is the distance from the edge of the discon- 
tinuity and K’ is a function of the angle (deHoop, 1958). This description is satis- 
factory as long as we consider motions close to the onset time r/c. In other words, 
the edge of the discontinuity acts as a two-dimensional line source in an infinite 
medium. 

On the other hand, for the degenerate wedge having wedge angle z, the trans- 
mission coefficient for Rayleigh waves is unity and the reflection coefficient is 
zero. In more mathematical terms, for an incident Rayleigh pulse 6(¢—x/c) the 
transmitted Rayleigh pulse is also 6(t—x/c). 

Hence, we represent the transmission and reflection of Rayleigh waves as the 
sum of two terms. The first of these corresponds to the normal transmission and 
reflection with unchanged wave-forms. The second corresponds to the Rayleigh 
waves radiated from a line source placed at the vertex of the wedge. It can be 
shown that, if the line source is excited as a delta function, the Rayleigh waves 
from this line source behave as A(ct/x—1)~'?(ct/x—1) at a great distance from 
the vertex where x is a coordinate along the surface. The amplitude A is an un- 
known function of the wedge angle. The wave shape from the line source is as- 
sumed to be independent of the wedge angle. 

The Fourier transform of the radiation from a line source at the vertex of a 
wedge is 


a(n, (= r) ei dn +f b(n, (= r) 
0 Up 0 Us 


where a and b are determined from the boundary conditions on the faces of the 
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wedge and v, and », are the P and S wave velocities. Passing to the limit of large 
argument wr/v>>1 by the use of the asymptotic expansions for the Hankel func- 
tions, the inverse transforms of the integrals that result have a behavior of the 


type 


C(O, 
(ot/r — — 1) 


for times near the arrival time r/v. 

Thus, the virtual source represented by the corner performs two operations 
for an incident surface wave having a delta function time dependence. These are 
a reproduction and a mild integration, the latter corresponding to the term 
(Vet/x—1)~'?H(V et/x—1) where Vz is the velocity of Rayleigh waves. We repre- 
sent the reflected and transmitted waves as the sum 


7) — 1) + — H(Vet/r — 1) (1) 


for times close to the onset time for an incident delta function. The subscripts R 
and T refer to the reflected and transmitted Rayleigh waves. The functions A 
and B are dimensionless coefficients indicating the partition between the two 
processes and a is Poisson’s ratio. 

The use of expressions of the type (1) are valid for times close to the arrivals. 
The errors are not significant if the actual wave-forms of the input pulses are of 
short duration compared with the travel times (Knopoff, 1958). 

The experimental input pulse is not a delta function but instead has a finite 
duration. We approximate the input wave function as the sum of three gaussian 
functions (Figure 3). This is not an inaccurate decomposition. We now find the 
response to the system to excitation by a single gaussian function. If the input 
wave function is exp | —8*(Vrt/r—1)?} then, by the convolution theorem, the re- 
ceived motion for both reflected and transmitted waves will be 


90, 


x 
= Ar,1r(Oo, + Br, of 
0 


Let 7 be the dimensionless arrival time of Rayleigh waves 


= B(Vat/r — 1). 


Fic. 2 (a & b). Typical observations of the incident reflected and transmitted Rayleigh waves for 
wedges of different angles. 
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Fic. 3. Representation of the wavelet of the incident wave as the sum of three gaussian functions. 


The wave function at the receiver is thus, 


x 


0 


Ar,rfo(t) + (4) 


where fo and f_1/2 are defined as the gaussian function and the convolution inte- 
gral respectively. 

The universal function f_1/2 has been evaluated numerically on SWAC! for 
the range r: —3(.1)3 (Table I). Outside this range for large positive or negative r 
standard approximations apply. For instance for large positive r. 


This result can be applied to any gaussian function of arbitrary width by suitably 
selecting the scale factor 8. Thus, for the given input function (Figure 3) the con- 
volution with ((Vet/r)?—1)' is given in Figure 4. The solid curves are the con- 
volutions of the individual gaussian functions and the dotted curve is the sum 
of these three contributions suitably spaced in time. 

For the particular input wave function go(t) given in Figure 3 the reflected 
and transmitted wavelets will be the functions. 


fr.r\t, 9, o} = [go(t) + ar, o)g-1/2(t)], (5) 


1 Standards Western Automatic Computer, an electronic digital computer at University of 
California, Los Angeles. 
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TABLE [ 


f(t) t t fit) 
00000 2067 1.0 .9868 
—2.9 —0.9 2587 1:1 9526 
—2.8 $2 9166 
—3.7 0003 3844 £.3 
0005 —0.6 4564 1.4 8444 
0008 .5328 1.5 
—2.4 0013 —0.4 6116 1.6 
—2.3 0021 6906 4.7 
—2.2 0033 7674 1.8 .7193 
—2.1 0052 1.9 
—2.0 0080 0.0 2.0 
0120 0.1 2.1 6495 
—1.8 0177 0.2 1.0089 2.2 6304 
—1.7 0257 0.3 1.0432 2.3 .6129 
0367 0.4 1.0655 2.4 5970 
0512 0.5 1.0758 25 
a —1.4 0702 0.6 1.0749 2.6 5690 
0946 0.7 1.0640 2:7 5566 
: —1.2 1225 0.8 1.0446 2.8 5450 
1.1 1622 0.9 1.0185 2.9 5343 
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CONVOLUTION SUM 


Fic. 4. Convolution of each of the three gaussian functions with (ct/x—1)~"? and the 
convolution of the incident wavelet. 
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Fic. 5. The sum go(t)+ag_12(t) with a as a parameter. 


where g_1/2 is the function of Figure 4. We have written arg,7= Br,r/Ar,r and 
the time scale is measured such that 8=1. The function f/A has been calculated 
for values of a= ~(A=0), +4, +2, +1, +3, +}, 0. These are shown in Figure 5 
and represent possible output wavelets for the appropriate ratio coefficient a. It 
can be seen that the shape of the pulse changes with a. We note the appearance 
of an additional loop on the right for a<0 and note the disappearance of the 
loop on the left for a< —2. 

Identifying the first negative loop as having amplitude a, and numbering the 
loops from 1 to 4, the ratios of the amplitudes a;/a2, a;/a2 and a;/a_; may be plotted 
as a function of a (Figure 6). This curve is now the basis of the interpretation. 

The procedure for the calculation of a from the data is as follows. Number the 
loops, measure the amplitudes, take the ratios and enter in Figure 6. The ratios 
should all yield an identical value of a. There is some ambiguity in determining 
the first negative loop as the multiplier function Az,r7 may be negative. However, 
any error in number will show up in terms of ratios of amplitudes that do not 
yield a consistent value of a; the two (or three) values of a corresponding to the 
amplitude ratios must all have the same numerical value for a given reflection 
or transmission. 

Once the values of ag,r are known, the calculation of the quantities Ar, Br, 
Ay and Br is elementary. Auxiliary curves of the amplitudes a; as functions of a 
are useful in this last interpretation (Figure 7) where the curves have been nor- 
malized to the values in the incident pulse, a=0.? It should be emphasized that 


2 By comparing the values of a; for the value of a already obtained with the measured values of 
a;, A can be determined. The value of B is obtained from the definition of a. 
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Fic. 6. Relative amplitudes of the loops in the received wavelet as a function of the parameter a. 


Fic. 7. Amplitudes of the loops in the received wavelet as a function of the parameter a. 


5 
4 
4 
‘i 
Gs 
-3 0 2 3 4 


1212 L. KNOPOFF AND A. F. GANGI 


2 T 1 T T T T T T T T T | ad 
tan a “4 


2 


2,.-! e— reflection - 
tan ay o — transmission 
adi. 1 1 1 1 L 1 1 1 4 
ce) 40 80 120 160 200 240 280 320 360 


WEDGE ANGLE (degrees) 


Fic. 8. Experimental values for a(@). 


the values of az, a7, Ar, Ar, Bre and Br thus calculated hold for any shape of 
pulse as long as the pulse duration is short compared with the travel time. 


Reflection and transmission coefficients 


By the method described in the preceding section the results for ag, 7r(o, 0.266) 
are presented in Figure 8. The curves of a,/az are very flat for large positive or 
negative a and the wave-form does not change much with large changes of @ in 
this region, hence an inaccuracy in @ may be reported that may be large but yet 
is of little significance. In order to avoid this complication we have compressed 
the scale —x<a<x to —1<2/mr tan™a<1 and plotted the function 
tan 0.266). 

The results in Figure 8 together with the values of the reflection and trans- 
mission coefficients Ar, Br, Ar, Br (Figure 9) form the complete interpretation. 

It should be noted that the interpretation curves for Figures 6 and 7 are not 
universal. They depend severely upon the shape of the input pulse (Figure 3). 
This condition requires that the input pulse be kept constant for the use of the 
interpretation curves. Mild variation in the experimental input wave-form was 
the major source of the inaccuracy in the interpretation. Interpretation curves 
must be recomputed for each input function. 


Energies 


deBremaecker (1958) has reported results of reflection and transmission of 
Rayleigh waves at corners for incident Rayleigh waves. The analysis therein pre- 
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Fic. 9. Experimental values for Az, Br, Ar, Br. 


sented takes no account of the changes of wave shape upon reflection. The quality 
of the seismic pulse in the experiments is poor. Comparison with deBremaecker’s 
work may be made for the energy reflection and transmission coefficients. 
Consider a gaussian wavelet Ao exp —(t/lo)? where fo is the breadth of the 
wavelet and Ao is the peak amplitude of the wavelet. The energy in the wavelet is 


f Ao” exp (—20?/to?)dt = (6) 
We thus compute the energy in a wavelet with N loops as 


N 
>> 
(7) 


where A; is the amplitude of the peak of the ith loop and Af; is the width of the 
loop, taken here to be the time interval between zeros of the ith loop. This may 
be normalized to the input wavelet. The energy in (7) is an approximation neg- 
lecting the ‘‘overlap” or interaction of the individual gaussian functions. 

Energy reflection and transmission coefficients for the wedges are presented 
in Figure 10 and compared with deBremaecker’s results. We note that although 
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Fic. 10. Energy reflection and transmission coefficients. Results for aluminum (¢=0.266) are 
compared with those obtained by deBremaecker (¢ =0.17). 


the present data have been taken at an angle interval of 5°, no significant differ- 
ences are observed in the structure of the curves obtained by deBremaecker for 
an angle interval of 10°. Significant shifts in the peaks both in angle and in mag- 
nitude, are observed. We note that the Poisson’s ratios for the two experiments 


are different. 
The reflected energy for wedge angles greater than 180° is less than 1 percent 


of the incident energy. 
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ON A SUCCESSIVE APPROXIMATION METHOD FOR 
INTERPRETING GRAVITY ANOMALIES* 


Z. FRANKENBERGER DANESt{ 


ABSTRACT 


A new method for quantitative interpretation of gravity anomalies is presented. The disturbing 
body is represented by a finite number of vertical prisms arranged on a pre-determined, regular 
grid. The horizontal dimensions of the individual prisms are small enough that they can be approxi- 
mated by vertical-line mass elements at the axis of t he prisms. Formulas for gravity due to one prism 
are derived and, for the case of Gulf Coast salt densities, plotted on Figure 1. Gravity due to the 
whole body is an algebraic sum of the contributions of all ‘prisms at the appropriate depths and dis- 
tances. 

This method makes possible the direct interpretation by successive approximations by introduc- 
ing proper geologic limitations. All the numerical work can be conveniently done on a high-speed 
digital computer. The method is especially suitable for features with predominant vertical dimension 
such as salt domes and igneous plugs. It gives at least the same, and possibly higher, degree of ac- 
curacy as the graphical dot chart methods and, carried out on a digital computer, should be about 
two thousand times faster. 


SYMBOLS USED 


a= theoretical logarithmic density contrast of salt against sediments at the 
depth of one ft in the Gulf Coast 

b=theoretical increase in logarithmic density of the Gulf Coast sediments 
between 1-ft and 10-ft depths 

d=distance of axes of two adjacent prisms; spacing unit 

g=observed gravity anomaly due to the mass M 

gravity of a vertical line element 

m=mass at the point (x, y, 3) 

n=number of successive approximations 


r= V/(x—x0)*+(y— yo)? +2 

x=horizontal coordinate, East, in feet 

y=horizontal coordinate, North, in feet 

2=depth, in feet 
depth to bottom of prism 
27 = depth to top of prism 

A=cluster of prisms approximating the disturbing mass M 

’ = F(x, y, 3) =0=formula of the surface of the anomalous mass M 


a da 
G= lim ( -f : ) See equation 6 
rT r 


H=half width of observed anomaly 
H = half width of anomaly of a line element 


* Manuscript received by the Editor April 6, 1960. 
t Boeing Airplane Co., Seattle, Washington. 
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L=vertical line element 
M =anomalous subsurface mass 

N=number of prisms approximating M 

P=reference point, gravity station 

R=residual gravity map; radius of cylindrical shell 

S= horizontal cross sectional area of an elementary prism 

V=volume of the anomalous mass M 

= gravity constant = 310~’ c.g.s.u. 

6=g—g=difference between gravity of M approximated by line elements 

and actual gravity 
e= the ratio between the tolerable discrepancy and the total gravity anomaly 


1; 


—— ; See equation 20 
8 2°+9/8d? 
p=V(x— x0)? yo)? 

po= distance at which 6=0 

o=density contrast a point (x, y, 3) 

a =closed surface around M 


INTRODUCTION 


A mass anomaly M causes at a point P an anomaly g in the vertical component 
of the earth’s gravitational field: 


f yzdm f zodV 
= - =¥ 
Vv 


M 
ozdxdydz 
[(x — 40)? + (y — yo)? + 


The constriction on the limits is given by 
F(x, y, 2) = 0. (2) 


To perform the integration in equation (1) is a difficult, and often impossible, 
task even for simple F and c. Therefore, graphical, numerical, and analog methods 
have been developed. The first and second methods give satisfactory results but 
are too time consuming. Analog methods are fast, but the accuracy is limited. 

With the development of fast digital computers, numerical methods become 
practicable and may even replace the trial-and-error way to define F. 

To present one way to the solution of this problem is the purpose of this paper. 


THEORY 


The vertical component of gravity, dg, at a point, P (2), yo,0) due to a mass 


‘ 
mish 
: 
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element dm, is 


Integrating over the whole mass M gives equation (1) where the limits satisfy 
equation (2). 
If the density contrast is a function of depth only, 


= a(z), (3) 


and if the horizontal dimensions are constant and negligible compared to the dis- 


tance from P, approximately 
sf dz. (4) 


4 22 dg 


In all practical cases, o can be approximated by.a function of z which permits 
integration. However, it is an advantage to use such a function for which the 


limit 
da 
G = tim - f ) (6) 
r 


exists. Then all prisms can be assumed to extend downward to infinity and graphs 
or tabulations can be made for zr only; and, from (5) and (6) 


da 
r 


Gravity of a finite prism is then obtained as a difference of gravity of two semi- 
infinite prisms with the depths at the top and bottom of the finite prism, 
respectively. 

In the Gulf Coast, where the density contrast of the salt of the domes against 
the surrounding sediments is given (approximately) by the formula 


o = at b logi 2, (8) 


Integrating by parts, 


formula (7) gives 


yd’ logio 27 
4 = 


+ byd?| = 
V(x — Xo)? + (y Yo)? + V/(x — xo)? + (y = yo)? + (9) 


Vv (x— x0)? +(y— yo)? 


‘ 
4 
: 
ne : 
an 
t 
° 
an 
(7) 
ax 
\ 
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This formula (9) is plotted on Figure 1 for a= 1.023, b= —0.2897, d= 2,000 ft, and 
sr =500 ft; 1,500 ft; 2,000 ft; 3,000 ft; 7,000 ft; 10,000 ft; 15,000 ft; 20,000 
ft; 30,000 ft; and 40,000 ft respectively. 

An obvious advantage of the logarithmic formula (8) is in the possibility of 
subsequent density adjustments. Both parameters a and 6 appear as multiplica- 
tive factors, so that if each bracket is evaluated separately, density adjustments 
involve only multiplications and additions of two terms for each station. 

Next let us check how closely g is approximated by g. To that end, let us 
follow the subsequent reasoning. 

First, let us assume that V=1 and M>O. (For negative density contrasts, 
such as those of deep salt domes, where M <0, the changes are obvious.) 

From Gauss’ Theorem it follows that 


ff = 4ryM. (10) 
Il 
ff (11) 
IL 
ll 


For a point at the axis only, all dg’s are added algebraically, whereas dg’s are 
added vectorially; therefore 6, must reach its maximum value, dmax, at the axis. If 
bmax does not exceed the limits of tolerable errors, other 6’s, smaller than dmax, 
certainly do not either. Let us then carry out our proof for dmx. 

Somewhere off the axis, at the distance po, 6=0 and beyond that 6<0, reaching 
a minimum. If x extends to infinity, 6 converges to zero asymptotically. 

Let us first assume a semi-infinite cylindrical shell of radius R and constant 
density contrast. For a point at its axis, equation (1) yields 


otRdadz 2rycR 
ano (8? + + R*)"? 


From equation (5), 


But also 


so that 


2ryaR 
(14) 
thus 
2ryoR 1 
§6= g (15) 
2T 
i+— 
Zr” 
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Postulating that 


< €g 


gives approximately, for sufficiently small ¢, 


le (16) 


as the upper limit of R. The condition is certainly satisfied by another body with 
the same vertical axis and entirely within R. If, therefore, the diagonal of a square 
prism is equal to 2R, 


< 


and 


«> (17) 


4 


Table I shows the maximum possible errors introduced by this method. 


TABLE I 
MaximuM PossIBLE ERROR AS A FUNCTION OF RATIO 27/d 


(Actually, the situation is much more favorable. In the example on Figure 5, 
very satisfactory results are obtained with z7=d). 

It remains to show that po is actually less than the prism spacing , so that the 
discrepancies of the adjacent prisms do tend to cancel one another. To show that, 
let us consider a square prism of cross-sectional area S=d*. If approximated by 
a vertical line element at its axis, the gravity of this prism at point P at the axis 
of the adjacent prism (at distance d) is 


(18) 


Now let us divide the same prism into four equal prisms of S’ = jd* (Figure 4), 
replace each of them by a vertical line element at its particular axis, and calculate 
their combined effect at P: 


1 1 1 
9 ( 1 /1— 4 ) 


3 
i 
: 
a 
€% 
2 6.2 
3 2.8 
4 1.6 
10 0.2 
Ve 
1 
V2 + 
one 
> 
HENS 
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where 


(20) 


——— . (21) 


Obviously, 2’ gives a better fit to the actual g than g does. If, at P, g<g’ the 
gravity discrepancies will tend to cancel each other and 6<émax. Therefore, let 
us prove that at P, <2’: From equations (19) and (21) 


V1 — V1 + 


Postulating the 2’ = g requires that 


+ 8&* + 20& + = 0. (23) 


Equation (23) has only two real roots: &=0, which requires that d=0 and is there- 
fore trivial; and + —0.05, which has no physical meaning, as it requires com- 
plex d or z. Therefore, 2/2’ cannot become unity for any real finite d, and if 
2< in one case, it has to be soin all cases. Now for z=0 and d= 1, equation (20) 
and (22) yield 


1 
(24) 
and 
— = 1.024621 ---, (25) 
therefore, 
(26) 


quod erat demonstrandum. 

Attention shall now be paid to the placement of prisms as suggested in the 
preceding chapter, step d. Table II and Figure 3 show the sequence of the prisms 
according to the increasing distance, p, the number of prisms in that ring, and 


| 

1 d? 
— 
— 

Rewriting equation (18) gives 

8 

7 
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the number of prisms, N, in and inside that ring. Figure 3 shows a little more 
than the first quadrant of the placement of the prisms. 


METHOD 


The method is based on the assumption that the anomalous mass M can be 
imagined as composed of a finite number, N, of vertical prisms. The horizontal 
dimensions of each prism are small enough, so that within the limits of accuracy 
the gravity of one prism can be approximated by that of a vertical mass line ele- 


TABLE II 
SPACING OF PRISMS ON A SQUARE GRID 


Square of 


Square of 

Consecu- Number number Consecu- distance Number =number 
_Uve from center of Prisms of prisms in | _ tive from center of Prisms of prisms in 
Number and inside | Number /. and inside 

of Ring (in spacing in Ring ring | of Ring (in spacing in Ring ring 

units) units) 

1 0 1 1 26 52 8 169 

2 1 4 5 27 53 8 177 

3 2 4 9 28 58 8 185 

ft 4 4 13 29 61 8 193 

5 5 8 21 30 64 12 205 

6 8 4 25 31 65 8 213 

7 9 4 29 32 68 8 221 

8 10 8 37 33 72 4 225 

9 13 8 45 34 73 8 233 

10 16 + 49 35 74 8 241 

11 17 8 57 36 80 8 249 

12 18 4 61 37 81 4 253 

13 20 8 69 38 82 8 261 

14 25 12 81 | 39 85 16 277 

15 26 8 89 | 40 89 8 285 

16 29 8 97 41 90 8 293 

17 32 4 101 | 42 97 8 301 

18 34 8 109 43 98 4 305 

19 36 + 113 | 44 100 12 317 

20 37 8 121 45 101 8 325 

21 40 8 129 46 104 8 333 

22 41 8 137 47 106 8 341 

23 45 8 145 48 109 8 349 

24 49 4 149 49 113 8 357 

25 50 12 161 50 116 8 365 

51 117 8 373 


ment at the prism’s axis. It can be expressed as a function of depth and horizontal 
distance by a simple analytical formula. Gravity of the whole mass, M, is then 
computed as an algebraic sum of the contributions of all NV prisms at their proper 
distances and depths. Figures 2 and 5 show the suggested shape and location of 
a salt dome in the South Timbalier Area, offshore Louisiana, divided into 151 
such prisms. The horizontal cross-section of each prism is 2000 by 2000 ft. 

The inverse problem; namely, to find F(x, y, 3) to fit g(x, y, ¢), does not have 
a unique solution. However, limitations can be imposed, so that the solution 
does become unique or at least constrained to a certain range. Such limitations 
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may be, for example, a given depth to the mother salt of a salt dome and a postu- 
late that F(x, y, z) be of a certain type; or a given depth to the top (as revealed 
by drilling, for example) and other restrictions on F(x, y, 2). 

In the following, a method of successive approximations for deriving the shape 
of a salt dome is described under the following assumptions and postulates: 


1. The bottom of the dome is flat at the depth, zz. 
2. F(x, y, 3)=0, which determines the upper surface of the dome, is a single- 
valued function of x and y (no salt overhangs). 
. Density contrast between salt and sediments is a logarithmic function of 
depth, c=a+b logio z, independent of x and y. 


The solution then consists of the following steps: 

(a-1) Find the point, P;, on the residual gravity map, R, with the highest nega- 
tive value, g:. 

(b-1) Determine the average half width, H,, of the anomaly around P;. This may 
be done on a digital computer; so that values on the map in four directions; 
say north, east, south, and west; are compared with 3g:. If exactly this value 
is found, its distance from P, is used with the weight }. If not, the distances 
of the nearest two adjacent values, one larger than 3g, and the other 
smaller, are used with the weight } each. 

Find a vertical line element ZL; with the bottom at zz, giving approximately 
the same half width, 7/;. (If H; is larger than the largest H,, use this maxi- 
mum value of H;.) Its top is at 71. The gravity value at the axis of this 
line element is £1. 

N=g,/Z, is the number of elements L; necessary to account for g:. There- 
fore, cluster V, elements L; at and around P,, according to the pattern in 
Figure 3. If \; is not in the figure, use the next lower number. This cluster 
represents the first approximation, A;, to the function, F. 

Evaluate the gravity of A; for all points, P, and subtract from the residual 
values. The map of the differences, Re, is the second residual map and is a 
measure of the approximation, A;. If A; is satisfactory, the process is 
finished. If not, the operation is repeated on Ro. 

Find the point, P2, on map, Re, with the highest negative value, go. 
Determine the average half width, H2, of the anomaly around P, in a way 
analogous to that described in b-1. 

Find a line element, L2, with the same or almost the same half width 77, 
in a way analogous to that in c-1. However, the bottom of Lz may be at 
the top of another prism rather than at 2, if such a prism was required 
by previous approximation. The top of L2 is at zr2. The gravity value at 
the axis of Le is 

N2= g2/Z2 is the number of elements, L2, necessary to account for go. There- 
fore, cluster V2 elements, L2, at and around P, in a way analogous to that 
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(441 45] 46] 
41 | 39736 | 37 | 38739441 | 42] 
(8) | (8) | (4) | (6) | (48) | 


i _| 30| 31 | 524-34 | 36 [44 | 49 
(8) | | (OB) | (BAL (@N 


28] 251 24] 25 | 27| 26430 N35 | 39 | 43] 49 
(8) | (8) U2) | (4) | (12) (8) | 02) (8) 
-23| 21 | 20| 20] 21 | 23+ 26 | 29 | 337] 39 | 44 
(8) | (84-6) (4) | (8) | (8) 
16 | 15 | 15 | 16 | 18) 22 29 | 55] 47 
(8) | (8)-4-(8) | (12) | | (8) (B) | (8) | (8) 
147 13] 11 | 13 22 Ne 42] 50 
(8) | (8) (8) | (12) 14) | (8) | 12) (8) 
7 | 8 | 9N 121] 14 | 18 | 23) 28) 34] 41) 48 
(8) | (8) | 4) | (8) | (8) |  @\ (8) | \@ 
/9 Ne | 13 |\27 39 | 
(8) | (4.1.8) | (4) | LA4) | | | (8) | | 
2|¥ 8 | 15 | 20] 25] 31 sel} 4s 
(8) | (8) | (4) 44) (8) | (8) (8) | (8) | | (8) | 
a 7 10 | 14 19 | 24 > 44 | 
5 5 [8] 15] 20] 25] 31 45 
(8) | (8) | | (8) | (8) | AB) | (8) | | (12) | (8) | 
Ks | 4 | 9 21 | 27| 39] 46 
\ 8) | | (84 (8) (8) | | /@)/| 16) (8) 
9 7 14/| 18 A 23 28 | 34] 41 


Fic. 3. Spacing of prisms on a square grid, upper number=consecutive number of ring, lower 
number in parenthesis =number of prisms in the particular ring. 


in d-1 but piling them on top of elements previously postulated where 
necessary. The combination of A; and the new N»2 elements represents the 
second approximation, A», to the function, F. 

(e-2) Evaluate the gravity of A, for all points, P, of the residual map and sub- 
tract from the original residual values. The map of the values thus derived, 
R;, is the third residual map and is a measure of the approximation of A». 


This operation can be iterated until a satisfactory approximation is obtained. 
An analogous procedure can subsequently be applied on gravity maxima, and 
pieces of the prisms can be removed. Alternatively, both ways can be used if 
attention is focused to the point with the highest absolute value of gravity, | g.| , 
adding or subtracting the prisms according to the sign of gi. 

In the ideal case, after m steps, the residual map R,41 would have zero values 
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a —4 


Fic. 4. Gravity effect of a square prism at the axis of an adjacent prism. 


at all points P. As this case will hardly ever be achieved, in practice the procedure 
has to be stopped when a certain degree of approximation is obtained. The cri- 
terion depends on the nature of the problem. 


APPLICATION 


Part of the method has already been programmed for the IBM-704 computer 
by Dr. E. Weinberger at the Gulf Research Laboratory, so that gravity effects 
of a salt dome of a known or assumed shape and density contrast can be deter- 
mined. Dr. Weinberger also evaluated equation (9) in a tabular form. The method 
then was tested on the above mentioned dome shown on Figures 2 and 5. Three 
independent evaluations were performed and are reported in Table III: 

1. Evaluation by means of conventional dot charts. 

2. Evaluation by means of the table evaluating equation (9). 

3. Evaluation by the digital computer. 

The same parameters were used in all three evaluations. Results, in general, 
show good agreement among all three methods. Differences average a few 
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| 
— +20,000 — 


LAMBERT COORDINATES 
STATE OF LOVISIANA 


| 


2,360,000— 


| 25} 


Assumed Structural Contours on Salt 
Depth in Kilofeet 
Top View of an Elementary Prism 
Depth to Top in Kilofeet 
_ Gravity Evaluation Point 
~. Gravity in milligals from Prism Analysis 
Gravity in milligals from Dot Chart Analysis 


Over the body of the assumed salt dome, gravity was evaluated for the point at the axis of the particular prism. Bottom 


of all prisms assumed at 30,000’. Depth contour —17.5 identical with — 22.5. 


Fic. 5. Gravity of assumed shape of salt dome. South Timbalier area, Offshore Louis. 
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TABLE III 


COMPARISON OF RESULTS OF THE PrisM METHODS AND Dot CHART METHOD 
Salt Dome Anomaly 8-6-G, South Timbalier Area, Louisiana 


Gravity Difference 
Calculated Gravity (milligal) 
La. Coordinates of (milligal) —_—_——_—_— 
Prism Prism Prism 
Prism Man. Comp. Comp. 
Chart Manually Computer Minus Minus Minus 
Dot 
Chart 


—0.07 
—0.03 
—0.05 
—0.01 
+0.01 
—0.04 
—0.07 
+0.03 
+0.14 
+0.08 
+0.14 
—0.03 
+0.04 
—0.03 
+0.03 
+0.01 
+0.01 
—0.07 
+0.03 
—0.01 
+0.02 
—0.02 


| 


| 
= 
= 


| 
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2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
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* Error Suspected. 


hundreds of a milligal and in most cases can be ascribed to errors from round- 
ing off. Differences reaching or exceeding a tenth of a milligal, found in seven 
cases, are obviously due to occasional errors, as only one evaluation differs, while 
the other two are in good agreement. Of these errors, the dot chart method seems 
to involve three; the prism method manually, four; the computer, none. 
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HEAT FLOW IN PENNSYLVANIA AND WEST VIRGINIA* 


WILLIAM B. JOYNERT 


ABSTRACT 


In order to obtain heat-flow values for six wells in Pennsylvania and West Virginia a simple 
technique was developed for estimating thermal resistivities from well-sample logs. This technique 
made use of average resistivity values for various categories of sedimentary rocks. The averages were 
calculated from available resistivity measurements on rock specimens, many of which did not come 
from the region in which the wells were located. The results indicate that, allowing for error in the 
estimates, the heat flow at the wells lies between 1.1 and 1.5 microcal/cm? sec. The problem of de- 
veloping a more accurate technique for estimating resistivity is discussed. Such a technique would be 
valuable in extending our knowledge of the earth’s heat flow. 


INTRODUCTION 


A major source of information about the interior of the earth is measurement 
of the rate at which heat flows out to the surface. A knowledge of heat flow is also 
important in seeking the causes of crustal deformation. To determine this 
quantity it is necessary to know both the vertical temperature gradient and the 
thermal resistivity (the reciprocal of the conductivity) at the same place. Un- 
fortunately, there are few places where both these requirements have been met. 
In 1954 Birch reviewed the status of geothermal investigations and published a 
map showing the distribution of reliable heat-flow determinations. Since that 
time much more work has been done at sea and at least three new measurements 
have been made on land, but the data are still meager. 

Temperature measurements made by Van Orstrand et al of the United States 
Geological Survey are available for hundreds of wells in this country. There are 
no cores from these wells for determinations of thermal resistivity, but detailed 
sample logs have been made for some of them. The present study, originally done 
in 1953, is an attempt to estimate the resistivity from the sample logs and thereby 
obtain heat-flow values for six wells in Pennsylvania and West Virginia. It is 
also hoped that this work will contribute to the development of a general tech- 
nique for estimating the resistivity of sedimentary rocks. Such a technique would 
be of considerable help in geothermal studies, since direct measurements of re- 
sistivity are possible only where suffrcient cores are available. 

The wells are listed in Table 1, and Figure 1 shows their locations. Both 
temperature measurements and sample logs were available for wells no. 1 and 2. 
For no. 3 through 6 only temperature data were available, and these wells were 
analyzed with the aid of sample logs from other wells, designated by letters on 


* Manuscript received by the Editor May 12, 1960. 

+ Humble Oil and Refining Company, Oklahoma City, Oklahoma. 
Published under the auspices of the Committee on Experimental Geology and Geophysics 
and the Division of Geological Sciences at Harvard University. 
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TABLE 1 
List oF WELLS 

Well Name Operator County State 
1 Lehigh Portland Cement Co. No.1 Jas. H. DuffandJ.T.Galey Butler Penn. 
2 Lewis Maxwell No. 11-F Columbian Carbon Co. Doddridge W. Va. 
3 Earl A. Hill No. 1 Sylvania Corp. Potter Penn. 
4 E. N. Clair No. 1 G. H. Rogers et al. Potter Penn. 
5 J. H. Lake No. 1 Hope Natural Gas Co. Marion W. Va. 
6 M. O. Goff No. 1 Hope Natural Gas Co. Harrison W. Va. 


Figure 1. All the wells lie in the Appalachian Plateau province, where the topo- 
graphic relief is moderate and the beds essentially flat-lying. 


ESTIMATING THERMAL RESISTIVITIES 


The key problem lay in estimating the thermal resistivities of the various 
strata from their descriptions in the logs. This was accomplished by first com- 
piling average resistivity values for different categories of sedimentary rocks. The 
averages were calculated from measurements made on rock specimens from 
various localities in the United States and Europe. Descriptions based on micro- 
scopic examination were available for all the specimens. As many categories as 


possible were set up. 
The resistivity measurements are given in Table 2 and summarized in Table 
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Fic. 1. Index map showing the locations of wells used in the study. 
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TABLE 2 
DATA ON SPECIMENS 


Thermal Temperature 


Resistivity of Porosity Density 
) incm °C Measurement in Percent in gm/cm* — 
sec/cal in °C 
A. QUARTZOSE SANDSTONES 
Georgia 80 30 0.5 2.64 Birch (unpublished) 
Massachusetts or 
New York 87 30 1 2.64 Birch (unpublished) 
New York 94 30 7 2.46 Birch (unpublished) 
Colorado 95 30 5 Birch (unpublished) 
Baltic 96 30 6 2.49 Birch (unpublished) 
? 125 30 26 2.02 Birch (unpublished) 
Pennsylvania 92 45 0.5 2.64 Clark (1941) 
B. IMPURE SANDSTONES 
Michigan 168 30 19 2.10 Birch (unpublished) 
Pennsylvania 125 30 20 2.41 Birch (unpublished) 
California 222 30 16 2.23 Birch (unpublished) 
California 149 30 7 2.50 Birch (unpublished) 
Pennsylvania* 120 2.63 Birch (unpublished) 
C. SANDY SHALES 
Pennsylvania* 189 30 - 2.76 Birch (unpublished) 
Pennsylvania* 185 30 7 2.76 Birch (unpublished) 
California 260 30 4 2.53 Birch (unpublished) 
D. SHALES 
Pennsylvania* 333 30 sad 2.71 Birch (unpublished) 
Pennsylvania* 278 30 = 2. Birch (unpublished) 
Pennsylvania* 333 30 - 2.55 Birch (unpublished) 
Pennsylvania* 286 30 = z.49 Birch (unpublished) 
Pennsylvania* 250 30 ” 2.73 Birch (unpublished) 
E. LIMESTONES 
Vermont 151 45 1 2.67 Clark (1941) 
Bavaria 162 45 3 2.30 Clark (1941) 
Indiana 197 45 13 2.61 Clark (1941) 
F. CALCAREOUS SHALES 
Pennsylvania* 204 30 We 2.53 Birch (unpublished) 
Pennsylvania* 238 30 wie 2.49 Birch (unpublished) 


* From the A. W. Bennet No. 1 well in Sullivan County. 
** Not measured. 


3. Four were done by Clark (1941) and the rest by Birch (unpublished). The de- 
scriptions by which the specimens were classified were made by J. F. Rominger 
or the writer. In these descriptions the mineral content was estimated; Rosiwal 
transverses were not made. The “‘quartzose sandstones” contained more than 95 
percent quartz (in all cases chalcedony and chert were lumped together with 
quartz); the “impure sandstones” had 70 to 90 percent quartz and the remaining 
material was primarily fine-grained matrix or feldspar; the “‘sandy shales” con- 
tained 20 to 35 percent quartz and no appreciable calcite; the “shales” had less 
than ten percent quartz, plus calcite; the “limestones” were essentially pure cal- 
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TABLE 3 
AVERAGE THERMAL RESISTIVITIES FOR DIFFERENT ROCK CATEGORIES 


Average Resistivity 
at 30°C in cm °C Range of Resistivity 


Quartzose sandstones 7 95 80-125 
Impure sandstones 5 157 125-222 
Sandy shales 3 211 189-260 
Shales > 296 278-333 
Limestones 3 162 144-188 
Calcareous shales 2 


221 204-238 


cite; and the “‘calcareous shales’’ contained from 15 to 35 percent calcite and less 
than five percent quartz. Of course, the statements above do not take account of 
quartz in particles too fine to be seen with the petrographic microscope, but the 
fine-grained rocks were tested with hydrochloric acid for the presence of calcite. 

During the resistivity measurements all of the specimens were compressed 
and all of the porous specimens but one were wet (Francis Birch, personal com- 
munication; Clark, 1941). Differences between the conditions under which the 
samples were measured and those prevailing in the wells at depth, except for 
temperature, are probably not sufficient to cause appreciable error. 

The average resistivities in Table 3 were corrected to the temperature in the 
wells before they were used. The corrections were less than the uncertainty,in the 
resistivity estimates and did not have an important effect on the results. The 
following formula was used: 


Resistivity at T°C=Resistivity at [1+ K(T—30) ], 


where 7 is the centigrade temperature and K the temperature coefficient. For 
shaley rocks K was given a value of 0.001, for limey rocks and quartzose rocks 
0.003. Resistivity temperature curves (Birch and Clark, 1940, p. 554) for a slate, 
a Solenhofen limestone, and a quartzose sandstone, respectively, were used to 
obtain the coefficients. The original resistivity measurements on the three lime- 
stones had been made at 45°C, and the values given in Table 3 were corrected 
to 30°C by the above formula. The rest of the measurements, with one exception, 
had been made at 30°C. 

It is difficult to evaluate the uncertainty in this method of estimating re- 
sistivities. The uncertainty is probably least in the “shales” and “calcareous 
shales.” The specimens for these two categories were all taken from the A. W. 
Bennet No. 1 well in Sullivan County, Pennsylvania. While this well is many 
miles away from the wells in which the temperature measurements were made, 
the specimens did come from the stratigraphic intervals in which those measure- 
ments were taken. 

Of the three specimens of “‘sandy shale,’’ two came from the A. W. Bennett 
well. One of these, however, is a red shale of the Catskill facies and is not simi- 
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lar to any of the rocks in which the temperature measurements were made. 
Estimating the uncertainty for the “‘sandy shales” is particularly difficult because 
there are only three samples. This is unfortunate since the heat flow values for the 
upper intervals in the wells are strongly influenced by the resistivity adopted for 
“sandy shales.” 

In each of the other categories the specimens came from scattered localities 
and represent a wide range of porosity. The spread of resistivity values within 
each category should therefore be a rough guide to the uncertainty. With two 
exceptions, all the specimens in these categories have resistivities within 25 per- 
cent of the mean. 

Everything considered, the writer believes that 25 percent is a safe upper 
limit for the error in all the categories. An independent check on the error will 
be presented in the discussion of the results. 


WELL-SAMPLE LOGS 


Detailed well-sample logs provided the basic lithologic information. The logs 
are much too lengthy for reproduction here, but they have all been published 
elsewhere (see Table 4). 

The rocks with which we are concerned lie entirely within the Devonian. 
Between the top of the Devonian and the top of the Onondaga limestone (Middle 
Devonian) are four to six thousand ft of shales and sandy shales. This interval is 
referred to as the ‘‘Devonian shales” by Martens (1939, p. 25), who describes 
the rocks as follows: 


“In a general way it may be said that the upper part of the Devonian Shales... 
is gray to greenish-gray and sandy and contains many interbedded thin layers of 
fine sandstone and siltstone. Going down from the top the amount of sandy ma- 
terial becomes less and the color darker. The basal portion of these shales is 
black. ...’’ In order to distinguish the green, sandy, upper portion from the 
lower, black shale portion we will use the term “‘Chemung-type” for the former 
and ‘‘Marcellus-type” for the latter. These names come from formations charac- 
teristic of the types. Below the Devonian Shales lies the Onondaga limestone, 


TABLE 4 
SouRCES FOR WELL-SAMPLE Locs 


Wells Source 

1 Martens, 1939, p. 48-53 

2 Martens, 1939, p. 103-111 
A Fettke, 1949a 

B Fettke, 1949b 

€ Fettke, 1938, p. 119-124 
D Fettke, 1938, p. 137-140 
E Martens, 1939, p. 109-111 
F Martens, 1945, p. 449-455 
G Martens, 1945, p. 471-477 
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which is 20 to 150 ft thick. (The Onondaga is probably the equivalent of what in 
West Virginia is called the Huntersville chert.) Below that in order come the 
Oriskany sandstone and the Helderberg limestone, each about 200 ft thick. 

The logs for wells no. 1 and 2 (the wells for which we have both temperature 
measurements and complete logs) were made by examining samples with a 
binocular microscope and treating them with dilute hydrochloric acid (Martens, 
1939, pp. 6-7). Both no. 1 and 2 were drilled with cable tools, and it is therefore 
unlikely that contamination of samples by caving from above has led to any 
serious errors in the logs (Martens, 1939, p. 5). In any case, contamination could 
not have influenced the results from the Marcellus-type shale intervals because 
the logs describe them as homogeneous. The heat-flow values from the Marcellus- 
type intervals in wells no. 1 and 2 do not differ significantly from the mean values 
for those wells. In the other wells only Marcellus-type intervals were used. 

All the sandstones within the Devonian Shales were assigned the thermal re- 
sistivity of ‘impure sandstones” (Table 3). There is some direct evidence to sup- 
port this interpretation. From a microscopic examination Martens (1939, p. 24) 
described six sandstone samples of the Chemung formation. Their quartz content 
ranged from 60 to 80 percent. The “impure sandstones”’ of Table 3 range from 
70 to 90 percent. The Bradford Third Sand, which Fettke (1949a, p. 2) calls 
“sandstone,” lies within the Devonian Shales. Kyrnine (1940, p. 55) studied the 
Bradford Third Sand in detail and described it as follows: 


60-65% individual quart grains 


13% feldspar 
30% rock fragments (quarizites, schists, phyllites, slates) 
2-3 fairly large mica flakes 


5-10% “clay paste.” 


There is no guarantee, of course, that the terms “‘sandy shale,” ‘‘calcareous 
shale,” etc., were used by those who compiled the well-sample logs in precisely 
the same sense as they were used in making up Table 3. But, since the categories 
are quite broad, the error from this source would probably not be large enough 
to justify raising the estimate of uncertainty above the 25 percent previously 
decided upon. 

Categories intermediate between those listed in Table 3 were necessary in 
order to make full use of the well-sample logs. Therefore, ‘‘very sandy shale”’ was 
given a resistivity halfway between that of “impure sandstone” and ‘sandy 
shale”; “slightly sandy shale” halfway between ‘“‘sandy shale” and “shale”; 
“slightly calcareous shale” halfway between “‘calcareous shale” and “‘shale’’; and 
“argillaceous limestone” halfway between “‘calcareous shale” and “‘limestone.” 
“Sandy limestone” was assumed to be the equivalent of 30 percent ‘“‘quartzose 
sandstone” and 70 per cent “limestone,” and its resistivity was an average of the 
“limestone” resistivity weighted by 0.7 and the ‘‘quartzose sandstone”’ resistivity 
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weighted by 0.3. Likewise, “‘calcareous sandstone”’ was assumed to be 70 percent 
‘“‘quartzose sandstone” and 30 percent “‘limestone.”” Chert was treated as the 
equivalent of macrocrystalline quartz. 

In using the logs it was necessary to introduce certain arbitrary rules. ‘‘Shale 
with a little interbedded sandstone” was taken to be 90 percent shale and ten 
percent sandstone. ‘‘Some interbedded sandstone,” or just ‘‘interbedded sand- 
stone,”’ was interpreted as 20 percent sandstone, ‘‘much” or ‘considerable inter- 
bedded sandstone” as 30 percent sandstone, ‘‘shale, partly sandy” as 50 percent 
shale and 50 percent sandy shale. Similar interpretations were made of analogous 
statements involving other rock types. ‘‘Silty shale’? was considered to be 
“slightly sandy shale” and “‘siltstone”’ to be ‘“‘sandy shale.” ‘‘Shale, partly sandy 
and silty’ was taken as 60 percent “shale” and 40 percent ‘‘sandy shale.” 

The well-sample logs were used in the following manner. For each interval in a 
well between adjacent temperature measurements a resistivity was determined by 
weighting the appropriate resistivity value for each bed by thickness and averag- 
ing. The temperature gradient was then divided by this resistivity to obtain the 
flow of heat through the interval. To avoid corrections for topography or climatic 
history only intervals from depths greater than 1,000 ft were used. 

It should be added that, in order to guard against bias in the process of 
estimating resistivities, the author refrained from analyzing the temperature 
measurements until all decisions affecting the average resistivity values had been 
made. 

RESULTS FROM THE FIRST GROUP OF WELLS 


Tables 5 and 6 show the results for wells no. 1 and 2. The spread of the indi- 
vidual values is encouragingly small. In Figures 2 and 3 the gradients and heat: 
flow values for wells no. 1 and 2 are plotted against depth. The gradients show a 
striking correlation with rock type; the heat-flow values show a corresponding but 
greatly reduced correlation. For intervals dominated by the Chemung-type rocks 


TABLE 5 
Heat FLow 1n WELL No. 1 


Depth of Interval Gradient in Average Resistivity Heat Flow in 
in feet 10°* °C/cm in cm °C sec/cal microcal/cm? sec 
1000-1500 233 224 1.04 
1500-2000 237 206 1,55 
2000-2500 313 254 We 
2500-3000 331 253 1.31 
3000-3500 371 275 £35 
3500-4000 353 299 1.18 
4000-4500 375 291 1.29 
4500-4809 189 200 0.95 
4809--5000 191 145 1.32 


mean heat flow (weighted by thickness) 1.20 microcal/cm? sec 
standard deviation 10% 
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TABLE 6 
HEAT FLow IN WELL No. 2 


Depth of Interval Gradient in Average Resistivity Heat Flow in 
in feet 10-* °C/cm in cm °C sec/cal microcal/cm? sec 


3500-3850 224 1 
3850-4000 243 1 
4000-4500 295 1 
4500-5000 313 gs 
5000-5500 270 : 1, 
5500-6000 361 1 
6000-6500 411 1 
6500-7000 197 1 
mean heat flow (weighted by thickness) 1.22 microcal/cm? sec 
standard deviation 11% 


the heat-flow values are in every case lower than the means. It is probable that 
the resistivities used for “impure sandstones” and ‘‘sandy shales” were too high 
for these particular rocks. 

Any error in the mean heat flow due to errors in temperature measurements is 
probably negligible compared with the uncertainty in estimating resistivities. 
That uncertainty was judged to be 25 percent and, on that assumption, the heat 
flow in wells no. 1 and 2 is 1.2+0.3 microcal/cm? sec. 

From other considerations, however, it appears that the heat flow cannot be 
less than 1.1 microcal/cm?* sec. The average gradient exceeds 360 10-* °C/cm 
over a 1,500-ft interval in well no. 1 and over a thousand-ft interval in well no. 2. 
A heat flow less than 1.1 microcal/cm*® sec would imply average resistivities 


THERMAL GRADIENT °C/CM 
190_200 1000 


THERMAL GRADIENT IN °C/CM 
100 200 7392 400 


CHEMUNG - TYPE 


CHEMUNG - TYPE ae 


TRANSITION 


IN| FEET 


FEET 


MARCELLUS - TYPE 


DEPTH 


DEPTH 


MARCELLUS - TYPE 


Top of Onondaga | 
HEAT FLOW 5000 , GRADIENT MEAT FLOW .. 
0 05 1.0 5 0 05 1.0 15 


HEAT FLOW IN MICROCAL/CM?® SEC HEAT FLOW IN MICROCAL/CM? SEC 


7000 


GRADIENT 


Fic. 2. Thermal gradients and heat flow Fic. 3. Thermal] gradients and heat flow 
in well no. 1. in well no. 2. 


1236 
‘ 
05 
10 i 
fecal 32 
24 
08 
41 
= 
07 
ae 
de 
2000 
4000 
ear 
3000 
| 
! 


HEAT FLOW IN PENNSYLVANIA AND WEST VIRGINIA 1237 


TABLE 7 


HEAT-FLOW VALUES IN WELLS No. 1 AND 2 FroM Groups OF INTERVALS DOMINATED BY 
DIFFERENT LITHOLOGIES 


Combined percent sandy Heat Flow in 


Well Depth of Interval in feet _ shale, impure sandstone, . 
and very sandy shale microcal/cm? sec 
1 1000-1500 62.4 1.04 
1 1500-2000 80.2 1.15 
2 3500-3850 61.7 1.05 


mean heat flow (weighted by thickness) 1.08 microcal/cm? sec 
Heat Flow in 


Well Depth of Interval in feet Percent shale microcal/cm? sec 
1 3500-4000 100.0 1.18 
1 4000-4500 86.2 1.29 
2 6000-6500 70.6 1.41 
mean heat flow (weighted by thickness) 1.29 microcal/cm? sec 
Combined percent chert, 
Well Depth of Interval in feet calcareous chert, cherty Heat Flow in 
limestone, and limestone microcal/cm? sec 
1 4500-4809 46.6 0.95 
1 4809-5000 87.4 1.32 
2 6500-7000 58.6 1.07 


mean heat flow (weighted by thickness) 1.08 microcal/cm? sec 


greater than 330 cm °C sec/cal in these intervals. But the highest resistivity of the 
five shale specimens (Table 2) was 333 cm °C sec/cal. These specimens were taken 
from a well in the same region as wells no. 1 and 2, and they came from the 
same stratigraphic interval as that in which the gradients greater than 360 
<10-*°C/cm were observed. 

A further check on the error was made by comparing heat-flow values from 
intervals dominated by different rock types. Three groups of three intervals 
apiece were used. In the first group were placed the three intervals with the high- 
est combined percentage of “sandy shale,” “impure sandstone,” and ‘‘very sandy 
shale,” in the second group the three intervals with the highest percentage of 
“shale,” and in the third group the three intervals with the highest combined per- 
centage of ‘‘chert,” ‘‘calcareous chert,” ‘‘cherty limestone,” and “limestone.” 
The results from the third group of intervals are completely independent of the re- 
sults from the other two. The only relationship between the first group and the 
second is that two of the three specimens used in obtaining the resistivity of 
“sandy shale” and all of the specimens of ‘shale’ came from the same well, the 
A. W. Bennet no. 1. The “‘shale” specimens, however, were black and nonsandy 
Marcellus-type rocks from 6,536 ft or deeper, whereas of the two ‘“‘sandy shales” 
one was a red shale of the Catskill facies from a depth of 2,360 ft and the other a 
green Chemung-type rock from 2,345 ft. 

Table 7 gives the results from the three groups of intervals; their consistency 
strengthens the conclusion that the error of the mean values for wells no. 1 and 2 
is not greater than 25 percent. From these arguments it is believed that the true 
heat flow lies between 1.1 and 1.5 microcal/cm? sec. 
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RESULTS FROM THE OTHER WELLS 


Temperature measurements but no sample logs were available for wells no. 3 
through 6. The average resistivities for intervals in these wells were obtained from 
the sample logs of nearby wells. Intervals were correlated from one well to an- 
other by distance above and below the top of the Onondaga limestone (or 
Huntersville chert). The Marcellus-type shales, because of their lateral uni- 
formity, seemed likely to give the best results, and so only intervals which lay 
mostly within those shales were used. 


TABLE 8 
Heat Fiow 1n WELLs No. 3 THROUGH 6 


Depth of Interval Gradient in Average Resistivity Heat Flow in 
in feet in 10-* °C/cm in cm °C sec/cal _microcal/cm? sec 


4000-4500 368 
4500-5165 402 
mean (weighted by thickness) Ls 


4000-4500 325 
4500-5000 364 
5000-5287 368 
mean (weighted by thickness) 


6000-6500 386 
6500-7000 284 
mean (weighted by thickness) 


6000-6500 
6500-7000 
7000-7310 


mean (weighted by thickness) 


6 
6 
6 
6 


Mean heat flow for wells no. 3 through 6 (weighted by thickness) 


In wells no. 3 and 4 average resistivities were determined from the sample log 
for well A, in well no. 5 from the log for well F, and in well no. 6 from the log for 
well E. As a check on the error involved in using the sample log of another well to 
termine the average resistivity of an interval, the average resistivities for inter- 
vals in well no. 4 were recomputed using the logs from wells B, C, and D and for 
intervals in well no. 6 using the logs from wells F, G, and no. 2. In no case did the 
recomputed value differ by more than 12 percent from the original. 

The results from wells no. 3 through no. 6 are presented in Table 8. Agree- 
ment with the results from the first two wells is satisfactory. 


DEVELOPING A MORE ACCURATE TECHNIQUE FOR ESTIMATING THERMAL 
RESISTIVITY IN SEDIMENTARY ROCKS 


The ideal would be a technique for estimating resistivity within ten percent 
from observations on samples possibly supplemented by electric, neutron, or 
velocity logs. Birch and Clark (1940) have accomplished this for igneous rocks 
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and Herrin and Clark (1956) have for rock salt. To achieve the desired accuracy 
for other sedimentary rocks more measurements on rock specimens would be 
necessary. Some information on porosity would also be needed. Linear regression 
analysis (Figure 4) on the three “limestones” and the six purest of the seven 
“‘quartzose sandstones’’ (Table 2) shows that a ten-percent change in porosity 
produces a resistivity change of about 25 percent. The porosity information could 
probably be obtained from analysis of electric logs or, preferably, velocity logs. 
It might be better to ignore porosity and work instead with the correlation be- 
tween resistivity and acoustic velocity. 
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Fic. 4. Thermal resistivity as a function of porosity for three pure limestones 
and six pure quartzose sandstones. 


The resistivities of coarse-grained rocks composed of quartz, calcite, and/or 
dolomite could probably be estimated within ten percent from their mineral com- 
position and porosity. More measurements on specimens would be needed, but 
the data on Figure 4 show the feasibility. The intercepts of the straight lines at 
the zero-porosity axis represent the effective resistivities of calcite and quartz at 
zero porosity. Using those effective resistivities an attempt was made to predict 
the resistivity of a nonporous, sandy limestone which contained 40 percent quartz 
and 60 percent calcite. An average of the effective resistivities weighted by vol- 
ume percent was computed, and it fell within five percent of the directly measured 
resistivity. This is essentially the method used by Birch and Clark (1940) for 
igneous rocks. 

Using a similar system the writer tried to predict the resistivities of the ‘“‘im- 
pure sandstones” (Table 2) from their quartz content and porosity. The effort 
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was a failure. A study is needed of the factors governing resistivity in rocks con- 
taining shaley material. X-ray diffraction analysis might help in such a study. 

With an accurate technique for estimating resistivity in sedimentary rocks it 
would be possible to obtain heat-flow determinations at a reasonable cost by 
making temperature measurements in selected wells after they had been aban- 
doned. Exploration for oil and gas, particularly in recent years, has been expand- 
ing to many new areas and provides an opportunity for good geographic coverage 
of the world. 

CONCLUSIONS 


This study showed that the heat flow at the wells investigated probably lies 
between 1.1 and 1.5 microcal/cm? sec. More important, it demonstrated the feasi- 
bility of getting heat-flow values accurate to better than 25 percent by estimating 
thermal resistivities from well-sample logs. 

With additional measurements on rock specimens and with some porosity in- 
formation, probably from velocity logs, it should be possible to reduce the error 
to ten percent or less for coarse-grained rocks composed only of quartz, calcite, 
and/or dolomite. To achieve such accuracy with the rocks containing shaly ma- 
terial would require a study of the factors governing their resistivity. A technique 
of estimating the resistivity of sedimentary rocks within ten percent would make 
heat-flow determinations in abandoned wells possible without the numerous cores 
now needed. 
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MAGNETIC AND TELLURIC CURRENT 
DISTURBANCES IN ALASKA* 


E. M. WESCOTT 


ABSTRACT 

Much of Alaska lies in or near the auroral zone, where magnetic and telluric current disturbances 
are large compared to those at lower geomagnetic latitudes. Pronounced and different diurnal varia- 
tions of disturbance are evident from analysis of K indices from Barrow, College, and Sitka magnetic 
observatories. The quietest period for magnetic or electrical surveying is the late morning and after 
noon hours locally. 

Recordings of telluric current disturbances can be used as a sensitive indicator of ionospheric 
activity. The large amplitude and rich frequency spectrum of the telluric currents suggests their use 
as a prospecting tool in Alaska. 


INTRODUCTION 


Basic research concerning high latitude ionospheric phenomena is conducted 
at the Geophysical Institute of the University of Alaska. Certain aspects of this 
research may prove to be of immediate practical interest to parties considering 
magnetic or electrical prospecting for oil or minerals in Alaska. In and near the 
auroral zone hardly a day passes without the occurrence of magnetic and telluric 
disturbance that would be considered a storm at lower geomagnetic latitudes. 
Since much of Alaska lies in the zone, exploration programs must be planned with 
knowledge of ionospheric phenomena. 


DEFINITION AND LOCATION OF THE AURORAL ZONE 


The auroral zone is generally defined as that zone or region in either hemi- 
sphere where the maximum number of aurorae occur, and is thus an average over 
some period of time. The area of maximum auroral activity may move north or 
south on various nights. The auroral zone is known to vary with the seasons, and 
depends upon sunspot activity. The center of the northern zone varies somewhat 
with longitude, ranging between 65° and perhaps to more than 70° geomagnetic 
latitude. 

During the IGY a number of all-sky cameras were set up to record the aurorae 
in Alaska. The films were analyzed for presence and motion of aurora by Davis 
and Kimball (1960). Figure 1 is a map of Alaska with the geomagnetic coordinates 
including a curve of the normalized presence of aurora projected along the 256° E 
geomagnetic longitude line. On the basis of these data the auroral zone projected 
parallel to the geomagnetic latitudes, lies between College and Barrow centered 
at about 67°. Whitham and Loomer (1957) who discuss the magnetic activity at 
Resolute Bay, Baker Lake, and Meanook in Northern Canada, show a map with 


* Manuscript received by the Editor February 10, 1960. 
¢ Geophysical Institute, University of Alaska, College, Alaska. 
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Fic. 1. The normalized presence of aurorae along the 256.56° geomagnetic longitude line. 
I g g 


the center of the zone north of Barrow. Barrow is at a lower geomagnetic latitude 
than Baker Lake, and Sitka is at about the same latitude as Meanook. 


MAGNETIC DISTURBANCES 

The diurnal variation of magnetic disturbance was obtained by analysis of 
the quasi-logarithmic A indices from the U.S. Coast and Geodetic Survey mag- 
netic observatories at Barrow, College and Sitka.' The lower limit for K=9 at 
Sitka is 1,000 y, compared to 2500 y at Barrow and College. In order to obtain 
ne a true average of the diurnal magnetic activity, the K indices were converted to 
a field intensity in gammas by means of the K to ak index conversion. The eight 
daily values were averaged on a monthly basis. The monthly averages were also 
grouped by equinoctial and solstitial periods e, 7, and d. 

The K indices are readily available, and thus are handy for determining mag- 
netic disturbance patterns, but they are by no means ideal. The A indices are an 
indication of world wide sudden commencements or solar events which may occur 
at any time distorting the normal daily disturbance pattern. Since the index rep- 


1 For a description of magnetic indices see Appendix 1. 
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Fic. 2. The diurnal variation of magnetic disturbance for the months of 1959, 


resents the maximum range of activity of H, D, or Z in a three-hour period, one 
cannot determine which component was the most active. No effort was made to 
refine the data, but probably the seasonal average curves are little affected by 
random events, 
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Fic. 3. The dirunal variation of magnetic disturbance for the equinoctial and solstitial periods, 1959. 


The effect of smoothing over a three-hour period, makes it very difficult to tell 
whether the diurnal patterns occur on local or geomagnetic time. Whitham and 
Loomer (1957) refer the times of maxima in their diurnal curves to geomagnetic 
noon and midnight, although a replot of these data gives a better fit with local 
noon and midnight. Because the present study is based on a limited number of sta- 
tions with small longitude separation, the diurnal curves are referred to local 
time. In spite of the shortcomings, the K index does give a good measure of mag- 
netic activity for a study of magnetic survey problems. 

The diurnal curves for Barrow, Figures 2 and 3, show a persistent morning 
peak during the period 08-11 hours, 150 WMT, corresponding to the predominant 
daytime peak characteristic of magnetic stations inside the auroral zone (Mayaud 
(1956)). There is also a nighttime peak of activity centered near local midnight, 
which is similar to the normal nighttime maximum characteristic of stations be- 
low the auroral zone such as Sitka and Meanook. The nighttime maximum varies 
in amplitude with the seasons, predominating in the winter, and diminishing in 
the summer. The development of the nighttime peak agrees with the results of 
Mayaud (1956) and Whitham and Loomer (1957), and suggests that the auroral 
zone moves north during the winter months. 

The diurnal curves at College have a predominant nighttime maximum near 
local midnight, and activity is high until about 0600 150 WMT. There is very 
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little indication of a daytime peak. In general the activity at College is less than 
at Barrow and particularly in the morning hours. 

At Sitka the diurnal curve has a single peak which occurs during the night- 
time. The amplitude of disturbance is in general about half that of College. 

Activity tends to be greatest during the winter months and least during the 
summer months at all three stations as demonstrated in Figure 3. This result has 
been verified on the basis of four years of telluric current data at College (Hessler 
and Wescott 1960). In the auroral zone zone between College and Barrow, one 
may expect larger disturbances. 


MAGNETIC SURVEYS IN ALASKA 


Morley (1953) makes the observation that if corrections from a ground con- 
trol station are applied to aeromagnetic profiles made during a storm, the flight 
line must be less than 87 miles from the ground station to have any meaning. Be- 


TABLE la 
No. oF Days 1n 1959 K <2 (1007) 
Hours 150 WMT 02-05 05-08 08-11 11-14 14-17 17-20 20-23 23-02 
College 95 123 184 228 198 179 139 93 


TABLE 1b 
No. oF Days 1959 K <4 (3507) 


05-08 08-11 11-14 14-17 17-20 20-23 23-02 


Hours 150 WMT 02-05 
College 218 251 314 349 349 323 274 204 


176 


Barrow 212 


207 206 287 326 306 241 
cause the form of a magnetic storm varies with distance from a station, a map 
produced by correcting profiles made during storms would probably be so inac- 
curate as to be of little use in oil prospecting. For accurate magnetic surveying it 
seems best to fly only during the quiet periods. 

Table 1 shows the number of days in 1959, when the magnetic activity was 
less than 100 y, and 350 ¥, in the three hour periods. Data are given only for Col- 
lege and Barrow, as the activity at Sitka is much less overall. 

During the quietest hours, about one day in four shows activity less than 100 y 
at Barrow. One day in two is this quiet at College, while at Sitka most of the days 
were as quiet at this time. At all stations the afternoon activity was generally less 
than 350 y. On the basis of the average diurnal pattern of disturbance, and the 
number of quiet days, the best period for conducting surveys is the afternoon 
hours. 

TELLURIC CURRENT DISTURBANCES 


Telluric currents are continuously recorded by the Geophysical Institute at 
College and Barrow. There is a high degree of correlation between earth currents 
and magnetic disturbance. At College the telluric currents are generally oriented 
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Fic. 4. Scatter diagram and least squares regression line. Telluric current activity versus 
magnetic activity at College, April 1958. 


along a line N 35° W, so for scaling purposes only the N-S component is consid- 
ered. The telluric current records are scaled for range of activity in mv/km dur- 
ing the same three hour periods as the K indices. 

To obtain a measure of the correlation between electric and magnetic ac- 
tivity, the daily average telluric current range in mv/km was compared with the 
magnetic equivalent daily amplitude A,, in gammas. Figure 4, shows the scatter 
diagram mv/km vs A..(y) for April 1958. The correlation coefficient is 0.975. 

The diurnal disturbance variation of telluric currents at ‘College is shown for 
the three seasons of 1959 in Figure 5. As one would expect the curves are quite 
similar to the corresponding magnetic curves. 

The earth current records show rapid fluctuations of considerable amplitude. 
During the largest storm at College since 1956, potential gradients greater than 
5,000 mv/km were recorded. 


ELECTRICAL PROSPECTING 


Due to the large, rapid variations in telluric currents near the auroral zone, 
prospecting techniques requiring measurements of natural or artificial d.c. po- 
tentials may be seriously hampered during storms. From knowledge of the ex- 
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Fic. 5. Diurnal variation of telluric current activity for the solstitial and equinoctial periods, 1959. 


pected size of telluric potential disturbances, and the expected accuracy of po- 
tential measurements one can determine the feasibility of working during active 
periods. 

The telluric current disturbances of appreciable amplitude at College have 
periods of about six seconds and longer (Hessler and Wescott 1959). They are 
often nearly sinusoidal in form, although irregular wave shapes predominate, 
especially at the longer periods. The large amplitudes, regular occurrence, and 
rich frequency spectrum of the telluric currents suggests their utilization in 
telluric or magneto-telluric prospecting in Alaska. Also large areas of flat sedi- 
mentary rocks fulfill the requirements of uniform telluric current sheets. 

Aside from being a possible nuisance, or aid to electrical prospectors, the 
telluric currents offer a sensitive indicator of ionospheric activity. All that is re- 
quired for measuring them is a pair of electrodes in ground free from artificial 
sources, an amplifier and recorder. Such a recorder can be set up almost anywhere 
and operated at less expense than a recording magnetometer as an indication of 
general ionospheric activity. 

CONCLUSIONS 

With this present state of knowledge it is not possible to predict magnetic 
storms very far in advance. Often storms will occur about 48 hours after the ob- 
servation of a large flare on a favorable location on the sun’s surface. At College 
medium to large storms are often observed to repeat for several nights at about 
the same time. If a very large storm occurs at night, the next afternoon is quite 
often disturbed. At each magnetic station a diurnal variation of disturbance is 
clearly evident, with significant differences between stations. The quietest period 
occurs locally during the late morning to afternoon hours. 

Magnetic and electrical surveying can best be planned with a knowledge of 
ionospheric activity. Recordings of telluric currents have proven to be a sensi- 
tive indicator of activity. The simplicity of recor: . 1g, makes them a natural aid in 
planning survey operations near the auroral zo. -. 
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APPENDIX 1. MAGNETIC INDICES 

Current information on the time variations of geomagnetic activity is given 
by three parameters for a particular observatory: 

1. the three-hour range index, K, 

2. the daily equivalent amplitude, As, and 

3. the daily character-figure, C. 

The three-hour range index, K, is derived from the records of the three or- 
thogonal field components, for instance H, D, and Z. For each three-hour period, 
beginning at Greenwich midnight, the range of the disturbance variation in the 
most disturbed component, is determined with a quasi-logarithmic A scale from 
0 to 9. Each observatory specifies the lower limit for K =9, and then uses a simple 
multiple of a normal table to obtain the range limits for the A indices. The 
normal table is: 


K=0 1 2 3 4 5 6 7 8 9 
Range = 0-5 5-10 10-20 20-40 40-70 70-120 120-200 200-330 330-500 5007. 


Thus for College, where the lower limit for K = 9 is 2500 y, the range limits would 
be multiplied by 5. 

To convert a K index reported by an observatory into an equivalent range 
ak, another standard table is used: 


K=0 1 2 3 4 5 6 7 7 9 
ak=0 3 7 15 27 48 80 140 240 400 (2). 
To find ak in gammas for any observatory, this table is multiplied by a factor ob- 
tained by dividing the lower limit for K=9 by 250. For instance at College the 
iower limit for K=9 is 2,500 y, so the factor is 10. Thus for K=4 at College, 
ak = 270 y. The equivalent daily amplitude, Ak, is the average of the eight values 
ak per Greenwich day. 

The daily character-figure C is assigned to each Greenwich day on the basis of 
C=0, if it is quiet, C=1, if it is moderately disturbed, and C= 2, if it is heavily 
disturbed. The method of assigning these characters is left entirely to the discre- 
tion of the directors of each observatory. 


REFERENCES 


Bartels, J., 1957, The geomagnetic measures for the time variations of solar corpuscular radiation, 
described for use in correlation studies in other geophysical fields: Annals of the I.G.Y., v. 4 
parts 4-7, p. 227-236. 


fa 
1249 
= 
f 


1250 EUGENE WESCOTT 


Davis, T. N., and Kimball, D. S., 1960, Incidence of auroras and their north-south motions in the 
northern auroral zone: Geophysical Institute, University of Alaska Scientific Report n. 4, 
N.S. F. Grant n. Y/22.6/327. 

Hessler, V. P., and Wescott, E. M., 1959, Correlation between earth-current and geomagnetic dis- 
turbance: Nature, v. 184 supplement n. 9, p. 627. 

——— and Wescott, E. M., 1960, Earth current disturbance data, College, 1956-1958: Geophysical 
Institute, University of Alaska, Scientific report n. 2, Contract n. AF 19(604)-3075. 

Mayaud, P. N., 1956, Acivité magnétic dans les regions polaires: Ann de Géophys. v. 12, p. 84-101. 

Morley, L. W., 1953, The areal distribution of geomagnetic activity as an aeromagnetic survey prob- 
lem near the auroral zone: Trans. A.G.U., v. 34, p. 836-840. 

Whitham, K., and Loomer, E. I., 1956, Characteristics of magnetic disturbances at the Canadian 
arctic observatories: Pub. Dom. Obs., v. 18, n. 12. 

1957, Irregular magnetic activity in northern Canada with special reference to aeromagnetic 
survey problems: Geophysics, v. 22, n. 3, p. 646-659. 


§ 
4 
| 
i 
4 
| 
tt 
i 
|. 
| 
| 
| 5 
: 
be 
| 
7 
4 
é 
a 
ey 
ih 
at 


GEOPHYSICS, VOL. XXV, NO. 6 (DECEMBER, 1960), PP. 1251-1259 


INDEX OF WELLS SHOT FOR VELOCITY* 
(NINTH SUPPLEMENT) 


V. U. GAITHERT 


ABSTRACT 


This is the ninth supplement to the original Index published in the 1944 issue of GEOpHysICs 
Information is listed on 340 well velocity surveys not heretofore reported in any previous publication 
of the Index. Corrections and/or additional information on 5 previously listed surveys are also 
tabulated. 


INTRODUCTION 


This is the ninth supplement to the original ‘‘Index of Wells Shot for Veloc- 
ity,”’ published in the October 1944 issue (Volume IX, Number 4) of GEopuysics. 
The publication of this supplement was made possible only through the coopera- 
tion of the entire geophysical industry, for which the editors of Geopuysics and 
I wish to express our sincere thanks. 

The surveys tabulated in the first portion of this supplement are those shot 
since the Eighth Supplement (published in the December 1959 Yearbook of 
Geophysics) ; also older surveys not reported in previous publications. In the lat- 
ter portion, shown in italics, are corrections and additional information on previ- 
ovsly listed surveys. 

It is requested that any errors which may be found in this index be reported 
to the Business Manager of the Society for correction in subsequent supplements. 

The writer is grateful to Mrs. Kathryn B. Smith for her assistance in recon- 
ciling and assembling information for this supplement. 


* Manuscript received by the Editor September 22, 1960. Copies of this Supplement are obtain- 
able from the SEG Business Office at fifty cents each. Copies of the Fourth, Fifth, Sixth, 
Seventh and Eighth supplements are also available at fifty cents each. 

+ Continental Oil Company, Ponca City, Oklahoma. 
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DEPARTMENTS 
PATENTS 


O. F. RITZMANNT 


MAGNETIC PROSPECTING 
U.S. No. 2,938,164. M. Hansburg. Iss. 5/24/60. App. 6/25/57. Assign. U.S.A. 
Automatic Earth Field Balance Control System for Magnetometers. A system for cancelling the 
ambient field affecting a magnetic detector in which the detector output is filtered by a low-pass 


filter and fed to an amplifier with degenerative feedback from the detector, so that signals passing 
through the filter are suppressed and only higher frequency anomalies are indicated. 


U.S. No. 2,942,177. J. Neumann, H. Weiss, and F. Assmus. Iss. 6/21/60. App. 5/6/55 and 5/1/56. 
Assign. Siemens-Schuckertwerke Aktiengesellschaft. 


Method and Means for Measuring Magnetic Field Strength. A device for measuring ferromagnetic 
properties of material by placing a wafer of semiconductor material with high carrier mobility edge- 
wise against the material to be tested and observing the change in transverse electrical resistance 
of the wafer or its Hall-effect voltage. 

U.S. No. 2,942,179. F. B. Coker. Iss. 6/21/60. App. 11/12/53. Assign. United Geophysical Corp. 


Magnetometer Circuit. A magnetometer having a flux-gate element with compensating winding 
that is connected in series with a resistor in a feedback loop arranged to compensate the field to 
which the element is exposed, the voltage across the resistor being measured by a self-balancing 
potentiometer-recorder. 

U.S. No. 2,942,180. F. B. Coker. Iss. 6/21/60. App. 11/12/53. Assign. United Geophysical Corp. 


Magnetic Field Measuring System. A magnetometer having a pair of flux-gate elements in an un- 
balanced bridge circuit whose output produces positive and negative pulses that are amplified and 
the difference in amplitudes between successive pulses of opposite polarity indicated. 

RADIOACTIVITY PROSPECTING 
U.S. No. 2,935,614. C. F. Teichmann and C. G. Ludeman. Iss. 5/3/60. App. 5/17/54. Assign. 

Texaco Development Corp. 

Radioactive Prospecting. An airborne radioactivity surveying system in which a camera periodi- 
cally photographs the terrain and a radioactivity detector collimated to receive radiation from the 
same area is connected to a counting rate meter and amplifier to make a visual record of observed 
radioactivity on the same or on a correlated film. 

U.S. No. 2,939,010. S. A. Scherbatskoy. Iss. 5/31/60. App. 12/23/53. 
Directional Radiation Detector. A directional detector for hard gamma rays having three scintil- 
lating crystals in line each with a photomultiplier tube connected to a triple coincidence circuit. 
SEISMIC PROSPECTING 
U.S. No. 2,935,697. L. B. McManis. Iss. 5/3/60. App. 7/18/57. Assign. Pan American Petroleum 
Corp. 
Seismic Amplifier Gain Control. A seismograph avc system in which the output of the seismic 


+ Gulf Oil Corporation, Patent Department. 
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channel is amplified, rectified, and used to control diode bridge attenuators in two preceding stages 
of the channel, the degree of attenuation for the early stage of the channel being higher than for a 
later stage. 


U.S. No. 2,935,723. W. J. Stark and R. L. Wilcox. Iss. 5/3/60. App. 11/9/56. Assign. Jersey Produc- 
tion Research Co. 


Multiple Geophone Connector. An electrical connector having an insulated support with a num- 
ber of polarized and interconnected terminal posts to which geophone wires and jumper cable can 
be connected. 


U.S. No. 2,935,728. A. R. Morgan. Iss. 5/3/60. App. 2/15/46. Assign. Radio Corp. of America. 


Underwater Object Locator. A directional object locator having a piezoelectric transducer excited 
with frequency modulated radio frequency pulses and the frequency of the reflection compared with 
that of the source. 


U.S. No. 2,937,915. R. A. Peterson. Iss. 5/24/60. App. 12/19/55. Assign. United Geophysical Corp. 


Recording System. A channelized galvanometric recording system in which each trace is optically 
repeated by light spots of lower intensity on each side of the main oscillograph beam, so that when 
the main beam exceeds the boundary of the channel at high amplitude the repeat traces still show 
the signal within the channel. 


U.S. No. 2,938,195. S. J. Begun, G. P. Chan, and A. P. Dank. Iss. 5/24/60. App. 7/8/54. Assign. 

Clevite Corp. 

Multichannel Magnetic Recording. A system for changing the time phase between the tracks of a 
multichannel magnetic record by re-recording the record with the tracks on separate drums that are 
adjustably mounted on a common shaft so that the time phase of the drums can be changed on play- 
back. 


U.S. No. 2,938,592. C. J. Charske and C. L. Hubbard. Iss. 5/31/60. App. 3/23/53. Assign. Jersey 
Production Research Co. 


Seismic Velocity Well Logging Apparatus. An acoustic well logging system having a source and 
two detectors embedded in a rubber sonde, the source being energized by a condenser at relatively 
long intervals and the difference in arrival time at the two detectors converted to a d-c signal that 
is recorded at the surface. 


U.S. No. 2,938,670. V. W. Bolie. Iss. 5/31/60. App. 2/28/57. Assign. Collins Radio Co. 


Means for Measuring Correlation. A correlation computer in which the two functions are applied 
to parallel resistors in series with an output resistor whose power is measured by a heat detector that 
controls a power source whose output is indicated. 


U.S. No. 2,939,079. P. L. Willmore. Iss. 5/31/60. App. 10/21/57. Assign. Canada. 


Method of Calibrating an Electromagnetic Seismograph. A system for calibrating a seismometer by 
connecting it in one arm of a bridge circuit and connecting the seismometer galvanometer as the 
bridge-balance detector, balancing the bridge with the seismometer coil clamped, and then observing 
the galvanometer deflection with a-c of various frequencies applied to the bridge when the seismome- 
ter is unclamped. 


OSS, No. 2,939,107. J. D. Ball. Iss. 5/31/60. App. 7/1/57. Assign. Jersey Production Research Co. 


Recordation of Seismic Signals. A calibrated reflection seismograph avc system in which the shot 
energy is determined by a rectifying and integrating circuit whose input is cut off after the first 
arrivals and whose output is applied to a potentiometer that controls the seismograph amplification, 
the potentiometer being mechanically varied with time by a connection to the recording drum. 
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U.S. No. 2,940,061. R. G. Piety. Iss. 6/7/60. App. 10/24/55. Assign. Phillips Petroleum Co. 
Signal Recording System. A seismic recording system in which the signals are amplified and 


differentiated and the amplified signal applied to positive and negative gate circuits that are con- 
trolled by the differentiated signal, the outputs of the two gate circuits being separately recorded. 


U.S. No. 2,940,062. J. P. Lindsey. Iss. 6/7/60. App. 2/6/56. Assign. Phillips Petroleum Co. 

Tuning System, A seismograph filter system having a noise rejector circuit that rejects damped 
wave trains and followed by a tuning network that is sensitive to waves having the form of the second 
derivative of the normal density function. 


U.S. No. 2,940,536. C. B. Bazzoni, L. G. Ellis, and A. C. Winterhalter. Iss. 6/14/60. App. 3/2/53" 

Assign. Sun Oil Co. 

Method and Apparatus for Anaysis of Seismic Records. An apparatus for producing a corrected 
record from a magnetic seismogram in a re-recording process in which the pick-up heads are adjusted 
so as to provide fixed corrections and also moved to correct for normal stepout and for stepout due 
to a dipping reflecting bed. 

U.S. No. 2,941,184. H. H. Moody. Iss. 6/14/60. App. 4/27/56. Assign. Continental Oil Co. 

Method for Plotting Seismic Data. A method of automatically transferring seismic data from a 
seismogram to a time or depth section in which the traces are scanned in pairs and each time a pre- 
determined condition of amplitude or slope reversal is simultaneously recognized on both traces a 
pulse is generated that is plotted on a synchronized tape. 


U.S. No. 2,941,615. T. H. Holmes. Iss. 6/21/60. App. 5/1/56. Assign. Texaco Inc. 

Device for Loading Tandem Charge Arrays. An anchoring device for the bottom of a shot hole 
having a hollow cylindrical container with outwardly flaring prongs that engage the wall of the 
hole and with a smaller transversely mounted cylinder around which a firing line with spaced charges 
can be threaded so as to permit adjustment of the charges to the desired depths. 


U.S. No. 2,942,195. W. C. Dean. Iss. 6/21/60. App. 5/15/58. Assign. Gulf Research & Development 

Co. 

Electrical Filter Circuits. Seismograph filter circuits made up of extended ladder sections with a 
running ground and with each section having a response of the Laguerre type, the outputs from a 
number of sections being combined to provide a readily computable overall transient response 
characteristic. 


U.S. No. 2,942,462. H. M. Buck. Iss. 6/28/60. App. 6/5/57. Assign. Jersey Production Research Co 


Device for Measurement of Pressure or Temperature in Wells. A go-devil for recording a charac- 
teristic of well fluid in a flowing well having a sensing device recording on a drum that is advanced 
stepwise each time a finger detects a tubing joint when the device is dropped in a well and finally 
releases an annular flexible packer so that the flow of well fluid will propel the device back up to the 
surface. 


U.S. No. 2,942,782. R. B. Blizard, R. Proskauer, and V. Vacquier. Iss. 6/28/60. App. 6/11/54 and 

1/26/55. Assign. Sperry Rand Corp. 

Sonar Ray Tracer. An analog computer for computing a family of ray paths of sound in a medium 
from measured values of sound velocity as a function of depth which is entered into a cathode ray- 
phototube type of function generator, the computed ray paths being displayed on a c-r tube screen. 
U.S. No. 2,943,258. E. N. Shawhan. Iss. 6/28/60. App. 8/12/54. Assign. Sun Oil Co. 


A pparatus for Impedance Measurements. An impedance measuring circuit for use with a capaci- 
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tance transducer having a long connecting cable and in which the exciting oscillator is switched at 
a low frequency between the transducing capacitor and a series-connected reference capacitor and the 
difference between the two parts of the cycle detected. 


U.S. No. 2,943,268. C. Dodge. Iss. 6/28/60. App. 7/30/57. Assign. Texaco Inc. 


Automatic Gain Control Amplifier Circuit. A seismograph avec system having two feedback loops 
each with a pair of silicon diodes and a resistor in a voltage divider circuit and each loop having a 
rectifier and filter whose output is used to control the impedance of the diodes. 


WELL LOGGING 
U.S. No. 2,934,947. H. M. Buck. Iss. 5/3/60. App. 8/8/55. Assign. Jersey Production Research Co. 


Flowmeter. A well flowmeter having a rotor on a vertical shaft with flow deflectors above and be- 
low the rotor to direct the flow against the rotor blades and with the shaft carrying a coded contactor 
that completes a circuit to the surface. 


U.S. No. 2,934,955. P. H. Monaghan. Iss. 5/3/60. App. 7/22/57. Assign. Jersey ?roduction Research. 

Co. 

Apparatus for Measuring Temperatures. A borehole temperature indicator having a number of 
vials of insulating materia] with contacts at their upper ends and a solid insulating meltable material 
of high density at the top of the vial and a conducting liquid of low density at the bottom of the vial 
so that when the heavier material melts the lighter liquid will close an indicating circuit to the surface. 


U.S. No. 2,935,615. M. E. True. Iss. 5/3/60. App. 6/24/57. Assign. Jersey Production Research Co. 


Well Logging Tool Fluid Dis placer. A radioactivity logging tool in which the measuring apparatus 
is inside an expansible packer that is inflated by gas pressure under control of a four-way solenoid 
valve actuated from the surface. 


U.S. No. 2,935,871. R. G. Pearson. Iss. 5/10/60. App. 4/10/56. Assign. Geolograph Service Corp. 


Well Logging Apparatus. A drilling rate recorder having two recording drums one of which re- 
cords traveling block movements when not drilling and the other records progress of drilling and is 
disengaged for other operations. 


U.S. No. 2,936,375. A. S. McKay. Iss. 5/10/60. App. 12/27/55. Assign. Texaco Inc. 


Radioactivity Bore Hole Logging. A neutron-gamma-ray well logging system in which a caliper log 
of the hole is simultaneously made with the voltage from the caliper tool effecting the recording of a 
reference curve for a borehole of the observed diameter to calibrate the gamma-ray curve. 


U.S. No. 2,937,278. G. V. Copland. Iss. 5/17/60. App. 3/20/57. Assign. Welex, Inc. 


Directional Scintillation Well Logging Apparatus. A radioactivity detector having an axial scintil- 
lating crystal and phototube with the upper portion of the axial crystal surrounded by a scintillating 
crystal and the lower portion of the axial crystal surrounded by another crystal with phototube, the 
phototubes being connected to an anti-coincidence circuit so as to record only horizontally directed 
radiation. 


U.S. No. 2,937,333. E. Boucherot. Iss. 5/17/60. App. 12/26/56 and 12/23/57. Assign. Societe de 
Prospection Electrique Procedes Schlumberger. 


Bore Hole Telemetering Apparatus. An electric logging signalling system in which a resistor equal 
to the load is connected to the conductors when the circuit is interrupted so as to eliminate transient 
effects. 


U.S. No. 2,937,527. D. F. Mason, G. Thodos, and L. F. Stutzman. Iss. 5/24/60. App. 4/9/56. Assign. 
Strike Investrrent Co. 
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Apparatus for Determining the Production Potential of Wells. A system for transmitting bottom- 
hole pressure or temperature to the surface by having the transducer move a Selsyn generator con- 
nected by a cable to a repeater at the top of the well. 


U.S. No. 2,938,117. K. H. Schmidt. Iss. 5/24/60. App. 3/23/56. Assign. Petroleum Service and Re- 
search Corp. 
Analysis Determinative of Gas or Oil Producing Strata. A system of mud logging in which gas is 
extracted from the circulating mud stream and passed through a multicomponent infrared absorption 
analyzer and the concentration of components separately recorded. 


U.S. No. 2,938,708. J. J. Arps. Iss. 5/31/60. App. 9/19/57. 

Simultaneous Drilling and Electrical Logging of Hydrocarbon Contents of Formations. A system of 
electric logging while drilling in which the self-potential of formations is measured by means of elec- 
trodes on an insulated drill collar and the salinity of the drilling mud adjusted to be the same as that 
of the fluid in the formation being drilled. 


U.S. No. 2,940,177. J. M. Bricaud. Iss. 6/14/60. App. 2/2/56. Assign. Schlumberger Well Surveying 

Corp. 

Apparatus for Exploring Bore Holes. A well surveying apparatus having a pool of conducting liquid 
in an insulated container with electrodes in its outer wall and also having a suspended magnet, the 
electrodes being scanned by a motor-driven contactor and the resistance measured and the azimuth 
of the magnet determined by a coil driven by the same motor. 


U.S. No. 2,940,370. J. L. Yandell. Iss. 6/14/60. App. 12/11/58. 

Bottom Hole Camera. A photographic device run on the end of well tubing and having a tank with 
lights and camera on its lower end and having by-pass valves so that clear water may be discharged 
from the tank for taking a picture or the tank by-passed for raising and lowering the device in the 
well. 


U.S. No. 2,941,405. P. F. Southwick. Iss. 6/21/60. App. 4/10/56. Assign. Gulf Research & Develop- 
ment Co. 

Mud Sampler and Pressure Indicator. A device for obtaining well fluid samples at known depth 
having a series of chambers each having a seal that is broken at a known pressure to allow fluid io 
enter the chamber and a float valve that closes the chamber when full, the frangible seals of the suc- 
cessive containers being connected in an electric circuit that signals the operation to the operator by 
a change of resistance. 


U.S. No. 2,041,784. J. L. Martin. Iss. 6/21/60. App. 7/5/55. Assign. The Atlantic Refining Co. 


Logging While Drilling. A system of electric logging while drilling using insulated sections in the 
drill string and a down-hole recorder connected to metallic sections. 


U.S. No. 2,942,111. A. E. Worthington. Iss. 6/21/60. App. 6/28/55. Assign. California Research 
Corp. 


Contoured Logging Sonde. A gamma-gamma density logging sonde in which both the incident and 
the scattered beams are collimated and the sonde has protuberances at the collimating openings so as 
to insure that these openings are always in contact with the borehole wall even thouga the latter may 
be rough and irregular. 


U.S. No. 2,942,112. D. P. Hearn. Iss. 6/21/60. App. 12/28/56. Assign. Well Surveys, Inc. 


Multiple Well Logging System. A system for transmitting a number of signals over a single con- 
ductor by means of a modulated carrier on which is superimposed positive and negative pulses, the 
signals being represented by the modulation, the positive spikes, and the negative spikes. 
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U.S. No. 2,942,461. E. R. Filley. Iss. 6/28/60. App. 12/16/55. Assign. Texaco Inc. 


Well Bore Hole Logging. A thermostated boot for calibrating a radioactivity logging device, the 
boot being heated either electrically or by circulating engine exhaust gas. 


U.S. No. 2,943,197. L. de Witte. Iss. 6/28/60. App. 10/14/54. Assign. Continental Oil Co. 

Method of Well Logging. A method of analyzing well fluid to locate points of entry of salt water 
by directing a stream of neutrons along the axis of the well and detecting only gamma rays of a par- 
ticular energy due to neutron capture by chlorine or sodium. 


U.S. No. 2,943,200. J. A. Rickard. Iss. 6/28/60. App. 10/19/56. Assign. Jersey Production Research 

Co. 

Stable Isotope Method for Well Logging. A radioactivity logging method using a neutron source 
and gamma-ray detector that is connected to a pulse-height discriminator that selects pulses indica- 
tive of neutron capture by two isotopes of carbon or sulphur and recording the ratio of their counting 
rates to determine the isotopic ratio. 


MISCELLANEOUS 


U.S. No. 2,940,039. L. Yost, W. R. Wilson, and R. A. McCreary. Iss. 6/7/60. App. 6/10/57. Assign. 
one-half to A. O. Smith Corp. and one-half to California Research Corp. 

Well Bore Electriccl Generator. A tachometer generator for a turbine drill having an annular per- 
manent magnet rotor on an extension of the turbine shaft and a sealed stator whose winding is con- 
nected to the drill pipe and to an insulated electrode, the a-c signal being picked up at the surface 
from connections to the drill pipe and a spaced ground electrode and its frequency indicated. 


U.S. No. 2,942,176. E. R. Brownscombe and E. M. Connell, Jr. Iss. 6/21/60. App. 5/7/57. Assign. 

The Atlantic Refining Co. 

Method and A pparatus for Determining Characteristics of Earth Formations. An apparatus for de- 
termining the formation factor of a parallel-sided rock sample in which the sample is placed between 
the ends of two coaxial capillary tubes filled with the rock-saturating liquid and the electrical resist- 
ance across the sample is measured. 


U.S. No. 2,943,319. W. R. Hunsicker. Iss. 6/28/60. App. 5/3/55. Assign. Seismograph Service Corp. 
Radio Location System. A continuous-wave phase-comparison type of radio location system in 
which the fixed stations use the same radiating antenna for both coarse and fine position indicating 
signals. 
U.S. No. 2,943,320. J. E. Hawkins. Iss. 6/28/60. App. 10/4/55. Assign. Seismograph Service Corp. 
Radio Location System. A continuous-wave phase-comparison type radio location system in which 
three fixed stations radiate signals of different frequencies and the signals are heterodyned in pairs 
and the beat frequency used to modulate the third signal from which one side band of the modulation 
is eliminated. 


U.S. No. 2,943,321. M. A. Karpeles. Iss. 6/28/60. App. 1/6/58. Assign. International Telephone and 
Telegraph Corp. 
Course Guidance System. A radio navigation computer in which polar coordinate position informa- 
tion and wind direction and velocity information obtained from a TACAN radio beacon is converted 
to a null-type indication to guide the pilot along a preset straight line course. 
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Patent No. Subject* Patent No. Subject* Patent No. Subject* 


2,934,944 496 2,936,388 2,938 ,117* 172, 196, 252 
2,934,947" 148 2,936,401 2,938,119 308 
2,934,951 148 2,936,419 2,938,121 
2,934,953 200 2,936,436 2,938, 123-125 
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2,936, 627-628 2,938 , 204-205 
2,936, 950-951 2,938,208 
2,936,953 Re. 24,829 
2,936,956 2,938,382 
2,937 ,005 2,938,383 
2,937 ,064 2,938 , 390 
2,935,617 2,937 ,239 2,938, 592* 
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2,935,740 2,937 ,292 2,938,964-965 224 
2,935 ,742-743 2,937 ,330 2,938 ,966 
2,935 , 867 2,937 ,333* 116, 516 2,939 ,010* 
2,935,871 2,937 ,368 228 2,939 011-014 
2,935,879 2,937 ,524 188 2,939,079 
2,935,960 2,937,525 460 2,939,082 
2,936,114 2,937 ,527* 288, 292 2,939,105 
2,936, 119-120 2,937 , 531-533 2,939, 106 
2,936, 134-136 2,937 , 562 2,939, 107* 
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* A key to the subject classification system will be found in Geopxysics, v. 12, pp. 256-264 (April, 
1947). 
a Abstracted on preceding pages of this issue. 
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Reviews of Resistivity Determinations from Electric Logs and Electric Analogue for Resistivity 
Logging by Hubert Guyod, and of Analysis Charts for the Determination of True Resistivity from Elec- 
tric Logs by Hubert Guyod and John A. Pranglin, appeared in the June issue, pages 701 and 702. Ob- 
jection to the opinions expressed in these reviews was received from the principal author, Dr. Guyod, 
and others. Independent reviews were accordingly invited by the Editor from Mr. G. E. Archie, the 
Associate Editor for WELL Loccrne. These reviews are presented herewith. 


Analysis Charts for the Determination of True Resistivity from Electric Logs, Hubert Guyod and John 
A. Pranglin, Hubert Guyod, Houston, 1959, 202 pp., $22.50. 
Resistivity Determinations from Electric Logs, Hubert Guyod, Hubert Guyod, 1957, $20.00. 
These manuals give a wealth of information on estimating the true formation resistivity from 
electric logs comprising normal resistivity, lateral resistivity, and SP curves which record only ap- 
parent resistivity data. The electrical resistivity of a rock is a function of its porosity, fluid saturation 


1268 


3 
ate 
Bh 
aig 
a 
| 
| 
| 
| 
: 
| 
| 
a 
= 
: 
iy 


REVIEWS 1269 


and other factors. In order to determine rock characteristics in this manner, therefore, one must be 
able to determine the true resistivity. The manuals do not explain how to estimate the reservoir 
characteristics once the true resistivity is known. 

The first book presents a method for determining the true resistivity and the extent of the in- 
vaded zone from the 16-inch and 64-inch normals, and the 18-foot 8-inch lateral curves. The method 
is applicable when R, (invaded zone resistivity) is greater or smaller than R; (true formation resistiv- 
ity). An analogue computer was used to obtain most of the data. The book contains over 250 charts 
that permit determining R, and, in many cases, R,; and D, (depth of filtrated invasion) for single iso- 
tropic beds from two to 50 feet thick, and for comprehensive ranges of resistivities “R;,” “R:,” “Ry” 
(shale resistivity), and R,, (mud resistivity) hole sizes and invasion diameters. Instructions are given 
for the application of data to certain nonuniform formations. A discussion of the resolving power of 
resistivity curves and of the effects of inaccuracies in the various pertinent factors Ro, Rm, R., € (thick- 
ness), d (borehole diameter) is presented. 

The second book was originally prepared as a field manual for six major oi] companies. This 
manual gives the background necessary for the practical use of resistivity curves and the field data for 
determining the resistivity of single beds and of multiple-layer formations. 

The two books should, therefore, be used together. Typical logs are given for single formations, 
multiple- and sandwich-type formations. Familiarity with typical characteristic shapes of various 
resistivity curves as recorded through layers of these geometries is important for interpretation. 

The large number of variables encountered in attempting to arrive at the true resistivity of a 
layer makes it necessary for the specialist to have a full understanding of a range of simple cases to 
more complicated ones, as well as being fully aware of inaccuracies. Therefore, these manuals provide 
excellent background for the specialist who wants to take the time necessary to find the best estimate 
of the true resistivity of a layer in order to approach in-situ reservoir conditions. Although normal 
and lateral resistivity logs are being replaced by focused resistivity logs, older well logs recorded on 
nearly all wells since the 1930's will be more important in the United States in the future because of 
the emphasis on secondary recovery. 

G. E. Arcute 
Shell Development Company 
Houston, Texas 


Electric Analogue for Resistivity Logging, Hubert Guyod, Hubert Guyod, Houston, 1958, 137 pp., 
$50.00. 

This book summarizes work done with an analogue computer for six major oil companies. Cal- 
culated data for resistivity well logging utilize the relationship between apparent resistivities, Ra, 
true resistivities, R,, and other pertinent factors. It is generally impossible to state this relationship 
as an explicit expression giving true resistivity. Therefore, the usual method for determining the re- 
sistivity of beds from electric logs is to establish charts relating R, and R; and to use these charts to 
analyze the field data. In the past these charts were generally obtained by calculations, using a method 
which has the disadvantage that it is very time-consuming. 

The report comprises two parts: The first part describes at length the construction of the ana- 
logue and its operations; well bore and rock formations were simulated by electrical resistances. The 
second part is a thorough summary investigation of the response of normal and lateral resistivity 
curves on single and multiple beds. It discusses the apparent resistivity, shapes of curves, and methods 
of analysis. 

This book is not necessary for the field man; however, it should be extremely useful in research 
departments, universities and log analysts’ libraries, and valuable to anyone who may be confronted 
with similar problems that may be solved on an analogue computer. The principle, construction, and 


operation are discussed in detail. 
G. E. ARCHIE 


Shell Development Company 
Houston, Texas 
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High-Frequency Filtering as a Means of Eliminating Multiple Reflections, A.M. Epinat’eva and 
L. A. Ivanova, Bulletin (Izvestiya), Academy of Sciences, USSR, Geophysics Series, No. 3, March 
1959, pp. 244-252, translated and published by The American Geophysical Union, May 1960. 


The two Russians who collaborated on this experiment set out to design a seismic system that 
would provide better discrimination against certain multiple reflections than the conventional sys- 
tems previously used. Their success with a high-frequency system (peak = 105 cps) is worthy of note. 
In this reviewer’s opinion, however, they may be seriously in error in attributing their success entirely 
to the frequency response of the system used. 

The authors chose to conduct the majority of their experiments in an area where the basement 
depth is only 1600 to 1700 meters, excellent for high-frequency techniques. Four geophones per group 
spaced five meters apart along the profile were used in the high-frequency system, whereas the con- 
ventional medium frequency system (30 to 50 cps) utilized only one geophone per group. Recordings 
at long ranges (up to 1000 meters) were made, and the normal moveout observed on these long range 
records was used as the basis for distinguishing interbed multiple reflections from primary reflections. 
At the high frequencies used, the geophone groups will be especially effective in discriminating against 
the higher normal moveout of the interbed multiples. No such discrimination would be obtained from 
the single geophone used in the medium frequency system. (The geologic section is conformable with a 
dip of less than one degree, so dip moveout is very small.) Thus, the dimensions of the geophone 
groups may be as important a factor in their successful system as the frequency response. 

Obviously, very shallow reflections at long ranges would have the greatest normal moveout so 
would be most heavily attenuated by the geophone groups. This problem was neatly sidestepped by 
arbitrarily ignoring all events occurring earlier than a “reference reflection” at a depth of 750 to 850 
meters. 

In justification of using geophone groups, the following statement appears: “Grouping of geo- 
phones increased the intensity of the recordings 3 to 4 times. The trace appeared also more regular 
than with single geophones, and no addition or loss of waves occurred.” (???) 

It should be pointed out that the simple multiple of the “reference reflection” was neatly elimi- 
nated by the high-frequency system. Examination of the records presented appears to indicate that 
the geophone group dimensions are not a contributing factor. 

This article, incidentally, is quite difficult to read because the translator obviously is not well 
acquainted with seismological terms. 

Tuomas L. SLAVEN 
Western Geophysical Company of America 
Los Angeles, California 


Electric Prospecting by the Telluric Currents Method, A. M. Zagarmistr and M. N. Berdichevskiy, 
Petroleum Geology, Vol. 3, No. 1-B, 1959 (English edition published June 1960), pp. 40-49. 


The English edition of the original Russian paper is clearly written and informative. It is divided 
into three parts, (1) introduction section covering the general aspects of the telluric-current method 
together with its relation with other geophysical methods, (2) regional telluric-current investigations 
of reconnaissance nature in the poorly studied Kanta-Mansiysik depression of the western part of 
West Siberian Lowland, and (3) telluric-current exploration of local geologic structures within larger 
geological structural units. 

Russian interest in the telluric-current method of prospecting for oil and gas is made evident by 
the revelation that 111,200 sq km of Russia were covered by this exploration method during the three- 
year period, 1955-1957 inclusive. 

Very few details of the method as used in Russia are reviewed; however, a concise and informative 
general statement of the cosmic origin of telluric currents and the required technique of simultaneous 
measurements of two orthogonal components of the telluric field at two distant (10 km) points is given 

The second part of the paper deals with the initial geophysical reconnaissance coverage of wide 
areas of unknown regional and local geology. A specific telluric-current survey of the northern part of 
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the Shatilkov depression is briefly reviewed. It is pointed out that subsequent seismic work conducted 
in the same area has confirmed the information derived from the method of telluric currents. 

The third part of the paper reviews exploration work by the telluric-current method, within a 
small-scale portion of the Shatilkov depression. It is pointed out that the local structures were studied 
by merely increasing the density of the network of telluric-current observations. 

Several interesting sketch maps and data profiles, combining geological data with results obtained 
by the telluric method, are presented. The agreement between the telluric-current method measure- 
ments and geological information is emphasized. It is reported that the method gave good geological 
results in a majority of the cases. Apparently the method was less effective in areas where thick ce- 
mented deposits outcrop. 

The conclusion of the paper is that the use of the telluric-current method of exploration will or 
has become an integral part of the regional geophysical investigations of the broad areas in Siberia and 
Central Asia. 

This paper is of value to the geophysicist interested in the possibilities and limitations of the 
telluric-current method of exploration. 

Josuva L. SosKE 
Stanford University 
Stanford, California 


On the Application of Cagniard’s Method to Dynamical Problems of Elasticity, M. Mitra, Gerlands 
Beitrige Zur Geophysik, Vol. 69, No. 2, pp. 73-86, 1960. 


The paper by M. Mitra uses first the Laplace transform with respect to the time, then the Fourier 
transform with respect to the horizontal distance, and then finally the Laplace transform again with 
respect to the vertical distance or the distance perpendicular to the interface, thus reducing the prob- 
lem he considers to an algebraic one. The paper is an interesting variation of Cagniard’s method ap- 
plied to the problem of reflection and refraction of two-dimensional seismic pulses. 

C. Hewitt Dix 
California Institute of Technology 
Pasadena, California 


Space Technology, edited by Howard S. Seifert, John Wiley & Sons, Inc., New York, 1959, 1172 pp., 
$22.50. 


This remarkably excellent volume is based on a series of postgraduate lectures in space technology 
offered by the Extension Division of the Department of Engineering and Physical Sciences of the 
University of California. It is interesting to note that in spite of the highly analytical nature of this 
course, over four thousand students (including this reviewer) registered for this series at the various 
branches of the University. In his preface to the book, Dr. Seifert states: ““That such an intense inter- 
est was shown in the subject was not unrelated to the fact that Sputnik I was launched a few months 
before the course began.”’ 

The subject of space technology is discussed in five broad categories, subdivided into thirty-three 
chapters and written by thirty-eight authors, each of whom has gained recognition in the field of his 
specialty. The five categories discussed are as follows: 

Flight Dynamics: This section of eleven chapters is largely a space-age adaptation of the classical 
disciplines of analytical dynamics and celestial mechanics. After introducing the fundamental con- 
cepts and definitions suitable for the proper parametrization of rocket flights, the authors treat the 
theory and the problems of optimization of terrestrial, satellite, lunar, and interplanetary trajectories. 
Some attention is given to the common fallacy of the simple balancing of Newtonian fields and this 
leads quite naturally to a discussion of the velocity contours based on the very important Jacobi’s 
integral for the restricted three-body problem. Looking into a somewhat more distant future, a full 
chapter is devoted to the well-known relativistic time dilatation effect, leading to the inescapable con- 
clusion that for sufficiently high velocities of travel, astronauts will indeed age more slowly than those 
who have remained “‘at home.” This section closes with a discussion of the thermal problem of re- 
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entry into the atmosphere as well as the more general problem of achieving nondestructive landing 
on a given planet. 

Propulsion and Structures: The seven chapters of this section unfold an exciting drama of man’s 
ingenuity (though, as yet, no final solution) in the face of natural physical laws which are ostensibly 
unsympathetic with his cosmic goals. 

A designer would like the terminal velocity of his rocket to be as high as possible and at the same 
time to retain a large fraction of the original mass. An elementary consideration, however, of the con- 
servation of momentum between the rocket and its exhaust, leads to the following fundamental equa- 
tion of rocketry with which he must contend, namely, Terminal Velocity = Exhaust Velocity Xlog, 
(Mass Ratio), the last being the ratio of the initial total mass to the final total mass. Obviously, then, 
beyond a certain point a small increase in the terminal velocity has to be paid for by a terrible sacri- 
fice in the final mass, which of course includes the payload. This impasse leads the designer to con- 
sider the remaining parameter, viz., exhaust velocity. Here, more sophisticated considerations of 
momentum and adiabatic expansion of gases lead to a very cumbersome equation which, in part, 
states that the exhaust velocity is proportional to »/7/M, where T is the combustion temperature 
(°K) and M is the mean molecular weight of the exhaust gases. Now, then, if one can supply heat 
by means other than chemical, one can also use substances of low molecular weight such as hydrogen 
or helium, and thus gain a substantial overall advantage. While this opens new and unlimited vistas 
of approach, it also unfolds a myriad of unsolved problems and engineering difficulties. Dr. Seifert 
summarizes this whole situation in the following classic statement: ‘Since an effective rocket should 
spend mass frugally, it must squander energy prodigally.” 

The topics in this section include the following: fundamentals of nozzle flow and optimum nozzle 
geometry (which is by no means self-evident); liquid, solid, nuclear, and thermonuclear rocket pro- 
pulsion (with an excellent chapter on the fundamentals of magnetohydrodynamics) ; analysis of struc- 
tural problems; and, finally, application of the systems approach to the integrated design. 

Communication, Guidance, and Control: A variety of topics is covered in the six chapters of this 
section. Weight limitations place an upper bound on the permissible transmitter power in the space- 
craft. On the other hand, because distances may well get very large and in view of the inverse square 
law, the inherent noise of the receiver is an immediate problem, and sets the lower bound on the usa- 
bility of a given transmitted signal. This problem is treated within the framework of information 
theory which is very briefly sketched. Next, the theory, advantages, and disadvantages of the follow- 
ing modes of guidance are discussed: radio guidance, inertial guidance, radio-inertial guidance, and 
the new subject of midcourse and terminal guidance. An interesting subsection treats the topic of 
smoothing a stationary time series in the context of guidance problems. 

Man in S pace: This phase of space technology is perhaps more familiar to general readers because 
of the numerous articles which have appeared in national magazines. The four chapters of this section 
analyze man in terms of his psychological and physical responses to weightlessness, confinement, stress, 
radiations, and accelerations. The section includes a chapter on the design of space cabins. 

A pplications of Space Technology: Two chapters discuss the use of satellites for research in space 
physics, some aspects of special instrumentation, and a brief sketch of methods of data handling. 
The third chapter is a panel discussion: What the Future Holds. 

This reviewer has been asked on several occasions why the name of von Braun was not included 
among the authors of the book. Dr. von Braun was, in fact, one of the invited speakers. However, 
instead of appearing in person, as all the other busy lecturers have done, he chose rather to send a 
““warmed over’’ film version of his ‘““Exploration of Mars,”’ with himself as the narrator. A fair portion 
of the audience considered this an affront, as was evidenced by vociferous sounds of disapproval. The 
editor is to be commended for excluding this material from his book for, in this reviewer’s opinion, it 
lacked the two most important qualities present in all the other lectures: freshness and sophistication. 

It is warmly suggested that all earth-men whose thoughts are not earth-bound ought to have a 
copy of this unique and excellent volume on their reference shelves. 

Mark HoizMan 
Western Geophysical Company of America 
Los Angeles, California 
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Nomography, A. S. Levens, John Wiley & Sons, Inc., New York, 1959 (Second Edition), viii+-296 
pp., $8.50. 

To those who have used the first edition (1948) no words of introduction are needed. Improved in 
format, revised and expanded by about 90 percent, and including three new chapters on circular 
nomograms, projective transformations, and the relationship between concurrency and alignment 
nomograms, the new edition is an up-to-date treatment of a subject which is continually enlarging 
its field of application. 

A comparison with the first edition shows that although the basic structure is unchanged, addi- 
tions have been made to almost every chapter. Noteworthy features of the new edition include the 
application of nomographic methods to test the validity of experimental data curves, an expanded 
treatment of the “method of determinants,” a simple approach to projective transformations, and 
the Duality Principle (‘‘for every line of the concurrency nomogram there is a corresponding point in 
the alignment nomogram’’) and its application to the transformation of families of experimental data 
curves from concurrency to alignment charts. The Appendix which previously contained twenty-nine 
nomograms has been doubled and covers applications to engineering, statistics, physics, food tech- 
nology, medicine, etc. 

The book can serve admirably both as a text and asa reference. 

H. KAUFMAN 
McGill University 
Montreal, Canada 


Modern Probability Theory and Its A p plications, E. Parzen, John Wiley & Sons, Inc., New York, 1960, 
xv +464 pp., $10.75. 


Despite the increasing importance of the theory of probability to the physical, social and engineere 
ing sciences, there has been a dearth of texts combining a mathematically sound treatment with well- 
chosen technical applications. One of the best to date has been Feller’s Introduction to Probability and 
Its Applications. Parzen’s book is a worthy addition to the growing list of acceptable volumes. 

The first six chapters constitute a course in elementary probability theory. Beginning with the 
notions of “event,” “sample description space,” and the basic definition of probability as a function 
of events on a sample description space, the author carefully proceeds to develop the concepts of in- 
dependence and dependence, numerical-valued random phenomena, mean and variance of a prob- 
ability law, and the normal, Poisson and other distributions. The theory is profusely illustrated with 
well-chosen examples. The unhurried development is attested by the fact that not until chapter 7 is 
the concept of “random variable” introduced. 

The final two chapters are of a considerably more advanced nature, providing careful and rigor- 
ous derivations of the law of large numbers and the central limit theorem. A valuable feature is the 
inclusion of sets of exercises, theoretica! and nontheoretical (numerical) at the end of each section. 

Those of us who have been confronted with the rudimentary type of text on the one hand, and 
the formidable treatises of Doob and Loéve on the other, will be amply rewarded by a thorough study 
of this book. 

H. KAUFMAN 
McGill University 
Montreal, Canada 


An Outline of Geomor phology, S. W. Wooldridge and R. S. Morgan, Longmans, Green, and Co., Ltd., 
London, 1959, xix4-409 pp., 32 shillings. 


This is the eleventh edition of a book first published in 1937. It is designed as a textbook for ad- 
vanced geography students, and its purpose is to describe topography in terms of geological processes. 
The illustrations, 271 in number, plus the wealth of material within its 23 chapters, make this book 
most interesting reading. It is probably of considerable value to those who are encountering geology 
for the first time and to whom geology as such is of secondary importance. To geologists and geo- 
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physicists, many imperfections should be apparent. The original title was “The Physical Basis of 
Geography, An Outline of Geomorphology.” It is unfortunate that the authors reversed it, as it places 
an improper emphasis on the book as a text on geomorphology. 

A statement on the cover says that “. . . the bibliography has been thoroughly revised to include 
the most recent papers and texts.” Perhaps so, but the text itself does not show evidence of much re- 
vision. In more than a dozen places, it is apparent that recent work has been overlooked. An example 
of this ‘‘. . . recently there have been significant developments in geophysical prospecting, involving 
the use of detailed torsion balance surveys. . . .”” Furthermore, in an eleven-page discussion of the 
controversial subject of continental drift, there is no mention of evidence of drift and polar wandering 
from paleomagnetic measurements. 

The first ten chapters deal with the geological basis or background for the remaining 13 chapters 
which pertain to the configuration of the earth’s surface. Subject matter includes “The Origin of the 
Earth,” “Isostasy,” ‘“Mountain Building,” “Origin of the Forces Responsible for Earth Movement,” 
“The Cycle of Erosion,” “‘Denudation Chronology,” ‘Marine Erosion and Shore-Lines,”’ and ‘‘Glacial 
Control in the Physical Geography of the Earth.”’ The phrases in quotations are the titles of eight of 
the chapters. 

It is the reviewer’s opinion that this book, with weaknesses as described above, is a valuable addi- 
tion to the library. 

James AFFLECK 
Gulf Research & Development Company 
Pittsburgh, Pennsylvania 


The Geological Evolution of North America, T. H. Clark and C. W. Stearn, The Ronald Press Company, 
New York, 1960. vi+434 pp., $7.50. 


The authors have done an excellent job of presenting a logically developed, smooth reading intro- 
duction to historical geology. The book reads more like a novel than a text book. There are no refer- 
ences in the text itself though selected references are listed at the end of each chapter. 

“The historical development of the continent is present in terms of its three major structural 
units: the bordering geosynclines, the stable interior, and the Canadian Shield. A greater unity of 
treatment is, we believe, achieved by considering the development of these great structural units 
throughout geologic time, or through a major part of it.” 

The introductory chapter discusses geologic time and illustrates the principles involved in litho- 
logic and paleontological correlation. Theories regarding the origin of the earth and the solar system 
are briefly reviewed. 

As much of historical geology is the story of sedimentation and tectonics of geosynclines, the 
authors provide a simplified explanation of the geosynclinal theory. They define a geosyncline as “‘an 
elongated belt in the earth’s crust that is subsiding currently or that has subsided for long periods in 
the past so that it acts, or has acted, as a trap or collecting basin for sgliments removed by erosion 
from uplifted areas of the crust. Axial zones of maximum subsidence receive great thicknesses of sedi- 
ments until they become unstable and the belt is then crushed by foiding and elevated into a moun- 
tain chain.” 

After describing the geologic history of the stable interior (including the origin of such famous oil- 
producing formations as the St. Peter, Oriskany, Chester, Mission Canyon, and the Silurian and De- 
vonian reefs), the history and evolution of the marginal geosynclines—the Appalachian and Cordil- 
leran—are presented in considerable detail. The stratigraphy, lithology, and structural evolution of 
these marginal geosynclines are traced from the Precambrian to the Recent. 

These chapters are followed by an up-to-date review of the geology of the Arctic Islands. The 
complex geologic history of the Canadian Shield, or the continental nucleus, and its role in the growth 
of the North American Continent are discussed. The final topic of this section of the book is The His- 
tory of the Great Ice Age. 
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The last four chapters concern the evolution of life on earth as shown by the fossil record. An 
appendix contains a short summary of the Plant and Animal Kingdoms. 

The book is profusely illustrated with photographs and simplified diagrams which are extremely 
effective. In describing the history of the more important formations, the authors point out which 
ones are major hydrocarbon producers. 

C. E. WEAVER 
Continental Oil Company 
Ponca City, Oklahoma 


New Methods of Elevation Control Speed Reconnaissance Gravity Surveys in Northern Areas, 
Norman R. Paterson, The Canadian Mining and Metallurgical Bulletin, Vol. 53, No. 580, pp. 
614-622, August 1960. 


The paper gives a quantitative comparison of three methods of gravity surveying which do not 
require running levels between stations. The comparisons were undertaken to evaluate the methods 
which might be used in areas such as northern Canada where ground transportation is very difficult 
and where reconnaissance gravity surveys with lower accuracy than is obtained with usual field pro- 
cedures could give useful results. Transportation by helicopter can carry gravity meters readily to 
station locations under such conditions but if instrumental surveying for horizontal control and or- 
dinary spirit leveling are used to determine positions and elevations the economy of helicopter trans- 
portation is largely lost. The paper, therefore, is concerned with methods of operation in which loca- 
tions are determined from aerial photographs or maps and elevations by other means than spirit level- 
ing or are avoided. 

Elevations can be determined by photogrammetric methods with the elevations along the line 
of flight controlled by the Airborne Profile Recorder (APR) which is a precise radar altimeter with 
readings corrected by a hypsometer which determines variations from the normal pressure elevation. 
The paper discusses the types and magnitudes of errors which may be encountered. The example used 
covers an area of about 15 by 20 mi. northeast of Toronto with local relief of approximately 500 ft. 
The photography was at a height of 20,000 ft and scale 1:40,000 with the photogrammetric determi- 
nations made in a Wild Model A7 stereo-triangulation instrument. The results were compared with 
spirit level elevations at the same points. A histogram of the error distribution shows that 77 percent 
of the elevations had errors of less than 5 ft and 5 percent had errors greater than 10 ft. 

A test over the same area was made with barometric altimeters. The various methods of altimeter 
operation are discussed. The test was made using the altimeter in much the same way as the gravity 
meter, i.e., by assuming a linear pressure change between base readings and without using recording 
barographs. The author concluded that most of the errors are due to uncompensated instrument drift 
and hysteresis and to local pressure changes due to wind velocity. The instrumental errors (drift and 
hysteresis) possibly would have been reduced materially by carrying several altimeters and averaging 
the results. The histogram of the errors shows that 65 percent are better than 5 ft and 11 percent are 
greater than 10 ft, with a mean square error of 5.6. 

A test of a “gravity gradient” system carried out by the Robert H. Ray Company is included. 
This system depends on a triangle of three accurate gravity stations with distances of the order of 
1,000 ft between stations and calculating the magnitude and direction of the average gravity gradient 
from the gravity differences and geometry of the three stations. A gravity map can be calculated by 
integrating the average gradient from these measurements in the same way that gravity maps are 
computed from torsion balance data. The only leveling required is that between the stations of each 
group and elevations from group-to-group are not needed. This system, therefore, could be used by 
selecting favorable local sites for each group and using helicopter transportation between these loca- 
tions and no intermediate surveying would be required. A test of this system was carried out in Lee 
County, Texas and compared with gravity meter readings. On the basis of some 20 stations from 1 
to 2 mi apart, the average differences between the gravity from gravity meter and from the inte- 
grated gradients was about 0.65 mg. 
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The paper concludes with a table summarizing the precision, rate of progress, and costs for vari- 
ous systems. Since this table is a good summary of the entire paper, it is reproduced here. 


Gravity Approx. Final 
Survey ” Average Speed | Cost for 2000— 
| Final Mean le 
Specifications Elevation Control Method | Accuracy of Survey 5000 Square 
(station 7 (milli Is) (square Mile Survey 
spacing in en's, miles/day) |(dollars/square 
miles) | mile) 
3X1 | Photo—transit 0.05-0.10 10— 20 50-70 
1X2 Photo—transit “0.05-0.10 15-30 | 40-60 
1X2 | Photogrammetry. Existing | 0.50-1.00 | 80-100 | 10-15 
1:60,000. Existing A.P-R. and 
| ground control. 
1X2 | Photogrammetry. Existing | 0.30-0.50 80-100 | 20-30 
| 1:60,000. New A.P.R. 
1x2 | Photogrammetry. New 0.20-0.30 80-100 30-40 
| 1:40,000 and A.P.R. | 
1X2 | Altimeter. Base looping. 0.50-1.00 80-100 812 
1X2 | Altimeter. Base looping and 0.40-0.80 | 60-80 | 15-20 
| base recording. 
1X2 | Altimeter. Leap-frog and base 0.10-0.50 30— 50 30-50 
| recording. 
1X2 | Horizontal Gravity Gradients. 0.50-1.00 40— 50 20-30 


L. L. NETTLETON 
Gravity Meter Exploration Company 
Houston, Texas 


Die Karte der mittleren Freiluftanomalien fiir Gradabteilungen 6'X 10’ von Westdeutschland (Map of 
Mean Free Air Anomalies for a 6’X10’ Coordinate Grid of West Germany), Karl Gerke and 
Heinz Watermann, Publication No. 32, Institute of Applied Geodesy, Frankfort, Germany, 1959, 
18 pp.+map. (This is also Heft No. 46, Teil II of Reihe B, of the Deutsche Geodatische Kom- 
mission bei der Bayerischen Akademie der Wissenschaften.) 


As the title indicates, the authors have calculated mean free-air gravity anomaly values in West 
Germany and environs in a grid of areas about 7 X7 mi square. The values, which range from —52 to 
+94 milligals, are shown to the nearest milligal on a map of scale 1:1,000,000. The results are not 
contoured on the map but a color scale indicates areal ranges in values in 20-milligal steps. 

A brief text (in German with summaries in German, English and French) describes the details of 
the calculation procedure and precision of the results. The English summary is especially compre- 
hensive. References are given to source materials and related publications, mostly in German. There 
is no discussion of the geophysical applications of the results which presumably will follow in subse- 
quent publications. 

This paper is one of a series of gravity studies of central Europe by the Institute of Applied 
Geodesy in the German Geodetic Commission. A companion Bouguer anomaly gravity map was pub- 
lished in 1957. The associated map of average ground elevations, required in the calculations of the 
free-air anomalies, was published in 1959. (See accompanying review.) 
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The area reported covers all of central Europe from Belgium and the northeastern tip of France 
on the west, Poland on the east, Denmark and southern Sweden on the north, and small portions of 
northern Italy and Yugoslavia on the south. This series of publications furnishes very useful data for 
regional studies of isostasy, tectonics and structure in the area. 

StcMUND HAMMER 
Gulf Research & Development Company 
Pittsburgh, Pennsylvania 


Karle der mittleren Hohen von Zentraleuropa (Map of Mean Elevations in Central Europe), Alfred 
Schleusener, Publication No. 31, Institute of Applied Geodesy, Frankfort, Germany, 1959, 8 pp. 
+map. (This is also Heft No. 60 of Reihe B, of the Deutsche Geodatische Kommission bei der 
Bayerischen Akademie der Wissenschaften.) 


This paper consists essentially of a map showing average ground elevations in a grid of rectangular 
areas of about 125 square kilometers each. A brief text (in German, with summaries in English, French 
and German) describes the construction of the map and gives references to original sources. The map 
covers all of central Europe from Belgium and the northeastern tip of France on the west, Poland on 
the east, Denmark and southern Sweden on the north and small portions of northern Italy and Yugo- 
slavia on the south. 

The grid areas have been conveniently chosen to be nearly square (6 ft latitude X 10 ft longitude). 
The average elevations were computed from available topographic maps and represent the average 
topographic elevation in each grid area within an accuracy of +2 to +10 meters. In some of the 
bordering areas and especially in the high Alps the errors may be considerably larger. The grid ele- 
vations range from 0 to nearly 3,000 meters. 

The map was constructed as a necessary computational step in preparing an average free-air 
gravity anomaly map of central Europe which has been published as a companion paper by Gerke 
and Watermann. (See accompanying review.) Its publication will facilitate gravitational studies of 
geodetic and regional tectonic problems in the area. 

SicMUND HAMMER 
Gulf Research & Development Company 
Pittsburgh, Pennsylvania 
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BOOKS 


Physics of the Upper Atmosphere, edited by J. A. Ratcliffe, Academic Press, New York, $14.50. 

Soviet Research in Geophysics v. 1, 2 & 3, American Geophysical Union, Washington, D. C., v. 1, 
$8.00, v. 2, $6.50, v. 3, $6.00. 

Crystal-Structure Analysis, by Martin J. Buerger, John Wiley & Sons, Inc., New York, $18.50. 

Statistical Theory of Communication, by Y. W. Lee, John Wiley & Sons, Inc., New York, $16.75. 

Lead Isotopes in Geology, by R. D. Tussell and R. F. Farquhar, Interscience Publishers, New York, 
$9.00. 
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JAMES LOUIS CATT 


James Louis Catt died suddenly of a heart attack on Easter Sunday, April 
17, 1960, at his home in Tulsa, Oklahoma. Funeral services were held on April 19 
and interment was at Lyons, Kansas. The news of his untimely death came as 
a shock to his family, his associates, and his many friends. 

Jim, as he was better known, was born October 15, 1914, at Cashion, Okla- 
homa. His family later moved to Lyons, Kansas, where Jim attended Lyons 
High School, graduating in 1932. In 1935 he entered Sterling College at Sterling, 
Kansas. He graduated in 1939 with a Ph.B. in mathematics. 

Jim started to work in the training program of Seismograph Service Corpora- 
tion in January, 1941. In December, 1942, he entered the Army of the United 
States. He was assigned to the 654th Engineer Topographical Battalion where 
he eventually became First Sergeant of Co. A. He was discharged in December, 
1945, and returned to Seismograph Service Corporation in January, 1946. He 
was made party chief in 1946. He worked in most of the oil provinces of the 
United States and also spent some time in Canada in 1948 and 1949. 

Jim left Seismograph Service Corporation in 1955 to join the newly formed 
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Geo Prospectors, Inc., as a seismologist. He was elected to offices in Prospectors 
Drilling Corporation and Geo Exploration, Inc., by his associates. He held these 
positions until the time of his death. 

Jim is survived by his wife, Viola, and by two daughters, Janet and Marla at 
homie, by his sister, Miss Cora Catt, of Larned, Kansas, and by his mother, 
Mrs. Zett Catt, of Lyons, Kansas. 

His passing is keenly felt by his many friends and associates. His friendliness, 
ability, and integrity will be an inspiration to all of us who knew and worked 
with him. 

B. G. PRESTON 


E. E. ROSAIRE 


Esme Eugene Rosaire was born in Chicago, Illinois on June 17, 1897, to 
Claude and Julie (Carrol) Rosaire. After graduating from Tuley High School, he 
entered college at the University of Chicago, but his studies were interrupted 
by the war. In 1917 he enlisted in the Air Service, was commissioned a Second 
Lieutenant in 1918 and served in France with the 138th Aero-Squadron, Fifth 
Pursuit Group, until his discharge in 1919. He returned to college and received 
his Bachelor of Science degree in 1920; this was followed by the Master of Science 
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degree in 1921, and in 1926 the University of Chicago granted him the Doctor of 
Philosophy degree in Chemistry and Physics. 

From 1923 until 1925, Gene was a physical chemist in the Benches 
Branch of the Western Electric plant in Hawthorne, Illinois. His originality and 
inventiveness immediately became evident and was rewarded by the issuance of 
two U. S. patents. 

Gene joined the Geophysical Research Corporation in 1925 and served as 
Field Superintendent until 1932. He was the first man in the field for the GRC, 
and in 1926 his refraction crew discovered the Moss Bluff salt dome in Liberty 
County, Texas, the first salt dome discovery made in this country by an Ameri- 
can crew. Next, his crew discovered Port Barre and Bullard; later, Gene estab- 
lished a record of nine salt dome discoveries in twelve months. 

After resigning from GRC in 1932, Gene assisted in the organization of In- 
dependent Exploration Company; he served as its president and general manager 
until 1937, when he severed his connection to organize Subterrex, the first 
American company to offer geochemical parties on a contract basis. From 1942 
until his death he was a senior partner in Geochemical Surveys of Dallas, Texas. 

The study of the strategy and tactics of war quickly gained Gene’s interest 
and attention. He began with “On War” by von Clausewitz, ‘“‘Napoleon” and 
“Great Captains” by Dodge, and Mahan’s “The Influence of Sea Power upon 
History.”’ Gradually, he acquired a large collection of books which covered the 
lives and campaigns of the great captains in detail. During the middle thirties, 


he began to recognize the similarities between the successful conduct of war and 
success in prospecting for petroleum; beginning with a paper in 1935, “On the 
Strategy and Tactics of Exploration for Petroleum,” he delineated elements in 


99 


exploration which were analogous to “Offense,” ‘Defense,’ “Surprise,’’ and 
“Economy of Force.”” The Second World War intensified his interest; following 
the third paper on strategy and tactics in 1939, he published ‘Prospecting 
Principles” for private distribution in 1940; in this paper he established explora- 
tion on a broad basis and exhibited striking parallelisms between various aspects 
of prospecting and the conduct of the initial campaigns of World War II. 

A keen intellect coupled with high intelligence roams the field of knowledge 
and pauses here and there to cultivate small areas. Thus, in the early forties, 
Gene became interested in the study of society and devoted intensive attention 
to the works of Pareto, Michels, Mosca, Sorel, Spengler, Sorokin, Toynbee and 
Burnham. A paper on the application of the scientific method in the social field 
followed this study; calling attention to Pareto, he defined logical conduct as 
human action in which three components simultaneously apply: 1. The action is 
motivated by a deliberately held goal or purpose; 2. The goal is attainable; 3. 
The means or steps taken to reach the goal are appropriate to achieve that end. 
Numerous examples of non-logical conduct were given; Gene had an abundance 
of experience with scientists who were unable to apply the scientific method even 
in their chosen area of science, and he did not anticipate the early adoption of 
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the scientific method in the field of social activity, since the use of the scientific 
method would require discarding emotion in human conduct. 

Gene believed that scientists and engineers should participate in civic affairs. 
Accordingly, he assumed civic responsibilities. He joined several young lawyers 
who had returned from the war, and the reorganized the Taxpayers Associa- 
tion of Houston into the Tax Research Association of Houston and Harris 
County. 

During the late war years, Gene became active in the National Society of 
Professional Engineers and its state and local groups. He served as a national 
director for one year, in 1946, as a state director, and for several years he was 
active in the local San Jacinto Chapter as president, 1947-1948, and as program 
chairman. The TSPE awarded to Gene a certificate for outstanding service in 
1958. 

When the Engineers’ Council of Houston was organized, Gene took an ac- 
tive part and became its vice-chairman in 1946-1947; he helped in writing its 
constitution and in the preparations for two symposia, one on atomic energy and 
the other on natural resources, which were addressed by illustrious speakers with 
national reputations. Finally, he was a director of the Houston Engineers Club 
and served as the editor of its monthly publication, The Slide Rule, during 1947- 
1948. 

In Dallas, Gene was a member of the Northwood Country Club and of the 
Hillcrest Estates Association; he was president of the latter organization in 1949, 

Early in the year 1930, Gene joined the late D. C. Barton and others in a 
drive to persuade the geophysisicsts residing in the Houston area to form an 
organization. After several luncheon meetings, the drive was successful and the 
Society of Petroleum Geophysicists was organized with Barton as President and 
Gene as Vice-President. The new Society became a technical division of the 
American Association of Petroleum Geologists in 1932; the AAPG had published 
the geophysical papers presented at the San Antonio meeting in March, 1931. 
The American Physical Society had published the geophysical papers offered at 
its meeting in New Orleans in December, 1931. Thereafter the AAPG published 
the papers presented at the meetings of the Society of Petroleum Geophysicists 
in its Bulletin as the Transactions of the Society. The last publication by the 
AAPG was Volume V, which appeared in March, 1935. The AAPG determined 
editorial policy and frowned upon the printing of highly technical and mathe- 
matical papers; however, the objections were usually removed by correspondence 
and meritorious papers were published. 

The growth of the Society was slow; two years were required to double the 
membership. Then, the membership doubled in 1934 and trebled by March, 
1937. The increase in 1934 was due largely, if not entirely, to Gene’s efforts; he 
had been elected President in that year. Soon he recognized some impediments 
in the way of growth; he became convinced that the Society must become in- 
dependent, must have its own publication, and that its activities must include 
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the instructors of the subject of geophysics in the institutions of learning and 
the mining geophysicists; in brief, the scope must embrace all geophysical activity 
and membership must be open to all qualified workers in geophysics in the broad- 
est sense. At about this time, it was reported that the AAPG was considering 
an amendment to its constitution to provide that all members of technical divi- 
sions must be members of the AAPG, but it was never established which level of 
the AAPG did the considering, and no official communication on the subject 
was received. 

More and more members of the Society adopted the viewpoint that the name 
must be changed, a new constitution must be written, and independence from the 
AAPG must be achieved. There were protests by AAPG members and predictions 
of dire consequences if the Society were to leave the fold; there were advices from 
Society members with geological background to leave well enough alone. The 
rebels prevailed; in 1936 the name was changed to The Society of Exploration 
Geophysics and, in 1937, the AAPG released the Society. The membership then 
increased spontaneously by leaps and bounds without any concerted effort on the 
part of the officers. For years upon years, the SEG and the AAPG held their 
annual meetings together, at the same time and place, and most cordial relations 
existed between them. As for Gene, his objective having been achieved, he never 
again made any effort to influence the operation or the policies of the SEG; its 
management and operation were taken over by devoted and competent younger 
men serving the Society. The success of the Society attests to the effectiveness of 
their efforts. In 1951, the Houston Geophysical Society elected Gene to honorary 
membership. 

The Texas Acacemy of Science, the American Assocition for the Advance- 
ment of Science, and the Institute of Petroleum Technologists (British) elected 
Gene to fellowship. He was a member of the AAPG, the AIME, the American 
Geophysical Union, the American Society of Military Engineers, the U. S. 
Naval Institute, the Society of the Sigma Xi, the National Society of Professional 
Engineers, the Foreign Policy Association, and, of course, the SEG. 

Gene first touched my life in 1929; during the ensuing years the acqua’ 
ship ripened into friendship. He had a pleasant, stimulating personalit,g 
tility of interests, and a marked talent for analysis. “Let’s take it a 
examine the pieces,’’ he would say. Whatever captured his interest, studied, 
organized, and mastered. An evening with him was always a stimulag- ug experi- 
ence. We could discuss and debate, agree and disagree, but our borg.’ of friend- 
ship were strengthened. After his removal to Dallas, we were fo/j/d to com- 
municate by correspondence and this we did, at least once each morph. 

Gene married Johnie Laureme Ransone on June 24, 1928. His gridow, daugh- 
ter Carol Jean (Mrs. Ronald F. Brown), grandchildren Carol gfe «nd Ronald 
Trent, sister Pauline of Park Ridge, Ill., and brothers Carol G.4@ Mattoon, IIl., 
and Forrest C. of Burbank, Calif., survive him. He was a kind gd loving father 
and husband with a highly developed sense of family respon #ility. His family 
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is joined in mourning for him by a host of friends who have lost a good friend of 
sterling character with a keen sense of right and wrong, who never failed to take 
the right course. The art of exploration for minerals in the United States has 
lost one of its pioneers who had acquired intensive training in several branches 
of prospecting, established the broad principles of general, ‘environmental’ 


exploration, and retained his equanimity in the face of adversity. 
L. W. BLau 
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CARL MARION WERT 


Mr. Carl Wert, age 48, underwent surgery in an Oklahoma City hospital on 
July 4th, but did not recover. He passed away July 18th. Funeral services were 
held in Henryetta, Oklahoma, July 21st at the First Methodist Church, with 
burial in the Westlawn cemetary. 

Carl was a veteran employee of the Gulf Oil Corporation, and a member of 
the Society of Exploration Geophysicists since 1940. He was born in Mexico, 
Juniata County, Pennsylvania, on March 28, 1912. He attended the Westing- 
house High School and Carnegie Institute of Technology, both of Pittsburgh, 
Pennsylvania. Carl joined Gulf Research & Development Company in June 1929 
as an assistant in the Gravity Department. Approximately 1932, he transferred 
to seismic work and traveled over the Midcontinent states, receiving promotions 
as Gulf’s seismic work expanded. In 1947, as party chief, he embarked on a 
3-year foreign assignment with a crew which worked in Nicaragua and Cuba. 
While living in Fort Worth, Texas, from 1950 to 1957, he served Gulf in the 
capacity of geophysical supervisor. During this period, he was elected to the 
office of Secretary-Treasurer by the Fort Worth Geophysical Society for a one- 
year term. He was a member of the North Texas Geological Society. Mr. Wert 
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transferred to Wichita Falls, Texas, in November 1958 as Staff Geophysicist, the 
position he held at the time of his death. 

As an actively dedicated 32nd Degree Mason, Carl was a member of the 
Masonic Blue Lodge, and of the Valley of Guthrie Consistory. Also, he was active 
in the Methodist Church, serving as Sunday School teacher and President of 
the Men’s Club. 

Carl and Edna Turk, of Maud, Oklahoma, were married December 23, 1934. 
Typical of our migratory profession, their first son was born in Oklahoma, and 
the second in Texas. The older son, Peter, is married and serving in the U. S. 
Navy with an NROTC Commission. 

The second son, Turk, was killed in a hunting accident in February 1959. 
While a friend was unloading a gun, it accidentally discharged and ended the 
career of a very promising boy at the age of 15. Turk was an outstanding athlete, 
and in addition had received many awards as an honor student. The many 
friends of the Wert family assisted in establishing a Turk Wert memorial at 
Stripling Jr. High School, Fort Worth, Texas. It consists of a plaque on which is 
engraved annually the name of the student who does the most with his inherent 
abilities. The student receiving this award also receives a replica of the plaque 
which hangs in the school library. 

Other survivors besides his wife and son include Carl’s father, S. V. Wert, 
Columbus, Ohio; three brothers, Jerry Wert, Thompsontown, Pennsylvania; 
Wilbur Wert, Mifflintown, Pennsylvania; Dane Wert, Harrisburg, Pennsyl- 
vania; two sisters, Mrs. Faye McCormick, Detroit, Michigan; and Mrs. Carol 
Jean O’Conner, Germany. 

Those who knew Carl remember him as a true friend. He always left people 
with smiles on their faces, and feeling happy because of the contact with him. 
With his untimely passing, our profession losses a member dedicated to its appli- 
cation. 

W. B. RoBinson 
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LETTERS TO THE EDITOR 


THE EFFECT OF TEMPERATURE ON THE RESPONSE OF 
SCINTILLATION COUNTERS* 


A. F. GREGORY AND A. L. FOURNIER 


It has been known for some years that the output of scintillation crystal- 
photomultiplier tube assemblies is temperature dependent (Kinard, 1957; Connor 
and Husain, 1960). Under certain circumstances, the counting rate may vary 
with the temperature but this response variation is also dependent upon the 
spectrum of pulses being investigated and upon the minimum pulse size which 
can be counted. 

Temperature effects have not generally been considered in geophysical inves- 
tigations with scintillation counters. Recently, a significant variation in counting 
rate with temperature (—1 to —4 percent per degree centigrade) was observed 
by the writers during a field investigation of the gamma radioactivity of natural 
mineral sources. This effect is probably negligible for ordinary prospecting pur- 
poses but may be important in certain quantitative aspects of ground, airborne 
and well-logging surveys in which the instrument temperature is not controlled. 

In the work reported here, the variation of counting rate with temperature 
has been determined for four different detectors with a pitchblende source. Brief 
descriptions of the four instruments and their observed response-temperature 
coefficients are given in the accompanying table. As there may be large differ- 
ences in temperature sensitivity between two photomultiplier tubes of the same 
type (Kinard, 1957), the response-temperature coefficients probably will be 
peculiar to the crystal-photomultiplier assembly and the spectrum for which it 
was determined. 


TABLE 
Response- Observed 
Crystal multipli Detector Temperature Tempera- Source 
ultiplier Coefficient ture Range 
3 inch x14 inch diam, | Dumont | 6292 EA 136S airborne scintil- —4.0% °C. 56-72°F. 90 Ibs. ‘pitchblende (ap- 
Nal(TI) lation counter —2.2%/°F. proximately 15% 
at 5 feet and 42 feet 
from crystal 
3 inch X3 inch diam. Dumont 6363 Detecting element pre- —2.3%/°C. 20-25°C. ” vial of pulverized pite h- 
Nal(T)) pared for use with a 100 —1.3%/F. blende at 24 inches from 
channel kicksorter, C.D C. crystal 
model A.E.P.-2230 
4 inch x2 inch diam. RCA 5819 M.E.L. type A.E.P. 1903R —1.3%/°C. 66-86°F. 90 Ibs. pitchblende (ap- 
Nal(Tl) Mark II airborne scintil- —0.7%/°F. proximately 15% UxOs) 
lation counter at 174 feet from. crystal 
1 inch X1 inch diam. RCA 6199 PRI model 111 probably small, bet 61-66°F. “90 Ibs. pitchblende (ap- 
Nal(TI) data not definitive. proximately 15% Ui0s) 
Possibly 0.5 to at 5 feet from crystal 
1.0% /°C. 


* Published with permission of the Director, Geological Survey of Canada. 
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The EA-136S and PRI-111 detectors are self-contained, portable units. For 
~ these instruments, the counting rate variation was measured for different ambient 
air temperatures. The M.E.L. unit has crystal, photomultiplier and preamplifier 
enclosed in a light-tight box comprising $-inch aluminum sheeting with }-inch 
aluminum cover. In this instance, the temperature is that registered by a ther- 
mometer bulb located inside the box at the crystal end of the assembly and about 
; inch from the end of the box. The photomultiplier and crystal used with the 
C.D.C. kicksorter were housed in a box of thin aluminum and pine sheeting. This 
housing restricted but did not prohibit ventilation during the tests. The pre- 
amplifier was located in an adjacent but separate ventilated compartment. The 

thermometer bulb was placed next to the crystal container. 
The variation of counting rate with temperature for the M.E.L. unit (Figure 

1) is typical of the variations for the other instruments. 
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Fic. 1. Typical graph showing variation of response with temperature (June/59). 
Detector: M.E.L., type AEP 1903R, Mk II., #2 counter only. 4” X2” diam. NalI(TI) crystal, 


RCA 5819 photomultiplier. 
Source: 90 Ibs. pitchblende (approximately 15% U,sOs) at 174 feet from crystal. 
Response-temperature coefficient: —0.7 percent per degree fahrenheit. 
—1.3 percent per degree centrigrade. 


In measuring response-temperature coefficients, care should be taken to con- 
firm that the measured temperature is that of the crystal-photomultiplier assem- 
bly. Apparent divergency from the normal variation may result from the heat 
of direct sunlight falling on the metal housing of the detector. Such an effect was 
observed for the M.E.L. detector (see Figure 1). It is probable that the housing 
absorbed heat and reradiated it to the nearby thermometer. The crystal and 
photomultiplier in the center of the insulating air mass app«rently did not warm 
up to the same extent. Thus, upon replacement of the cardboard shade protect- 
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ing the housing from direct sunlight, the divergent part of the curve returned to 
follow the normal response-temperature variation. Similar effects were observed 
with the other detectors. 

Net counting rates were also compared with air density, background counting 
rate and time but no correlation was found. 

The observed variation of counting rate with temperature of. the scintilla- 
tion crystal-photomultiplier tube assembly may be important in ground, air- 
borne and well-logging surveys in which the radiometric data are contoured or 
compared as to magnitude, especially with high counting rates. For ordinary 
prospecting, in which radioactivity maxima are recognized relative to the back- 
ground of adjacent rocks, these temperature effects are not likely to be significant. 

REFERENCES 
Connor, R. D. and Husain, M. K., 1960. Some factors affecting the gain and resolution of photo- 
multiplier tubes: Nuclear Instruments and Methods, v. 6, no. 4, p. 341. 


Kinard, F. E., 1957, Temperature dependence of photomultiplier gain: Nucleonics, v. 15, no. 4, pp. 
92-97. 


ON THE CALCULATION OF AFMAG* ANOMALIES 


JAMES R. WAITT 


It has been proposed by S. H. Ward (1959) and C. S. Wright et al. (1960) 
that nonuniformity in the earth’s crust will tend to tilt the direction of the fluc- 
tuating components of the geomagnetic field from its normal horizontal direc- 
tion. In fact, Ward and colleagues have developed a practical scheme (AF MAG) 
for geophysical prospecting which makes use of this phenomenon in the frequency 
range from 100 to 400 c/s. Sir Charles Wright et al., on the other hand, have 
found that at much lower frequencies, .005 to 1 c/s, a substantial vertical com- 
ponent of the field in the vicinity of coast lines or other situations where the 
earth’s crust departs from horizontal stratification. 

It is the purpose of this short note to show how these anomalies may be cal- 
culated for certain idealized situations. Actually, the problem is closely related 
to those in electromagnetic prospecting using artificial fields. 

At least for frequencies in the ELF (extremely low frequencies, meaning the 
range 1 to 3.00 c/s), the origin of the fields are from lightning discharges (Pierce 
1960). In fact Holzer and others (1957) have shown that there is a high correla- 
tion between the world-wide thunderstorm activity and the integrated level of 
the noise in this frequency range. For frequencies below 1 c/s the origin of the 
fields is not yet understood. It is probable that they are related to current sys- 


* AFMAG—Audio Frequency Magnetics. 
t Newmont Exploration Ltd., Danbury, Connecticut. 
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tems at ionospheric heights which themselves are induced by solar particles 
(Campbell 1960). In any case, at the surface of a homogeneous smooth earth the 
fluctuating component would be horizontal or approximately so (Wait 1960a). 

For purposes of illustration, the earth is assumed to be locally flat, thus it 
can be represented as a homogeneous half space, of conductivity 7) and permea- 
bility wo, defined by y<0.' The surface of the earth is the xz plane. Within the 
earth and running parallel to the z axis is an infinite cylinder of radius a, of 
conductivity o and permeability u. The axis of the cylinder is located at y= —b 
in the piane x=0. In view of the nature of the problem, it is only the perpen- 
dicular or x component of the exciting magnetic field Ho which needs to be con- 
sidered. Furthermore, because of the assumed uniformity of the exciting field 
Ho, at the surface, it can be assumed to decay exponentially with depth. Thus 
the field which acts on the cylinder may be taken as Ho, exp (—~vo0b) where 
Yo = (ioouow)'’*. Such a field will induce a line dipole at y= —b. Higher-order poles 
can be neglected if | ya| <1 or if ab. It will be assumed that such a condition 
is met. The strength of this induced line dipole can be readily computed by mak- 
ing use of known results (Wait 1960b). However, obtaining the resultant fields 
of the line dipoles at the surface is not so straight-forward since the boundary 
conditions at the air-ground interface (i.e., y=0) must be met. In principle, the 


1 See Fig. 1. 
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solution can be obtained but it is not of much use until some of the integrals 
can be evaluated. (This problem is now being investigated.) As a rough approxi- 
mation, it is assumed that the resultant fields are those derived by magneto- 
statics (strictly valid if yo=0) but then multiplied by the correction factor 
exp (—vop) where p=(b?+.27)'” is the distance from the axis of the cylinder to 
the observer at the surface. 

A straight-forward derivation yields the following for the surface fields 


r| | [—-yo(x? + 6)! b 
— = — ———— | exp [—yo(x? + b?) vo 
(x? + 82)? I Y Ye J 


| exp + — yobd], 


+ 
and X =(opuw)'? a, J; is the modified Bessel function of order one and the prime 
indicates a derivative with respect te its argument. 
Numerical values of the complex factor T; have been given previously (Wait 
1960b). For a highly conducting body, such that Y¥<1, 7; approaches +1. For 
X<1 it is readily shown that 


which could be obtained directly from magneto-statics. For intermediate values 
of X the imaginary part of 7, is finite which means that there is a phase shift 
between the H, and H,. In any case, it is seen that right over the cylinder, the 
magnitude of H, has a maximum. It is given by 

H, | T,a° 


2b? P [= 


For low frequencies, such that (aouow)'/? 61, and if the depth 6 is, say, 3a, 
the vertical component would never be more than 6 percent of the undisturbed 
horizontal field. 

This approach may be used for other models such as a buried wedge or knife 
edge. While the method is certainly not rigorous it should be useful for obtaining 
first-order estimates. In the case of two-dimensional bodies it is necessary to 
consider the induced electric dipoles. These would be proportional to the electric 
(or telluric) field. 
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Hy, 2xba 
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Asociacion V. 1 de Geofisi 
(20), Caracas, Venezuela, char- 
tered January 5, 1959, officers 
elected in March. 

E. J. Assiter, pres. 

D. W. Danz, 2nd v-pres. 
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Colorado School of Mines Society 

of Student Geophysicists 

Joseph R. Anzman, secty. Depart- 
ment of Geophysics, Colorado 
School of Mines, Golden, Colo- 
rado. 

Meetings: Monthly, 2nd Monday, 
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Geophysical Society of Saint Louis 
University 
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Norman J. Guinzy, secty. 410 
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Meetings: Monthly, 2nd Wednes- 
day 7:30 P.M., meeting only, 
Institute of Technology. 
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University of Toronto Geophysical 
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B. R. Krause, secty. 502 Gilbert 
Ave., Apt. 107, Toronto, Ont., 
Canada 
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1515 Massachusetts Avenue 
Washington 5, D. C. 
Paul E. Klopsteg, President 
Paul A. Scherer, Treasurer 


Dael Wolfle, Executive Officer 


The American Association of Petroleum 


Geologists 


A Cooperative relationship 
Box 979 
Tulsa 1, Oklahoma 
.. G. Weeks, President 
A. H. Bell, Vice-President 


O. Koefoed, Editor 


L. Solaini, Past President 


National Research Council, Division of 
Earth Sciences 


SEG is a member society 
2101 Constitution Ave., N.W. 


Washington 25, D. C. 


John N. Adkins, Chairman 
Harry H. Hess, Past Chairman 
Edward B. Espenshade, Jr., Chairman- 


Designate 


William R. Thurston, Executive Secretary 


Petroleum Exploration Society of Libya 


H. T. Morley, Secretary-Treasurer 


G. E. Murray, Editor 


G. S. Buchanan, Past President 
R. H. Dott, Executive Director 


European Association of Exploration 


Geophysicists 


Mutually affiliated 
30, Carel Van Bylandtlaan 
The Hague, Netherlands 


Mutually A ffiliated 


P. O. Box 900 
Tripoli, Libya 


R. E. Doan, President 


P. F. Burollet, Jst Vice-President 
J. F. Romminger, 2nd Vice-President 


J. C. Wilson, Secretary 

P. Arni, Treasurer 

J. S. Cullison, Director 

R. S. M. Templeton, Director 
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MEMBERSHIP APPLICATIONS RECEIVED 


Applications for active membership have been received from the following candidates. This pub- 
lication does not constitute an election but places the names before the membership at large, in ac- 
cordance with Bylaws, Article III, Section 4. References are listed in parentheses following the names 
of each candidate. If any member has information bearing on the qualifications of these candidates 
he should send it to the president within thirty days. 

APPLICATIONS FOR ACTIVE MEMBERSHIP 
E. R. Critchett (W. J. Ricker, Jr.; N. Vander Sleen, E. E. Smith) 
Gilbert, Freeman (R. Bowman, M. K. Smith, M. M. Backus, E. J. Stulkew) 
John E. McCall (S. C. Stoneham, W. L. Besham, T. L. Slaven, R. E. Doe) 
Kenneth J. Parry (H. R. McCormack, J. M. Desmond, R. J. Koenig, J. Swistek) 
James H. Pennington (Chester Sappington, R. J. Hope, J. D. Laker, J. S. Spencer) 


APPLICATIONS FOR TRANSFER TO ACTIVE MEMBERSHIP 


Donald C. Bowman (J. F. Johnson, A. C. Reid, H. J. Jones, D. R. Oksa) 

Daniel F. Brennan (C. M. Moore, Jr., K. E. Burg, R. C. Dunlap, Jr., Neil W. Mann) 
James W. Buelow (L. Y. Faust, A. I. Innes, K. M. Lawrence, F. F. Campbell) 
Perry E. Byerly (W. H. Diment, R. M. Hazelwood, Rudolph Black, H. R. Joesting) 
Frank Clynch (J. M. Crawford, Frank Searcy, L. E. Whitehead, W. E. N. Doty) 
Gunter Wm. Conrad (C. N. Page, C. J. Lomax, A. E. McKay, W. A. Seal, Jr.) 
Theodore R. Dunn (S. B. Stewart, D. G. Coppedge, D. B. Darden, S. A. Hatcher) 
John O. Erickson (J. W. Keffler, Earl M. Wolters, H. L. Mendenhall, W. H. Courtier) 
Morris F. Gillett (E. L. DeLoach, C. H. Hightower, M. S. Hathaway, J. P. Woods) 
John C. Healy (J. Athanasion, Richard L. Fentem, Wm. M. Krum, J. B. Nichols) 
Herr, Wilfred C. (E. B. Miller, C. H. Gullick, L. E. Whitehead, L. K. King) 
Thomas C. Holt (T. O. Hall, M. H. Collum, R. J. Watson, Dale Turner) 

Jason S. Johnson (M. G. Frey, Alton McClung, E. J. Grivetti, E. J. Northwood) 
James B. Jordon (F. J. DiGuilio, B. B. Strange, N. P. Cramer, Dean Walling) 
Ernest A. Kaarsberg (Warren Westphal, J. T. Wilson, T. C. Poulter) 

John W. Keeler (S. R. Marsh, E. J. Laurent, J. D. Sollenberger, J. A. Dunn) 

Wm. E. Kelley (H. R. Hockoday, W. K. Rickel, J. T. Gilliam, R. C. Roesch) 
Robert F. Kraye (Lee B. Park, Robert J. Knott, Will Alex Grosse) 

Laidlaw, Wm. H. (R. J. Yuzbick, H. C. McCarver, L. V. S. Roos) 

Sam R. Lake (W. T. Lea, R. W. Buhrman, John W. Bell, Jr., Glenn E. Hunter) 
Lawyer, Leon C. (J. K. Pawley, R. B. Wilcoxen, J. C. Kern) 

Harold R. McCormack (R. J. Koenig, J. Swistek, J. A. MacMillan, J. J. O’Bertos) 
Donald K. MclIvor (J. B. Murphy, D. C. Skeels, C. J. Chapman, J. T. Hughson) 
Montgomery, John A., Jr. (J. H. Wright, Jr., Miles T. Rader, W. J. Pfeffer, G. B. Neville) 
John R. Muir (F. A. Hale, W. A. Robertson, G. C. Alvey, H. Agnew) 

John J. O’Bertos (H. R. McCormack, R. J. Koenig, J. MacMillan, J. Swistek) 
Chas. G. Oldham (F. G. Blake, N. A. Riley, G. D. Garland, L. P. Stephenson) 
Glenn Robertson (G. H. Samuels, O. G. Holekamp, C. H. Gregory, F. W. Oudt) 
Donal E. Ruminer (W. H. Roane, W. H. Courtier, F. H. Morris, J. M. Sanders) 
Henry B. Sawatzky (D. H. Parry, M. C. Woodard, A. Armstrong, R. H. Carter) 
James Wm. Schneider (C. B. Smith, S. B. Stewart, C. C. Sellers) 

Edward Thiel (J. P. Woollard, A. P. Crary, K. L. Cook, W. E. Bonini) 

Howard F. Thompson (C. Newton Page, Robert J. Knott, Wm. A. Seal, Jr.) 

Wm. M. Walton, Jr. (J. H. Wright, Jr., M. T. Rader, W. J. Pfeffer, D. C. Meyer) 
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ANNOUNCEMENTS 
SPECIAL ISSUE ON ENGINEERING GEOPHYSICS 

Plans are under way for a special issue of GeopHysics devoted to engineering geophysics. The 
editor has appointed Dr. Charles L. Drake, Lamont Geological Observatory, Palisades, New York, 
Chairman of a special editorial committee for the volume. Tentatively the December, 1961, issue 
has been chosen to feature engineering geophysics papers. 

Members of the society are invited to submit papers on engineering geophysics to Dr. Drake, 
whose committee will have the responsibility of reviewing all submitted papers for the special issue 
before their submittal to the editor. 

Dr. Drake has appointed Howard Pincus, Ohio State University, Columbus, to serve on the 
committee with him. Other members of the committee will be appointed soon. 


GEOPHYSICS ON MICROFILM 


Volume 21 through 24 of Geopuysics, 1956 through 1959, are now available, one reel per volume 
as follows: 
Volume 21 (1956) $4.45 
Volume 22 (1957) 3.95 
Volume 23 (1958) 4.60 
Volume 24 (1959) 4.60 

Complete with tables of contents, indexes, and advertising, these microfilmed volumes may be 
purchased direct from University Microfilms, Ann Arbor, Michigan. 


GEOPHYSICAL SOCIETY ESTABLISHES BEST PAPER AWARD 


The Canadian Society of Exploration Geophysicists announces the establishment of an annual 


“Best Paper’? award. One hundred dollars will be awarded for the best paper presented during the 
year at one of the several meetings sponsored by the Society. The author must be a resident of the 


Calgary area, but need not be a member of the Society. 

Judges, four in number, include R. B. Ross of British American Oil Company “td., D. C. Skeels 
of Imperial Oil Limited, A. Junger of Shell Oil Co. of Canada Ltd., and G. J. Blundun of Home Oil 
Co. Ltd. 

Funds for the awards will be provided by the publishers of Canadian Oil and Gas Industries. 


A JOURNAL OF THE AMERICAN INSTITUTE OF PHYSICS IS NOW A MONTHLY 
The Physics of Fluids, a journal of the American Institute of Physics, will be published monthly 
instead of bimonthly, effective with Vol. IV (1961). Inquiries should be sent to: American Institute 
of Physics, 335 East 45 Street, New York 17, N. Y. 
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PERSONAL ITEMS 


Harry W. Cartson, formerly with the California Company, is now a consultant in seismic in- 
terpretation. His new address is 726 Guaranty Bank Bldg., Denver, Colorado. 

Ep Fritz, with Pure Oil Company, has been transferred to Houston as Chief Geophysicist, 
Southern Division 

Jackson W. Brown, formerly Domestic Manpower Administrator for GSI, is now College Re- 
lations Coordinator with Texas Instruments, Inc. 

Ernest H. HumgE, Senior Seismic Computer with Humble, is on loan to Esso-Republic, Bordeaux, 
France. His assignment is geophysical interpretation of Esso-Republic permits. 

RoBerRT H. MANSFIELD is now Supervisory Geophysicist, USAF (AFTAC), at: HQ-USAF 
(AFTAC), Washington 25, D. C. 

Epwarp D. Lynton has a complete set of GEopHysics from Vol. I, No. 1 (Jan. 1936) to Vol. 
XX, No. 4, (Oct. 1955), including the bound Cumulative Index 1931 through 1953. He has priced 
them at $1.50 per copy but will consider any reasonable offer. They will be sold only as a set, not as 
individual copies. His address is 422 Beloit Ave., Berkeley 8, California. 

RoBERT VAN NosTRAND’s address is 23 Rue Galilée, Paris 16, France, not as erroneously quoted 
on page 725 of the June issue of GEopHysics. 

Dr. JULIAN F. Evans, who holds his Ph.D. degree from Rice University and has been employed 
by the Pan American Petroleum Corp. research division for the past 19 years, joined the University 
of Tulsa faculty Thursday as an assistant professor of mathematics. 

A native of Michigan, he received a bachelor of arts degree in English from the University of 
Oklahoma in 1926 and a master of science degree in physics in 1932. 

He received his Ph.D. degree in physics from Rice in 1939. From 1940-41 he taught physics, 
geophysics and mathematics at Texas Christian University. Dr. Evans then joined the Pan American 
geophysical research department. 

He is a member of Phi Beta Kappa and Sigma Xi honor societies and the Society of Exploration 
Geophysicists. 

Western Geophysical Company of America, a subsidiary of Litton Industries, has moved into 
its newly completed four-story headquarters building adjacent to its laboratories and shops at 933 
North La Brea Ave., Los Angeles, California. 

CarLTON D. SPEED, transportation chairman for the Houston Geological Society, headed a group 
of 26 Houston Geologists and their families who toured Western Europe before attending the Inter- 
national Geological Congress in Copenhagen on August 14. 

Epwin B. NeITzeEL, former Technical Supervisor of Atlantic Refining’s Geophysical Laboratory, 
has accepted a position as Engineering Group Head in the Collins Radio R and D IV Department. 
He will be working on computer data processing problems. 

Doy Le D. Mize, vice-president and Houston manager of Mandrel Industries, has been named 
executive vice-president in charge of all Mandrel operations. 

The announcement was made following the first meeting of the Board of Directors of Mandrel 
Industries held in Houston. 

Mr. Mize was appointed operations manager of Electro-Technical Laboratories in 1951. When 
Mandrel purchased Electro-Tech in 1956, Mr. Mize was appointed a vice-president. 

In 1958 Mr. Mize became head of the Electric Sorting Machine division, assuming responsibility 
for the entire Houston operations. 

In 1959, he was elected to the Mandrel Board of Directors. 

Prior to his association with Mandrel, Mr. Mize was technical and field operation supervisor 
for a leading geophysical company. 

Mr. Mize attended the University of Houston and majored in chemistry. He is a member of both 
the Society of Exploration Geophysicsts and the Houston section. 

A native of Iola, Texas, Mr. Mize, his wife and six children reside at 807 Chowning Rd., Houston. 
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PERSONAL ITEMS 1303 


DanteEL Lum, formerly geophysicist with the South Dakota Geological Survey, recently com- 
pleted gravity surveys on the islands of Bohol, Cebu, and Mindoro, Philippines, under contract 
assignment to Atha Exploration Co. of Los Angeles, California. Mr. Lum is presently residing in 
Honolulu, Hawaii. 

Peter C, Sunpt will be chief administrator and technical director of Mandrel’s two Houston di- 
visions, Electro-Technical Labs and Electric Sorting Machine Company. 

Mr. Sundt, who graduated from Lamar High School before entering Rice, began his career as an 
engineer with Hughes Aircraft Company. He came to Electro-Technical Labs in 1953 as a design 
engineer. 

He has since been chief engineer and general manager of Electro-Tech France with headquarters 
in Paris. 

He was called to Houston from France to accept the new position. 

Mr. Sundt is a member of the Houston Geological Society, Institute of Radio Engineers, and 
the Society of Exploration Geophysicists. 

He and Mrs. Sundt and their three children will make their home in Houston. 

KENNETH S. TOWNSEND has resigned his position as Station Manager for Birdwell Division of 
Seismograph Service Corporation to take a position as Geophysicist with Sunray Mid-Continent 
Oil Company in Tulsa. 

EvGENE McDermott, Chairman of the Executive Committee of Texas Instruments Incorporated 
in Dallas, Texas, has been named to the Stevens Institute of Technology Board of Trustees. 

Mr. McDermott is a 1919 Stevens graduate. He entered the oil exploration field in 1925 when he 
helped found the Geophysical Research Corporation. 

In 1930 he co-founded Geophysical Service Incorporated, thus sharing in the development of 
geophysical exploration techniques still considered the most reliable and economic means of locating 
oil deposits in the earth. 

Mr. McDermott served as President of St. Mark’s School for Boys in Dallas from 1945 to 1955 
and continues as a Trustee. He is Co-founder and Trustee of the Graduate Research Center of Dallas, 
member of the Executive Committee and Board of Trustees of the Texas Research Foundation, and 
a Trustee of the Hockaday School for Girls in Dallas. 

His other memberships include the Southwestern Medical Foundation, and the Visititing Com- 
mittee on Physics of the Massachusetts Institute of Technology. He also serves on such boards as the 
Dallas Art Museum and the Dallas Theater Center. Mr. McDermott is a past president of the 


Society for Exploration Geophysicists. 

Mr. McDermott is one of twenty-five Stevens trustees who serve five year terms. Three other 
trustees elected from the alumni body serve for three years. ' 

BERNARDO F. GROSSLING, formerly with California Research Corp., has joined the Inter-Ameri- 
can Development Bank, Washington, D. C. His present position is Acting Chief, Engineering Di- 
vision. 


FRANK GOLDSTONE, erroneously reported in August GeopHysics as having entered into partner- 
ship with Geoprosco of London, is in fact representing Petty Geophysical Engineering Company in 
their London office. 

Professor ENpERS A. RoBinson of the Mathematics Department of the University of Wisconsin 
will spend the academic year 1960-1961 at Statistiska Institutionen vid Uppsala Universitet, Jarnbro- 
gatan 18, Uppsala, Sweden. 

J. Brian Esy, an independent consultant geologist in Houston, Texas since 1931, is now asso- 
ciated with Bradco Oil & Gas Co. Mr. Eby received his A.B. and Ph.D. degrees from Johns Hopkins 
University. Before going into consulting work, he served as a geologist with the U. S. Geological 
Survey and with Shell Petroleum Corp. Dr. Eby has lectured at Johns Hopkins, Princeton, University 
of Texas, Texas A & M College and other colleges and universities. He is the author of numerous 
technical articles on geology, and has traveled in Europe, Africa, Canada and Mexico as foreign 
representative for U. S. oil publications. Dr. Eby’s professional memberships include AAPG, IPAA, 
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American Association for Advancement of Science, SEG; Texas Academy of Science, and Royal 
Society of Arts, London, England. 

Opie M. Hotmes (1937 to present), formerly with The Petty Geophysical Engineering Company, 
San Antonio, Texas, Rogers-Ray Geophysical Co. of Houston, Texas, and Seismic Analysis, Inc., 
Dallas, Texas, has recently gone into private practice as an Exploration Consultant, 4041 Rosa Road, 
Dallas 20, Texas. 

Cot. GeorGE B. Kaiser USAF (Ret.), has recently been appointed Manager, Planning Programs 
for Robert O. Vaughan, the Vice-President of Marketing at Southwestern Industrial Electronics. 

Col. Kaiser is a graduate of Texas A & M with a Masters Degree in Electrical Engineering and 
is a registered Professional Engineer in Texas. He is also an active member of the Institute of Radio 
Engineers, Society of Exploration Geophysicists, and American Association of Petroleum Geologists. 

Index Geophysical Surveys Corp., 2707 Milam Street, Houston, Texas, announces the recent 
addition of NEAL CLAyTON to the staff as a Geophysical Supervisor. For the past twelve years Mr. 
Clayton was District Supervisor for Sohio Petroleum Company in Lafayette, Louisiana. Mr. Clayton, 
a graduate of L.S.U. with a Masters in Physics and Mathematics, has published a number of articles 
in World Oil and Geopuysics during the past twelve years dealing with exploration problems. The 
new Supervisor is a native of Mississippi, and has had extensive experience in most of the active sedi- 
mentary basins of the United States in addition to the Texas and Louisiana Coasts. 

Effective August 15th, the name of ELECTRO-TECH FRANCE will be changed to MANDREL 
CIE. 

For reasons of overall convenience, the headquarters and manufacturing facilities of Mandrel 
Cie will be moved from Paris to Annecy (Haute Savoie), France. However, a sales office for Electro- 
Technical Labs. will be maintained at Paris. 

The Paris sales office will be headed by GEoRGE KARAMBERIS. 

The new addresses are as follows: 

Mandrel Cie Electro-Technical Labs.—Sales Office 
59, Avenue des Romains 54, Avenue de Versailles 

Annecy (Haute Savoie) Paris 16 

France France 


Davip W. Strancway has completed requirements for the Ph.D. degree in physics at the Uni- 
versity of Toronto. He is employed as a research geophysicist for Bear Creek Mining Co., Denver, 
Colorado. 
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CALENDAR OF MEETINGS 
1961 


February 
19-23 AIME Annual Meeting, St. Louis, Missouri. 


March 
20-23 IRE National Convention, New York City. 
April 
9-11 Society of Exploration Geophysicists, 14th Annual Midwestern Meeting, Skirvin Hotel, 
Oklahoma City, Oklahoma. 
24-27 AAPG-SEPM Annual Meeting, Denver, Colorado. 


October 
8-11 Society of Petroleum Engineers of A.I.M.E., Annual Meeting, Dallas, Texas. 
18-21 AAPG Mid-Continent Regional Meeting, Amarillo, Texas. 
25-27 Gulf Coast Association of Geological Societies, San Antonio, Texas. 


November 
2-4 GSA Annual Meeting, Cincinnati, Ohio. 
5-9 Society of Exploration Geophysicists, 31st Annual International Meeting, Denver, Colo- 
rado. 
9-10 Pacific Section AAPG, Los Angeles. 
13-15 API, Chicago, Illinois 


1962 
March 
15-16 Fort Worth Geophysical Society, 15th Annual Midwestern Exploration Meeting, Texas 
Hotel, Fort Worth, Texas. , 
26-29 AAPG-SEPM Annual Meeting, San Francisco, California. 
26-29 IRE National Convention, New York City. 
September 
17-20 Society of Exploration Geophysicists, 32nd Annual International Meeting, Calgary, 
Alberta. 
October 
7-10 Society of Petroleum Engineers of A.I.M.E. Annual Meeting, Los Angeles, California. 
31-Nov. 3 Gulf Coast Ass’n Geological Societies, New Orleans, La. 


November 
12-14 GSA Annual Meeting, Houston, Texas. 


1963 


March 
25-28 AAPG-SEPM Annual Meeting, Houston, Texas. 


October 
6-9 Society of Petroleum Engineers of A.I.M.E. Annual Meeting, New Orleans, La. 
24-28 Society of Exploration Geophysicists, 33rd Annual International Meeting, New Orleans, 
La. 


November 
13-16 Gulf Coast Ass’n of Geological Societies, Hot Springs, Ark. 
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INDEX TO GEGPHYSICS VOLUME XXV 


Abstracts from Geophysical Prospecting . 
‘Acoustic Velocity in Porous Media, James E. Berry.” ” Review by K. F. Kenn edy 
Affleck, James, review of, ““An Outline of Geomorphology,”’ S. W. Wooldridge R. S. 
Morgan 
Affleck, James, review of, “Contemporary Geodesy ” edited by Charles Whitten and Ken- 
neth H. Drummond. 
Airborne —- Meter—Description and Preliminary Results, H. T. Lundberg and J. H. 
Ratcliffe 
Airborne Radioactivity ‘Surveys in Geologic Exploration, R.M.Moxham . . 
“Airborne Radiometric Survey of Portion of the Olary Province,” J. L. Harris. Review by 
S. H. Ward 
“Amplitudes of Seismic Body Waves from Underground ‘Nuclear Explosions,’ ” Carl 
Romney. Review by Joshua L. Soske 
“Analysis Charts for the Determination of True Resistivity from Electric Logs, ” Hubert 
Guyod and John A. Pranglin. Review by G. E. Archie : 
“Analysis Charts for the Determination of True Resistivity from Electric Logs,” Hubert 
Guyod and John A. Pranglin. Review by R. A. Wilson . 
Angona, F. A., Two-Dimensional Modeling and its Application to Seismic Problems . 


Announcements 


“Application of Cagniard’s Method to Dynamical Problems of Elasticity, On the,’’ M. 
Mitra. Review by C. Hewitt Dix. 

“Application of Ray Theory to the Problem of Long- Range Propagation of Explosive 
Sound in a Layered Liquid,” C. L. Pekeris, I. M. Longmen, and H. Lifson. Review by 
J. W. C. Sherwood 

Application of Seismic Methods to a Ground-Water Problem in Northeastern Ohio, R. E. 

Archie, G. E., review of, “Analysis Charts for the Determination of True Resistivity from 
Electric Logs,” Hubert Guyod and John A. Pranglin 

Archie, G. E., review of, “Electric Analogue for Resistivity Logging,” Hubert Guyod : 

Archie, G. E., review of, "Resistivity Determinations from Electric Logs,’’ Hubert Guyod . 


Basement Mapping with Aeromagnetic Data—Blind River Basin, Nelson C. Steenland and 
. J. Brod . 


Basham, W. L., review of, “Airborne Gravity Meter- Description a and Preliminary Results,” 
2. Lundberg and J. H. Ratcliffe : 

Basham, W. L., review of, “C ontinuous Gravity Measurements on a Surface Ship with the 
Graf Sea Gravimeter,” J. Lamar Worzel_ . 

Biake, F. G., review of, ‘Physical Geography,” Arthur N. Strahler 

Blau, a W., Memorial, E. E. Rosaire : ; 

Born, W. T., A Review of ee Instrumentation . 

J. S. , see Dunlap, H. F. 

Braun, T. H., ’ Electrical Analog Approach to Dipmeter Computations , 

Brod, R see Steenland, Nelson C. 

Broding, R. A., Collection and Processing of Digitized Acoustic c Log Data 

Bryant, H. L. , Production Well Logging Techniques ae 


Calendar of Meetings 


Case History of the Bronte (Ellenburger) and Fields, Coke 
Rollin E. Phipps. 
Chombart, L. G., Well Logs i in Carbonate Reservoirs . 
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INDEX 


Classics of Geopuysics, Nelson C. Steenland . 

Cobb, Howard, Operations Research—A Tool in Oil Exploration ; 
Collection and Processing of Digitized Acoustic Log Data, R. A. A. Broding and J. 'L. Poole . 
Collins, Francis, Plane Compressional Voigt Waves , 


Committees . 


Comprehensive System of Automatic Computation in Magnetic and Gravity Interpreta- 

tion, A. Roland Henderson . 
“Contemporary Geodesy,’’ Charles W hitten and Kenneth H. Drummond, editors. Review 

by James Affleck 

Continuous Gravity Measurements on a Surface Ship with the Graf Sea Gravimeter, ue 
Lamar Worzel. Review by W. L. Basham . 

Continuous Signal Seismograph, John M. C rawford, William E. N. Doty, and Milford R. 
Lee . 

“Contribution to a T heory of the Directional Effect of Groups of Se ismographs i in a Case of 
Pulse Vibrations, I and II,” V. S. Issayev. Review by J. E. White. a 


Contributors. 


“Converted Waves by a Narrow Canal on the Surface of a Semi-infinite Elastic Medium,” 
Ziro Suzuki and Akira Ishigaki. Review by J. W. C. Sherwood : 
Crawford, John M., Continuous Signal Seismograph : 


Danes, Z. Frankenberger, Successive Approximation Method for 
Anomalies, Ona, . 

Density Logging, J. j. Pickell and J. G. Heacock : 

“Dielectric Constant and Electrical Resistivity of Natural-State Cores, ” G. V. Keller and 
P. H. Licastro. Review by David A. Lowitz 

“Die Karte der mittleren Freiluftanomalien fur Gradabteilungen ‘6X10 von Westdeutsch- 
land” (Map of Mean Free Air Anomalies for a 6 ft 10 inch Coordinate Grid of West 
Germany), Karl Gerke and Hginz Watermann. Review by Sigmund Hammer . 

Discussions and Communications - 

Dix, C. Hewitt, review of, “Application of Cagniard’ s Method to Dy namical Problems of 
Flasticity,”’ M. Mitra 

Dix, C. Hewitt, review of, “On the Propagation of Transient Elastic W aves,’ ’ Shozaburo 
Nagumo. 

Doll, H. G., Suggestions for Better Electric Log Combinations and Improved Interpreta- 

Doty, William E. N., see Crawford, John M. 

Dumanoir, J. L., see Doll, H. G. . 

Dunlap, H. F., Marine Seep Detection—A New Reconnaissance Exploration Method 

Duties of a Logging Coordinator in an Operating Company, H. M. Johnson 


“Earth, Its Origin, History and Physical Constitution, The,” Sir Harold ee: Review 
by P. C. Wuenschel_ . 

Eckharat, E. A., Memorial, William Schriever . 

Editorial, Nelson C. Steenland 

Editoriai—Silver Anniversary of Gropuysics, Nelson C. Steenland . 

Editors of Geopuysics—Feature Section 

Effect ef Pore Fluids on the Elastic Properties of Sandstone, R. L. Mann and I. Fatt . 

Electrical Analog ‘proach to Dipmeter Computations, T. H. Braun and G. Y. W sais ; 

Electrica! Resistivit» Survey for Ground W ater. An, Tsvi Meidav é 

“Electric Analogs ot Resistivity Logging,” Hubert Guyod. Review by G. E. Archie . 

“Electric Analogue Re slativity Logging,’ Hubert Guyod. Review by R. A. Wilson . 

“Electric Prospect? ‘ by the Telluric Currents Method,” A. M. arcsicevual and N. M. 
Berdichevsk:) Review by J. L. Soske . 
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Electrochemical Mechanism of Sulfide Self-Potentials, The, Motoaki Sato and Harold M. 
Mooney 

Electromagnetic Response ofa Conducting Sphere toa Dipole Field, On the, James R. Wait 
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WILLIAM CROWE KELLOGG 
Geological Engineer 
Kellogg Exploration Company 
Geologists — Geophysicists 
Electrical — Magnetic — Gravity — Radioactivity 


Air-Ground Surveys Interpretation 
3301 NorRTH MARENGO ALTADENA, CALIFORNIA 


LOUISIANA 


HENRY SALVATORI 
Western Geophysical Company 
of America 
1116 Pacific Mutual Bldg. 


523 W. 6th Street 
LOS ANGELES 14, CALIFORNIA 


MISSISSIPPI 


EWIN D. GABY 
Delta Exploration Company 
Jackson Mississippi 
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OKLAHOMA 


MISSOURI 


LERoy SCHARON 
Mining and Engineering Geophysics 


567 Brookhaven Court 


Kirkwood 


Phone 
St. Louis, Missouri YOrktown 6-4245 
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NEW YORK 


E. V. MCCOLLUM 
Geophysicist 


E. V. McCollum & Co. 

Namco, International 
515 Thompson Building 
TULSA, OKLAHOMA 


Geo Seis, Inc. 


ROLAND F. BEERS 
ROLAND F. BEERS, INC. 


Petroleum and Minerals Exploration Consultants 
round and Airborne Surveys 
Data Reduction and Analysis 
Troy, New York 


P.O. Box 1019 
AShley 2-6478; 2-2351 


OKLAHOMA 


GLENN M. McGUCKIN 
Seismograph Consultant 
Data Interpretation Current Supervision 
SALES: McGUCKIN (Patented) SEISMOGRAPH 
SECTION PLOTTER (DIP MIGRATOR) 
25 Years Experience: 
Supervising, Contracting, Consulting 
1304 Ann Arbor 


Phone JEfferson 4-0853 
NORMAN, OKLA. 


THOMAS J. BEVAN 
Geophysicist 


914 American Airlines Bldg. Phone CHerry 2-7508 
TULSA 3, OKLAHOMA 


ARNOLD H. BLEYBERG 
Petroleum Exploration Consultant 
Domestic & Foreign Services 


MID-AMERICA BANK BLDG, 


OKLAHOMA CITY 2, OKLA. PH.: CENTRAL 2-0913 


CRAIG FERRIS 
Geophysicist 

E. V. McCollum & Co. 

Namco, International 
515 Thompson Bldg. 
TULSA 3, OKLA. 


GeoSeis, Inc. 


RICHARD A. POHLY 
Gravity Surveys and Re-interpretation 


POHLY EXPLORATION COMPANY 
1135 E. 38th St. Riverside 2-2009 


TULSA $, OKLAHOMA 


BEN F. RUMMERFIELD 
Geologist and Geophysicist 
CENTURY GEOPHYSICAL CORPORATION 


503 Jenkins Building 1105 7th Ave., W. 
Tulsa 3, Oklahoma Calgary, Alberta 


JOHN J. RUPNIK 
Petroleum Exploration Consultant 
GEOPHYSICAL & GEOLOGICAL COORDINATION 
J. J. RUPNIK AND COMPANY 


730 BEACON BLDG., Telephone LU 4-6355 
TULSA 3, OKLAHOMA 


Contractor and Research and 
Consultant Developments 
V. L. JONES 
TERRAMETRIC EXPLORATION COMPANY 
GEOPHYSICAL AND GEOLOGICAL SERVICES 
P.O. Box 3731 Phone 
TuLsa 23, OKLA. DI 3-0012 


HUGH M. THRALLS 
Consulting Geophysicist 
Gibson 7-3921 


Box 9577 
TULSA, OKLAHOMA 


TOM D. MAYES 
Mayes-Bevan Co. 


Gravity Meter Surveys 
and Interpretations 
345 Kennedy Building 


TULSA, OKLAHOMA 


Please mention GreopHysics when answering advertisers 


i 

18 
Mates 

: 

@ 
Bris, 

fi 
Be 


GEOPHYSICS, DECEMBER, 1960 


TEXAS 


TEXAS 


JOHN F. ANDERSON 
ANDERSON & COOKE 
Oil Exploration Consultants 


Geological Consulting 
Seismic Surveys G Interpretations 
665 San Jacinto Bldg. Houston 2, Texas 


PAUL FARREN 
Geophysical Consultant 
Specializing in Seismic Interpretation, 
Review, and Supervision 


1528 Bank of the Southwest 


HOUSTON 2, TEXas FA 3-1356 


WALTER D. BAIRD 
Consulting Geophysicist 


NEIL P. ANDERSON BUILDING 
PHONE EDISON 6-8400 FORT WORTH 2, TEXAS 


KEITH R. BEEMAN 
Electronic Consultant 
Seismograph and Allied Equipment 


P.O. Box 13058 
Houston, Texas 


L. F. FISCHER 
Exploration Consultant 


Geophysicist 
Geologist 


624 First City Bank Bldg. Houston 2, Texas 


JOHN L. BIBLE 
Gravity-Magnetic-Surveys-Interpretations 


Bible Geophysical Co., Inc. 
1045 Esperson Bldg. Houston 2, Texas 


J. F. FREEL 
RESEARCH EXPLORATIONS, INC. 


$134 Westheimer Road 
Houston, Texas 


HART BROWN 
Brown Geophysical Company 


Gravity-Meter-Surveys 
Interpretations 


5300 Brownway Rd. Houston 19, Texas 


JOHN A. GILLIN 


National Geophysical Company, Inc 
Namco International, Inc. 


2345 West Mockingbird Lane 
Dallas, Texas 


R. A. CRAIN 
Texas Seismograph Company 


1502 Eighth St. 
WICHITA FALLS, TEXAS 


T. I. HARKINS 
Independent Exploration Company 


1973 West Gray 
P.O. Box 13237 
Houston 19, Texas 


J. G. HARRELL 
Geophysicists 


2810 Primrose St. 
Fort Worth, Texas 


Telephone 
TE 8-0284 


R. H. DANA 
Dana Explorations, Inc. 


1301 W. T. Waggoner Bldg. 
Fort Worth, Texas 
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TEXAS 


TEXAS 


J. O. HOARD 


HOARD EXPLORATION 
COMPANY 


Esperson Building Houston, Texas 


C. T. MacALLISTER 
Geophysical Consultant 
Seismic Interpretations and Field Supervision 


6327 Vanderbilt, 
Houston 5, Texas 


Telephone: 
MA 3-4181 


W. B. HOGG 
Geophysical Consultant 


619 Fidelity Union Life Bldg. 
DALLAS 1, TEXAS 


JOHN S. IVY 
Niels Esperson Building 


HOUSTON, TEXAS 


HAYDON W. McDONNOLD 
Geologist & Geophysicist 


Keystone Exploration Company 


2813 Westheimer Road 
HOUSTON, TEXAS 


A. E. “SANDY” McKAY 
Geologist and Geophysicist 
Continental Geophysical Company 
Namco International, Inc. 

FORT WORTH, TEXAS 
1409 Continental Life Bldg., Phone EDison 2-9231 


J. C. KARCHER 
Geophysicist 


Adolphus Tower 
DALLAS, TEXAS 


GEORGE D. MITCHELL, JR. 
Geologist and Geophysicist 


911 Mercantile Securities Bldg. Dallas, Texas 


MARTIN C. KELSEY 
Rayflex Exploration Company 


6923 Snider Plaza Dallas 5, Texas 


H. KLAUS 
Geologist and Geophysicist 
Klaus Exploration Company 


Gravimetric and Magnetic Surveys 
and Interpretations 


P.O. Box. 1617 Lubbock, Texas 


P. E. NARVARTE 
Consulting Geophysicist 
Seismic Interpretations 

Current Supervision and Review 


Frost National Bank Building 
San Antonio, Texas 


PAUL H. LEDYARD 
Mid-Continent Geophysical Co. 


Contract Seismograph Crews 
Seismic-Reinterpretations 
2509 West Berry Fort Worth, Texas 


L. L. NETTLETON 


Gravity Meter Exploration Co. 
Interpretation of Gravity Surveys 
Interpretation of Aeromagnetic Surveys 


3621 W. Alabama 
HOUSTON 27, TEXAS 


W. W. (IKE) NEWTON 
Geophysicist 
823 Corrigan Tower 
DALLAS, TEXAS 
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TEXAS 


TEXAS 


J. O. PARR, JR. 
Consulting Geophysicist 
Methods & Instrumentation 

for 
Seismic & Radiometric Surveys 


202 Janis Rae San Antonio, Tex. 


RAYMOND L. SARGENT 
Magnetometer Surveys 
Interpretations 


M & M Bidg. 
HOUSTON 2, TEXAS 


C. W. PAYNE 
Exploration Consultant 
Geology—Geophystcs 


812 Continental Life Bldg. 
FORT WORTH 


SIDNEY SCHAFER AND COMPANY 
Geophysicists and Geologists 


EXPLORATION CONSULTANTS 
Domestic and Foreign 


Sidney Schafer 2200 Welch Avenue 
Jack C. Weyand Houston 19, Texas 


H. B. PEACOCK 
Consulting Geophysicist 


9746 Audubon Place 
DALLAS 20, TEXAS 


J. C. POLLARD 
GEOPHYSICAL ENGINEERING 
e 


SEISMIC . . . GEOGRAPH .. . GRAVITY 
MAGNETIC SURVEYS 


2500 Bolsover Rd., P.O. Box 6557, Houston 5, Tex. 


HUBERT L. SCHIFLETT 
States Exploration Company 
Seismic, Gravity and Magnetic Surveys 


5313 Richmond Road Highway 75 North 
HOUSTON, TEXAS SHERMAN, TEXAS 


H. B. SMYRL 
Portable Seismograph, Inc. 


706 Frost National Bank Bldg. 
San Antonio 5, Texas 


ROBERT H. RAY 
GEOPHYSICAL ENGINEERING 
SEISMIC . .. GEOGRAPH ... 
MAGNETIC SURVEYS 


GRAVITY 


2500 Bolsover Rd., P.O. Box 6557, Houston 5, Tex. 


SAM D. ROGERS 
Rogers Geophysical Company 


Rogers Explorations, Sociedad Anonima 


3616 West Alabama Houston 6, Texas 


NELSON C. STEENLAND 
Gravity Meter Exploration Company 
Geophysicist 


3621 W. Alabama Houston 27, Texas 


R. C. SWEET 
Geophysicist 


1111 Bering Dr. Houston, Texas 


C. WHITNEY SANDERS 
Consulting 
Petroleum Geologist-Geophysicist 


U. S. and European Exploration 
1132 Bank of the Southwest Building 
HOUSTON 2, TEXAS 


Seismic Reviews ‘ Field Supervision 
KIRBY J. WARREN 
Geophysical Consultant 


909 Continental Life Building 
Fort Worth, Texas ED 2-9073 
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TEXAS 


CANADA 


E. DARRELL WILLIAMS 


GEORGE W. SANDER 

Geophysicist 

onsultin eophysicist 

Specializing in Radiometric Exploration for Oil 8 paper 
3114 PRESCOTT ST., HOUSTON 25, TEXAS 


174 Douglas Ave. N., 
PHONE MA 3-3991 ouglas Ave 


Phone 
Oakville, Ontario 


Victor 4-6345 


JOHN H. WILSON 


Exploration Consultant 


W. F. STACKLER 
1201 Sinclair Bldg. Consulting Geophysicist 
FORT WORTH, TEXAS 


Phone CHery 4-7303 
1937 25th Avenue S.W. 
CHARLES C. ZIMMERMAN 


CALGARY, ALBERTA 
Geologist & Geophysicist 


Keystone Exploration Company 
2813 Westheimer Road For Information 
HOUSTON, TEXAS 


WYOMING On inserting your card 


Exploration Geology 


Seismic Reviews 
Evaluations 


in the 
Seismic Supervision 
JOHN F. PARTRIDGE, JR. 


Professi ir 
Consulting Geologist—Geophysicist rofessional D ectory 


P.O. Box 258 Room 211 O-S Building 
Phones 2-6485 and 2-3328 CASPER, WYOMING 


CANADA Write 


Business Manager 
R. E. DAVIS 
Farney Exploration Company, Ltd. 


Box 1536 


830-8th Avenue West 
CALGARY, Alberta, Canada 


Tulsa |, Oklahoma 


JOHN O. GALLOWAY 


Petroleum Consultant 


TRANSLATORS-SCIENTISTS WANTED 

805 Eighth Avenue South West Proven ability to translate technical ma- 
CALGARY, ALBERTA terial into fluent English essential. Attrac- 

tive full time or free-lance arrangement. 

All languages of interest, particularly Rus 

THEODORE KOULOMZINE sian and Japanese. Send résumé to: 
Geologist & Geophysicist 

Koulomzine & Brossard Ltd. A-T-S, Inc. 

P.O. Box 880, VAL D'OR Que. Canada. Drawer 271, East Orange, N.J. 

Rm 905, 80 Richmond St. W. Toronto, Ont. Can. 

Rm 1014, 132 St. James St. W. Montreal, Que. Can. 


AMherst 2-9018 
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BETTER 
SEISMOGRAPH 
SERVICE 


TEXAS SEISMOGRAPH CO.,INC. 


1502 EIGHTH ° WICHITA FALLS, TEXAS 


An easily read, clearly illustrated text on 


Modern Geophysical Techniques 
EXPLORATION GEOPHYSICS 


by J. J. Jakosky, Se. D. 


In more than 1200 pages and 
with 715 illustrations, the 1957 
revised printing of Exploration 
Geophysics covers the entire field 
of exploration by modern geo- 
physical methods. It is concisely 
and clearly written by an inter- 
nationally known geophysicist, in 
close collaboration with 39 other 
leading authorities, 


Thirteen fact-packed chapters 
fully cover all contemporary 
methods; plus permit, trespass 
and insurance problems. A basic 
textbook for every geologist, 
geophysicist, engineer and phys- 
icist concerned with exploration, 
well logging and _ production. 
Adopted by many leading uni- 
versities, 


Send your money order or check for $12.50 for a copy of Exploration Geophysics on 
5-day approval. If you are not fully satisfied, merely return the book in its original 
condition and your money will be promptly refunded. 

TRIJA PUBLISHING COMPANY, 2502 W. COAST HIGHWAY, NEWPORT BEACH, CALIFORNIA 


Please mention GeopHysics when ansv ng advertisers 
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Take a NEW READING ON GRAVITY as an 
Effective, Economical Exploration Method 


If you have not used Gravity, or if your use 
has been limited in recent years, it will pay 
you to reappraise the benefits of this fast, 
low cost exploration method. Instrumenta- 
tion, operational techniques, and interpreta- 
tion have been greatly improved enabling 
many operators to provide specific and de- 
tailed geologic information. Originally, the 
Gravity Method was thought of primarily 
as a means of determining salt dome config- 
urations. Now, stratigraphic and structural 
information can be obtained such as: the 
areal extent and configuration of sedimen- 
tary basins; the maximum thickness and 
variation in thickness of basin sediments; 
regional structures within the basins; local 
and regional fault systems, including horsts 
and grabens, within the section; and local 
and regional thrust faults. 


To obtain the best results, use the world’s 
most proved gravity instrumentation .. . 
There’s a WORDEN Meter 
for Every Land Application 


The MASTER and PROSPECTOR model 
WORDEN Gravity Meters provide the most 


demanding operating specifications ever 
achieved in portable meters. They retain the 
true portability and flexibility of the quartz 
element design plus the exclusive “Universal 
Compensation” feature which extends the 
temperature compensated range up to 6600 
mgls. (world coverage). The MASTER is 
distinguished from the PROSPECTOR by a 
low power Temperature Stabilizer, which 
maintains a nearly constant internal tem- 
perature in spite of extreme outside thermal 
shocks. Both meters are available in Stand- 
ard and Geodetic models. 


The PIONEER is ideally suited for gravity 
programs in areas of limited latitude and 
temperature variations and where the need 
for less frequent base ties exists. Also avail- 
able in Geodetic model. 


The EDUCATOR gravity meter is intended 
to meet the needs of educational or training 
programs where the required tolerances are 
much wider than for commercial programs. 


Specifications for the respective WORDEN Gravity 


Meters are shown at right . . . contact TI's Gravity 
Department for additional information. 


| 
: 


When requesting a quotation, copy the 
specification categories in this column 
and specify the values desired from 
those found under the respective 
meters. 


Pome 


Minimum Recommended Operating 
Range—Specify to nearest 100 mgls. 
plus latitude and elevation indication 


3000 mgls. 


3000 mgls. 


2400 mgls. 


Limited . . . Will be 
at least 1600 mgis. 


Minimum Total (Reset) 
—Specify to nearest 100 mgls.* 


5200 mgls. 


5200 mgls. 


4000 mgls. 


Limited . . . Will be 
at least 2000 mgls. 


Small Dial Range 


800 Dial Divisions x 
Small Dial Constant 
(2200 on Special Order) 


800 Dial Divisions x 
Small Dial Constant 
(2200 on Special Order) 


800 Dial Divisions x 
Small Dial Constant 


800 Dial Divisions x 
Small Dial Constant 


Small Dial Constant—Specify constant 
desired* 


Standard .08 to .11 
mgl./Dial Division (From 
0.05-1.00 also available) 


08 to 


1./Dial Division (From. 


0.05.1; 00 also availabie) 


Standard .08 to 
1./Dial Division (From 
0.05-1.00 also available) 


0.10-1.00 mgl./ 
Dial Division 


Small Dial Calibration Linearity over 
full range 


1 part in 1000 


1 part in 1000 


1 part in 1000 


1 part in 1000 


Smali Dial Reading Accuracy 


0.1 of Small Dial 
Division 


0.1 of Small Dial 
Division 


0.1 of Small Dial 
Division 


0.1 of Small Dial 
Division 


External Temperature Control 


Temp. Stabilizer re- 
quires 1.00 watts 
power for a 100° F. 
differential. Four 

chargeable “D” Size 
Batteries maintain 90° 
F. differential for 10 
hours and are mounted 
in removable pack on 
side of meter. Any 4.5 
to 12 volt DC source 
may be substituted. 


Side and top reading, directly coupled (no gears 
for backlash) with screw which contacts quartz 
measuring spring. 


Self counting, read from top and directly coupled 
(no gears for backlash) with screw which contacts 
quartz measuring spring. 


Small Dial Mechanism 


Large Dial Range — Specify range tat 


desired 800 Dial Divisions x Large Dial Constant 


Large Dial Constant Standard 6.5 to 12.0 mgls./D.D. Geodetics not available 


Large Dial Calibration over full range 
(8 turns) 


Curve furnished with values accurate to one part 


in 1000—measured over full range (8 turns). Geodetics not available 


Size, Inches: 
Height 14” 1014" 


Weight, Pounds: 
Net Ibs. 5% Ibs. 


Including carrying case 15% Ibs. 12 Ibs. 


*Price of meter is determined meaty by the 
Ranges and Dial Constant Specified. 


sei Bulletin No. GM-206 
° TEXAS INSTRUMENTS 


INCORPORATED 


GEOSCIENCES & INSTRUMENTATION DIVISION 
3609 BUFFALO SPEEDWAY © HOUSTON 6, TEXAS © CABLE: TEXINS 
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Exclusively By 
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Three 
books 
of interest 


from 


McGraw-Hill 


GEOPHYSICS, DECEMBER, 1960 


INTRODUCTION TO 
GEOPHYSICAL PROSPECTING 


New Second Edition 


By MILTON B. DOBRIN, Triad Oil Com- 


pany, Ltd. Calgary, Alberta, Canada. 446 
pages, $9.50 


A thorough revision of a highly successful text. 
It is designed to present the principles of cur- 
rent techniques of geophysical prospecting for 
oil and minerals to students and technical per- 
sonnel employed in the fields of petroleum and 
mineral exploration. The book covers all the 
major methods of geophysical prospecting. For 
each method it discusses fundamental physical 
principles, instruments, field techniques, reduc- 
tion of data, interpretation, and examples show- 
ing results of actual surveys. 


THE WORLD 
OF GEOLOGY 


Harvard 


in January, 1961. 


Edited by L. DONALD LEET, 
University, 
FLORENCE J. LEET. Ready 


PHOTOGEOLOGY 


By VICTOR C. MILLER, Mil- 
ler & Associates, Inc., Denver, 
Colorado. Ready in Spring, 1961. 
and 
A text for graduate and advanced 
undergraduate courses in photo- 
geology or photogeologic interpre- 
tation covering briefly the con- 


A lively compilation of readings 
on geology for the general reader. 
The origin of the earth, life, the 
atmosphere ; the ocean depths; the 
interior and exterior of the earth; 
the glacial movements; the forma- 
tion of land features; and the con- 
struction of the earth are all dis- 
cussed by specialists. 


cepts of aerial photography and 
photogrammetric instrumentation, 
and then, in greater detail, stere- 
oscopy and procedures. Of es- 
pecial interest is the final third of 
the text which contains stereo-pair 
aerial photographic illustrations 
accompanied by individual texts, 
exercises, and pertinent topo- 
graphic or geologic maps, sketches 
and cross-sections. 


Send for copies on approval 


McGraw-Hill Book Company, Inc. 
330 West 42nd Street 


Please mention GropHysics when answering advertisers 
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POINTING THE WAY 
FOR THE DRILL 


< 


Your complete geophysical service company 


Independent Exploration Company service is more 
complete than ever. In addition to its 

unparalleled experience and ability in traditional 
geophysical survey techniques IX offers you 

a wide range of skill and highly specialized equipment 
for gravity meter surveys, weight dropping 

surveys, offshore and inland marine surveys. 


To get more usable data for your dollar invested 
in exploration, call in... 


INDEPENDENT 
EXPLORATION CO. 


Geophysical Surveys 


1964 West Gray ¢ Houston, Texas 


Cable Address: Independex 
PARIS / LONDON / ALGIER / CASABLANCA / MADRID / TRIPOLI 
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ASKANIA Magnetic Instruments 
for precise determination of magnetic anomalies 


Magnetic Field Balance Gf 6. 
After Schmidt, the classic field 
instrument for relative measure- 
. ments of the horizontal or ver- 
tical intensity 


Torsion Magnetometer Gfz. Out- 
tanding g the gneti 
instruments for relative meas- 
urements of the vertical inten- 
‘sity 

Universal Torsion Magnetometer 
UTM. For relative measurements 
of all three components of the 
earth magnetic field 


Magnetograph. The new devel- 
* opment for recording also fast 
‘' occurring variations of one com- 
ponent of the earth magnetic 
field (either D-I-H-Z or T) 


ASKANIA-WERKE 
U. S. Branch Office & Service Dept. e 4913 Cordell Ave., Bethesda, Maryland 
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MAKE A STRIKE 


WITH 


With the proven, accurate, dependable interpretations of our 


gravity i» 


you are more apt to make a strike. Get the benefits of prompt, 
precise geophysical service through the application of our latest 
proven scientific methods. 


E. V. McCOLLUM CRAIG FERRIS E 'V. McCOLLUM & CO. 
515 Thompson Bldg. Ph. LUther 2-3149 Pf} 
TULSA, OKLAHOMA 


Foreign Affliate: 
NAMCO International, Inc. 
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GRAVITY and MAGNETIC SURVEYS and INTERPRETATIONS 

for direct correlation with subsurface conditions by 
preparation of Density Logs and Differential 
Residuals. 


PHONE CApitol 2-6266 WA 
1045 ESPERSON BUILDING 
HOUSTON 2, TEXAS GEOPHYSICAL CO., INC. 


JOHN L. BIBLE, President 


BERT F. DUESING, INC. 


“Selling Atlas Explosives” 
AND 


Manufacturers of Blasting Agents 


Magazines and Plant Magazines 


BIG LAKE, TEXAS HASKELL, TEXAS 


phone phone 
Big Lake 500 UNion 4-2456 


TECHNICAL TRANSLATORS WANTED 


Proven ability to translate technical material into fluent English essential. 
Attractive full time or free-lance arrangement. All languages of interest, 
particularly Russian. Send résumé to: A.T.S., Inc., Drawer 271, East 
Orange, N. J. 


Yearbook of the S.E.G. 


Watch For If... 


Please mention GroPHysics when answering advertisers 
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(Continued from page 9) 


can be scribed up to a diameter of 22”. 
Fits in brief case as a portable drafting 
tool. $3.95 ROL-RULER CO., Riegels- 
ville, Pa. 


ROTARY ROCK BIT 


A new hard formation rotary rock bit, 
type VHLC, is now available to the geo- 
physical industry in sizes 444”, 442” and 
434”, according to Varel Manufacturing 
Company, 9230 Denton Drive, Dallas 20, 
Texas. This new bit has cup shaped teeth 
on the heel row which provides faster 
penetration in hard formations. The 
VHI-C is designed for fast drilling in 
hard abrasive formations. The “cup- 
shaped” teeth have cutting action on the 
back and both sides. The center teeth are 
more rugged to withstand heavier weights 
which are uesd in formations with high 
compressive strength. New VH1-C bits are 
already in the field and reports indicate 
extremely good reception from field crews: 


SPECIAL TRAVEL SERVICE 


TRAVEL UNLIMITED, 103 East 4th, 
in Tulsa, Oklahoma, has a special service 
developed for the companies of the petro- 
leum industry that suddenly find them- 
selves with overseas contracts. 

Give this travel agency a call at LU 
7-8408, and tell them which country you 


31 


are interested in. They will line up all the 
information necessary for travel to that 
country. They will help you get passports, 
visas, residence requirements, hotel reser- 
vations, airline and steamship reserva- 
tions. They also advise you as to climate, 
clothes needed, family living conditions 
and supply details on interesting sights 
and side-trips you might want to plan. 
George Reinke, Manager, emphasizes, 
“One call to our organization can help 
anyone going overseas, and the best part, 
is that our services and assistance are 
given at no cost to you or your company.” 


PORTABLE DRILL 


The new Minuteman Portable Drill an- 
nounced by Mobile Drilling, Inc. is a 
compact, multi-purpose, rotary unit de- 
signed for subsurface exploration and pro- 
duction drilling in soils, rock and con- 
crete. Precision engineered from light- 
weight, high-tensile alloys, the Minute- 
man is believed to be the only truly porta- 
ble drill that features big rig rotary drive 
for greater drilling capacity and unusual 
versatility. The Minuteman, for example, 
handles continuous flight augers from 3” 
to 12” in diameter and drives 6” diameter 
augers to 30’ depths. When equipped for 
core drilling, it drives “EW” core barrels 
to 200’. When equipped for masonry 
drilling, the Minuteman cuts 1” to 8” 
diameter cores from steel reinforced con- 
crete in standard lengths required for 
quality control testing. As a production 


(Continued on page 36) 
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Belt or tape—for accurate 
magnetic recording of 
seismic data ask your 
supplier for 


REG. U.S. PAT. OFF. 


In geophysical work, as in dozens 
of other sciences, “SCOTCH” BRAND 


Magnetic Products provide the reli- 
able response needed for accurate 
data acquisition. 


Makers of some of the most sensi- 
tive tapes used in instrumentation, 
“SCOTCH” BRAND has pioneered in 


developing the basic magnetic prod- 
ucts which leading manufacturers of 


seismographic systems include in 
their own catalogues. Belt or tape— 
AM or FM recording—for top out- 
put ask your seismographic supplier 
for “SCOTCH” BRAND—and compare 


the results. 


“SCOTCH’’ BRAND MAGNETIC PRODUCTS 


for geophysical recording 


Minne ane company 
RESEARCH 15 THE KEY TO 


“SCOTCH” is a registered trademark of 3M Company, St. Paul 6, 
Minnesota. Export: 99 Park Avenue, New York, N.Y. In Canada: 
London, Ontario. 
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“WEEN ACCURAGY: 


World ‘Wide Gravity. 


Exploration crews the world over are using 
more World-Wide Gravity Meters than ever 
before —for good reasons! Weighing only 
eight pounds, this economical, built-to-take-it, 
portable meter operates anywhere in the 
world in all kinds of weather with maxi- 
mum accuracy. World-Wide’s easy-to-read 
meter is thoroughly temperature compen- 
sated, requires no thermostats, no barometric 
temperature corrections. Sealed in a vacuum, 
the World-Wide Gravity Meter gives depend- 
able, trouble-free service in the most difficult 
prospect areas. 
Exclusive Features Of 
The World-Wide Gravity Meter: 

@ Easy reading counter, ,, operator reads the 


meter without removing from the tripod. 
® Recessed level bubble, . , eliminates 


level bubble creep. 


@ World-wide range on all meters, 


regardless of latitude. 
e Approximately 100 milligal range on 


counter, ,, minimizes resetting instrument 
in rugged terrain. 


World-Wide Gravity Meters are available on purchase or rental/purchase plans. 
Each instrument carries a full-two-year warranty. Write or wire for complete details. 


woman INSTRUMENTS, INC. 


3802 South Shepherd, Houston, Texas ° Cable Address: GRAVIMETER HOUSTON 
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EUROPEAN ASSOCIATION 
| OF 


EXPLORATION GEOPHYSICISTS 


The E.A.E.G. was founded in December 1951 


The aims of the Association are to promote the science of exploration geophysics by 
establishing gontacts and encouraging co-operation and fellowship between geo- 
physicists in Europe and elsewhere and by disseminating knowledge of the science 
through the agency of regular meetings and the publication of technical papers. 


MEMBERSHIP 


Active Members pay an annual membership fee of Neth. fis. 15.—, increased by 
Neth. fis. 0.50 collecting charges. In U.S. currency this charge amounts to $4.10. 


Prospective Members. Anybody interested in geophysics can apply for membership 
by sending in an Application Form, duly filled out. Forms will gladly be supplied 
by the Secretary-Treasurer of the E.A.E.G. but may also be obtained through the 
kind assistance of the Business Manager of the S.E.G. 


GEOPHYSICAL PROSPECTING 


Official Journal of the European Association of 
Exploration Geophysicists 


This journal is issued quarterly and contains articles written in English, French or 
German. English, however, is predominant and each article is preceded by an abstract 
in that language. 

Active members receive the journal free of charge. 


In accordance with the Terms of Affiliation with the Society of Exploration 
Geophysicists, members of that Society may enter a current subscription to the 
journal for the normal membership fee. 


The Subscription Rate for non-members is Neth. fis. 22.—(U.S. $5.80) per annum. 
Single copies are available at Neth. fls. 6—(U.S. $1.60). These rates include pack- 
ing and postage and are payable in advance. A limited quantity of previous issues 
is still available at the same price. 


Advertising rates will be sent upon request. 


All communications to be directed to: 


THE SECRETARY-TREASURER E.A.E.G. 
30, C. VAN BYLANDTLAAN © THE HAGUE ¢ NETHERLANDS 
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(Continued from page 31) 


boring rig it also permits fast, economical 
erection of fences, guard rails, parking 
meters and similar installations. In addi- 
tion, the Minuteman handles standard 
soil sampling tools. 

The Minuteman weighs only 150 
pounds, yet features a 6 HP engine, 8 
speed automotive transmission,  dia- 
phragm-type carburetor for drilling at any 
angle, 44” stroke with power or hand feed 
both in and out of the hole, wheel 
mounted base for easier handling and two 
point anchoring for maximum stabiliza- 
tion in soils, rock or concrete. Although 
designed for one man portability, the 
Minuteman may be adapted for truck or 
trailer mounting as well. 

The Minuteman can be supplied as a 
Basic Package suitable for pre-bid boring, 
dewatering, borrow exploration, horizon- 
tal boring and similar operations. The 
Minuteman Basic Package includes all 
auger tools needed for drilling 3” di- 
ameter holes 25’ deep and is priced at just 
$1,485 F.O.B. Indianapolis, Indiana. A 
Minuteman Core Drilling Package and a 
Minuteman Masonry Drilling Package are 
also available in addition to a complete 
line of soils, rock masonry and_hori- 
zontal drilling tools and accessories. 

A new 4-page brochure which gives 
Minuteman specifications, operating fea- 
tures and tool options may be obtained 
by writing Mobile Drilling, Incorporated, 
960 N. Pennsylvania St., Indianapolis 4, 
Indiana. 


MILLIVOLTMETER 


Model 411A RF Millivolemeter 
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A new RF millivoltmeter which quickly 
and accurately measures RF voltages from 
1 millivolt to 10 volts at frequencies up to 
1,000 MC is now available from Hewlett- 
Packard Company. Readings are virtually 
free from the temperature effects normally 
associated with sensitive RF meters. 

The new millivoltmeter, Model 411A, 
includes two linear voltage scales in a 1:3 
ratio. Thes linear scales, an outstanding 
feature of the instrument, help eliminate 
scale reading errors and provides maxi- 
mum voltage resolution. The meter also 
includes a db scale for readings from 
—42 to +33 db. 

Model 411A is a compact, versatile unit 
which, according to the manufacturer, 
provides greater accuracy chan any instru- 
ments of its kind. Accuracy from 1 MC 
to 50 MC is +3% of full scale; 50 MC to 
150 MC, +6% of full scale, and 500 KC 
to 1,000 MC, +1 db. Because of the milli- 
voltmeter’s high temperature stability, 
accuracy changes are negligible from 10° 
to 40° C. 

Voltage range of the Model 411A is 10 
millivolts rms full scale to 10 volts rms 
full scale in seven ranges. The unit pro- 
vides full scale readings of 0.01, 0.03, 0.1, 
0.3, 1, 3 and 10 volts rms. 

Model 411A is equipped with a pen 
type probe tip, 500 KC to 50 MC. Extra 
screw-on probe tips include units for a 
high frequency measurement, for measur- 
ing on as well as at the termination of 
coax transmission lines, and a capacity 
divider increasing the instrument’s volt- 
age capability to 1,000 volts. 

The Model 411A RF millivoltmeter is 
priced at $450. Its current availability is 
9 wecks. 

Hewlett-Packard Company, 1501 Page 
Mill Road, Palo Alto, California. 
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DO YOU HAVE A 
NEW PRODUCT 
OR SERVICE? 


Send it to: 
GEOPHYSICS 
810 Shell Building 
Tulsa, Oklahoma 
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One galvo block... ay. 


From SIE comes the first triple-recording oscillograph requiring only one galvo block. 
Choose variable-density, variable-area, conventional wiggle or yarjable-density- 
igg mbined. Simple lens-box change converts from variable-density to variable 

area. m Timing lines may be generated internally fror either flasher tubes or standard 

timing motor. Consistency of paper drive is assured by a synchronous motor with built-in 

400 cps power supply. m= Answer all your oscillograph recording needs with one 

instrument, the TRO-6. 


SOUTHWESTERN INDUSTRIAL ELECTRONICS 
10201 Westheimer + P. 0. Box 22187 + Houston 27, Texas “Sunset 2-2000. 


For information on the Petty Stacking technique write to: 


Petty 


GEOPHYSICAL 

ENGINEERING CO. TRANSIT TOWER, SAN ANTONIO CApitol 6-1393 
Accurate subsurface surveys are 
accomplished by Petty’s patented 
horizontal stacking technique. 
Verification of each subsurface 
point is accomplished with 
multiple ray paths. Maximum 
record clarity is assured and noise 
is greatly reduced, since only the 
reflection point is common to all. ) U.S. Patent No. 2,732,906 
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KING WINCHES 


King Winches keep you moving 
through the most difficult terrain 
you get action where there's 
no traction with dependable pull- 
ing power. King power winches 
have pulling capacities of 8,000 
to 19,000 lbs. 


Model 550 Koenig Full Cab 
and Model 151J King Winch 
on CJ-5 ‘Jeep’ Universal. 


. . for ‘Jeep’ vehicles, International*, Chev- 
rolet, GMC, Ford, Land Rover and others. 


COMPLETE, READY-TO-INSTALL KING 
FRONT-MOUNT WINCH ASSEMBLIES 
FEATURE: 
@ winch side arms to reinforce truck 
frame 
e bronze-bushed, 4-way cable guide 
rollers 
@ cable drum guard 
e@ heavy-duty pipe bumper 
needle-bearing, universal -joint 
spline-shaft drive assembly 
e Timken bearings on worm 


*King Winches for International 
trucks are available through In- 
ternational-Harvester dealers. 


= 


[ 
LL-STEEL caBs FOR on Jeep 
UNIVERSALS 
FULL AND HALF CABS 


Koenig cabs and King 
Winches for ‘Jeep’ 
vehicles are available 

through all author- 
ized ‘Jeep’ vehicle 
dealers. Write for free 
descriptive literature. 


KOENIG ALL-STEEL CABS HAVE THESE SUPERIOR 
FEATURES: 

@ PROTECTION @ SAFETY 

@ COMFORT @ CONVENIENCE 
Roll-down windows, full opening . . . full panel- 
board head lining and masonite door lining... 
safety glass throughout . . . all-steel welded 
construction . . . door locks, 


IRON WORKS, Inc. 


P.O. BOX 7726R © HOUSTON 7, TEXAS 


Please mention GropHysics when answering advertisers 


: 
39 
j - 
‘ 
iy 
i 
| 
Le. 
* 
| 
\ 
| 


GEOPHYSICS, DECEMBER, 1960 


Why another 


seismic amplifier? 


SIE has a new amplifier for you! Five 
engineers and their technician as- 
sistants have worked for over a year 
to develop it. We have gone to the 
field until we are mud-caked, web- 
footed and gritty-toothed. We have 
checked and re-checked. Now it is 
ready for you. And the first thing you 
will probably say is “Why?” 


Why, with all the various seismic am- 
plifiers available, has SIE produced 
the new PT-100? Why another am- 
plifier? 


The answer is not original. Ped- 
dlers the world around have been 
using it since t = 0. “The PT-100 is 
not just another amplifier.” Lest you 
dismiss this answer as another ped- 
dler claim let us hasten to explain! 


From Calgary to Mexico City, from 
Tokyo to London we have heard dis- 
cussions about distortion. In fact, we 
heard it so often that we came back 
to our air-conditioned pastel towers 
and began to look critically at dis- 
tortion characteristics of the best 
amplifiers available. 


What we discovered bugged us. 
You were not only right, you weren’t 
even strong enough in some cases. 
We began to worder how many tricky 
reflection “picks” had gone wrong at 
critical spots because of the worst 
kind of distortion—that YOU 
COULDN’T DETECT ON THE REC- 
ORD. And that’s about when we be- 
gan work on the PT-100. 


For the first time ALL the advan- 


tages of transistor circuits have been 
utilized in one amplifier. Net only 
does the PT-100 provide low-power 
consumption, compact design and 
light weight, but it sports an entirely 
new AGC system. And it’s this AGC 
control that is the all-important dif- 
ference, the reason “why”. 


PT-100 signal-handling capabilities 
are such that it will not overload in 
ANY area. And now you can record 
refraction with full AGC performance! 


Forget about critical adjustments. 
Don’t worry any more about distor- 
tion destroying the accuracy of your 
“picks’—you'll actually see improved 
reflection standout. At the same time 
enjoy all the other advantages of 
transistorized design. No more re- 
placing parts like vacuum tubes to 
keep the channels balanced — tran- 
sistors don’t wear out. 


All these “features” and a long 
list more are true bonuses. But most 
important to you if you want read- 
able, reliable records is PERFORM- 
ANCE. Guess we'll just come right 
out and say it. The PT-100 is a sig- 
nificantly new amplifier. 


So finally you challenge, “Prove 
it!” Only three steps are needed. 


1. Look over the “capsule” informa- 
tion here. 


. Write for the PT-100 brochure. It’s 
complete; we want you to know 
ALL about this amplifier. 


3. Call SIE. 
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PT-100 Distortion-Free Transistorized Seismic Amplifier 


SIGNAL CONTROL 

With unique light-sensitive photocells as 
losser elements, the PT-100 eliminates 
the usual amplifier AGC problems. The 
critical “diode-balancing” common to 
other types of AGC circuits is never re- 
quired. There is no direct electrical cou- 
pling between AGC and signal paths to 
cause distortion and oscillation. Signal 
voltages do not cause distortion when 
applied across losser elements as with 
diode bridge circuits. Drift and tempera- 
ture instabilities are eliminated. There 
is a useful operating range of from one 
microvolt to one volt input — 1,000,000 to 
1! Total harmonic distortion under all 
normal field operating conditions does 
not exceed 0.5% including heavy reflec- 
tion standout. External 4 kc signals may 
be used for programmed-gain operation 
during amplitude studies. AGC perform- 
ance down to 5 cps allows use of the 
PT-100 in refraction work with AGC con- 
trol. 


FREQUENCY RESPONSE 
The broad band response extends from 


Sik 


3 to 500 cps. Low-cut filters are placed in 
the geophone line to eliminate ground 
roll distortion. Ten logarithmically-spaced 
cut-off frequencies are provided—from 16 
to 135 cps for the low-cut and 23 to 235 
cps for the high-cut filter positions. One 
or two sections (18 or 36 db/octave) may 
be used. 


NOISE 

Noise is less than 0.3 microvolt for “filter- 
out” bandpass and drops to 0.1 microvolt 
for normal filter settings—as good or bet- 
ter than most vacuum tube amplifiers. 


OPERATION 

Operator-oriented master controls provide 
separate record and playback adjust- 
ments for both “early” and “final” gain. 
Trimming adjustments for each channel 
are also provided. Complete built-in am- 
plifier test units make testing and control 
adjustments rapid and easy. 


The PT-100 will provide you a new dimen- 
sion in recording from refraction through 
normal reflection up to high frequency 
use. 


DIVISION 


ELECTRONICS 


10201 Westheimer + P. 0. Box 22187 + Houston 27, Texas + Si/nset 2-2000 


CABLE: SIECO HOUSTON 


MEXICO CANADA 


Industrial El 


TWX: HO-1185 


EUROPE 


SIE Mexico 
Edison 40-1 


5513 Third Street S. E., Calgary. 
Phone: 21.5787 Phone: Chestnut 3-0152, 3-0937, 3-0964, 3-6922 


(Conedo) Limited SIE Division of Dresser AG 
Alberta, Canada leboch z 


Zurich, Switzerland 
Telephone: 32 84 87/89 + Telex: 52 6 83 
Cables: Dresserzur Zurich 
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GEOPHYSICAL INSTRUMENTS DESIGNED AND 
MANUFACTURED BY FORTUNE ELECTRONICS 


Better 
Performance, 
Longer. 


Cost you less! 


For complete details or for an office 
demonstration, write or call: 


¥ Wests 
ig 
_vortune Electron ies, 


| 


Parentis - Lacq 
Hassi-Messaoud - Zarzaitine 
Tiguentourine - Hassi R'Mel 


COMPAGNIE GENERALE DE GEOPHYSIQUE 


50, RUE FABERT - PARIS 7° - PHONE : INValides 46-24 


| 
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MARINE GEOPHYSICAL SERVICES CORP. 
2418 TANGLEY 
HOUSTON 5, TEXAS 
Phone — JA 6-4428 Cable — MARGEO 
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Do all bits look alike 
to you? 


Here’s Why 

Bb, D. 
Bits Are Better 

For Your Budget... 


“Blue Demon” Bits .. . the 
original replaceable blade bits 
developed to cut cost-per-foot 

in shot hole drilling . . . continue to 
pace bit progress as drilling condi- 
tions change. 


The addition of cuttings deflectors 
to the face of Hawthorne’s patented 
replaceable blades . . . to substanti- 
ally reduce plugging between forma- 
tion layers ... gives the driller 
considerable versatility with drilling 
weights . . . which in turn results in 
faster hole... 


Plus the 4-Blade feature devel- 
oped for faster penetration in harder 


formations . . . with more cutting 
surface on bottom to keep the bit 
from walking... providing a 
straighter hole to load, and less strain 
on your drill... 

Plus the highest quality of ma- 
terials and craftsmanship in every 
Hawthorne “Blue Demon” product 


Plus ample stocks of popular shot 
hole sizes maintained in every active 
exploration area... 

. . . are reasons why more shot 
hole drillers use “Blue Demons” than 
all other drag-type bits combined. 
Try a set on your next hard-rock 
location. 


WRITE FOR ILLUSTRATED CATALOG 


Cable Address: HAWBIT 
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PORTABILI 


that pays for itself! 


24-Trace Ali-Transistorized Seismograph 


The EXPLORER* Seismograph is setting economy and performance 
records in SUMATRA, PERU, COLOMBIA, MEXICO, FRANCE, 
CANADA, MISSISSIPPI, and SOUTH LOUISIANA. 

This Seismograph FIRST can do the same for you by... 


e Cutting recording and support personnel in half. 

e Halving time at and between recording set-ups, thus doubling 
production and proportionately reducing cost per profile. 
Giving you a wide operational range of frequencies from high 
resolution to refraction . . . wide selection of AGC, fixed, and 
TVG gain controls. 

Offering unequalled portability for “back-packed” or heli- 
copter operations. 


The EXPLORER Seismograph pays for itself in reduced operating 
costs and increased production ... let a TI seismic engineer supply 


you with proof. 
* A trademark of Texas Instruments 


TEXAS INSTRUMENTS 
INCORPORATED 


A GEOSCIENCES & INSTRUMENTATION DIVISION 
~ 3609 BUFFALO SPEEDWAY * HOUSTON 6, TEXAS * CABLE: TEXINS 
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C. N. PAGE A. E. "SANDY" McKAY 


Central Office 
Continental Life Bldg. 
Fort Worth, Texas 


Division Offices 


WEST TEXAS ~ MID-CONTINENT ROCKY MTS, 
Midland, Texas Tuisa, Okla. Denver, 
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The Allison 201 Continuously 
Variable Filter 


! 


196] 
Yearbook 


HERE'S A NEW 
PASSIVE NETWORK ot the 
FILTER IN THE 
SUB-AUDIO RANGE 

The new Allison 201 Filter goes 

into the sub-audio range, yet retains S 4 = ‘ G > 
the desirable characteristics of 
Allison Filters in the audio range. 
The low noise, low distortion and 
excellent transient handling capa- 
bilities of the 201 make it excel- 
lent for heart studies, geophysical 
work, low frequency vibrations, 
servo-systems and similar sub- 
audio frequency spectrum studies. 


ALLISON 201 SPECIFICATIONS 

Impedance — 600 ohms 

© Passive network 

e 30 db per octave attenuation rate 


©@ Independent high cutoff and 
low cutoff sections 


Low insertion loss 
Smooth pass band 
Negligible ringing 


Frequency coverage— 1 to 
256 cps 


Shipping weight: 35 Ibs. 
Price: $695.00 FOB Factory 


Proved dependable 
in years of service 


Allison 
Laboratories, Inc. 
11301 OCEAN AVENUE 

LA HABRA, CALIFORNIA 
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carries out 


GEOPHYSICAL INVESTIGATIONS 


FROM THE AIR 


for 


ORE PROSPECTING 
OIL EXPLORATION 


AND 


CIVIL ENGINEERING 


ABEM 


ON THE GROUND 


GEOPHYSICAL INSTRUMENTS 


UNDER WATER 
e@ Airborne EM and Magnetic Surveys 


e Ground EM, Magnetic, Electric, Seismic, and 
Gravimetric Surveys 


e@ Underwater Seismic Surveys 


IN MINES AND @ Underground Magnetic Surveys . 
BOREHOLES 


For further details about these geophysical services and 
instruments write to 


THE ABEM COMPANY 


DANDERYDSGATAN 11, STOCKHOLM, SWEDEN, 
or contact your nearest ABEM agent 


U.S.A. ENGLAND CANADA | FRANCE SOUTH AFRICA AUSTRALIA 

Geophysical Craelius Company Moreau Woodard S.A. Craelius, 92 | Norse Industries J. J. Masur & Co. 

Instrument and Limited, || Clarges | & Co., Ltd., 1880 Av. des Champs- | (Pty.), Ltd., 80! Pty. Ltd., John 

Supply Co., 1616 Street, London, O'Connor Drive, Elysées, Paris 8e. | Transreef House, Street, P.O. Box 48, 
Broadway, Denver WI. Toronto 16, ; Marshall Street, South Melbourne, 
| 2, Colorado. Ontario. Johannesburg. A 
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A New Lock at VARIAN’S 
MAGNETOMETER FAMILY 


The “Proton Free Precession” types find diverse applications: 
an “Optical Pumping” instrument extends research capabilities. 


y made 
by Russell Varian's “Proton Free 
Precession” concept has made 
Proton magnetometers the world 
standard in less than a decade. 
Using this principle, total mag- 
netic intensities are quickly meas- 
ured, without critical orientation 
or calibration requirements. Re- 
cent application of the ‘Optical 
Pumping” concept has enabied 
Varian Associates to develop a 
bidi vapor 

with the highest sensitivity avail. 
able today. This new instrument 

total genetic intensi- 


ties on a continuous basis. 


AIRBORNE 


The V-4914 Airborne\ Magnetome- 
ter (proton) is a rugged, light- 
weight and compact\version of 
Varian’s Station Magnhetometer. 
Sensitivity and accuracy are of 
observatory quality. Coarse- and 


be 

tape. The virtually “indestructibie 
“bird” contains no mechanical or 
electronic parts. Varian’s \ V-4914 
is easily accommodated in $ingle- 
engined aircraft. Range: 22, 
73,400 gammas. \ 


For further information on magnetometer sales or 
leases in many countries of the free world — For 
loan of the Varian 16 mm color and sound movie 
on proton magnetometry — For copies of Varian's 
new Geophysics Technical Memorandum Series 
— call, wire or write instrument Division 


OBSERVATORY 


Varian's V-4931 Modular Station 
Magnetometer (proton) is a highly 
versatile and accurate monitor of 
the earth's magnetic field. Strip- 
chart analog records and digital 
readout are supplemented by 
Printed-tape output which may be 
engaged, as desired, to provide 
permanent digital records. The 
V-4934 Station Magnetometer (ru- 
bidium vapor) is particularly suit- 
able for measuring and recording 


OCEANOGRAPHIC 


An accessory “fish” towed up to 
five hundred feet astern of the 
hydrographic vesse! is the sens- 
ing adjunct to Varian’s V-4931 for 
oceanographic survey work. Suit- 
able for broad-scale undersea 
mapping or offshore geophysical 
exploration, it furnishes informa- 
tion of observatory quality. Ana- 
log and digital readout and rec- 
ords, accuracy and operating 
range are identical! to the observa- 
tory model. 


PORTABLE 


bn Varian M-49, a complete pro- 
less 
20 ‘Ibs., furnishes direct read- 
ings in gammas every six sec- 
onds. Sensitive to better than +10 
gammas, it requires no calibra- 
tion or levelling and is so versa- 
tile it can make equally accurate 
field records from surveys on 
land, in the air, or over water. 
Eight plug-in tuner ranges, equal- 
ly spaced from 19,000 to 101,000 

allow op ion through- 
out the world. 


RESEARCH 


Examples of current projects in 
Varian's and ing 
laboratories eee proton and 
rubidium miniaturized 

S: for mis- 
sile, sateliite, and deep-space ap- 
plications and instruments being 
developed for specific military 
applications, such as ASW. Va- 
rian’s practical contributions to 


rd y are ig to 
advance the state of the art. 


VARIAN associates 


MR & EPR SPECTROMETERS, MAGNETS, FLUXMETERS, GRAPHIC RECORDERS, MAGNETOMETERS. MICROWAVE TUBES, MICROWAVE SYSTEM COMPONENTS, HIGH VACUUM EQUIPMENT, LINEAR 


PALO ALTO43, CALIFORNIA 


Please mention GEopHysics when answering advertisers 


— 
i 
a\ 
| 
| 
j 
| 
| 
| 
| 
| 
| 
| 
wk 
4 
| 
| tions. 
| micropulsa 
| 
| 
| 
ke 
| 
| 
| 
| 
ne 
| | 
| 
| 
2 
la | ad 
| + 
‘ 
| 
| 


GEOPHYSICS, DECEMBER, 1960 51 


COMPETENCE ... based on twenty-eight years 
of professional service. 


EISMIC 


XPLORATIONS 


SEI OFFERS .. . a group of experienced personnel—operational and staff. 
Equipment and technical standards controlled by our own laboratory. 


Leaders in modern methods of acquisition of seismic data and its 
interpretation in terms of geologic structure. 


SEISMIC EXPLORATIONS, INCORPORATED 
HOUSTON, TEXAS 
Midland Shreveport Denver 
Foreign Affiliate : Compagnie Reynolds de Geophysique, 


9 Rue du Marquis de Coriolis, Paris, France 
B. P. 14 Hussein Dey, Alger, Algeria 
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(Continued from page 36) 


VERTICAL FILES 


Hit-or-miss . . . try to find the plan you 
want in these dusty, space-consuming piles. 
If you do find it, chances are it’s dog- 
eared, wrinkled, hard to use . . . wasting 
your time and running up your overhead. 


Efficient . . . now, transfer those plans 
to PLAN HOLD Vertical Files: 1200 
sheets to only 24”, on the wall out of the 
way. Each sheet hangs clean, flat, ready 
for instant identification and use, putting 
profits in your pocket. 

Plan Hold Corporation, 5204 Chakemco 
St., South Gate, California. 


, DECEMBER, 1960 


CUTTINGS DEFLECTOR 


A new cuttings deflector has been de- 
signed into the popular “Blue Demon” 
4-Blade Insert Bit manufactured by Herb 
J. Hawthorne, Inc., Houston, for seismic 
shot hole drilling. 

Designed to deflect drilled cuttings 
away from the center of the bit, and into 
the return circulation stream, these de- 
flectors are forged into the face of the 
blade. 

Equally effective with either water or 
air as the circulation medium, the deflec- 
tors tend to prevent plugging, which 
makes the bit particularly efficient in 
drilling through alternating layers of hard 
and soft formation. The deflectors give 
the driller considerable versatility with 
drilling weights, which in turn results in 
faster drilling. 

The “Blue Demon” Replaceable Insert 
Blades are designed to be completely 
worn out and discarded, and replaced 
with a new set in the same bit body, pro- 
viding the driller with a_ factory-con- 
trolled, on-the-drill bit service. 

Developed for faster penetration in 
harder formations, with more cutting sur- 
face on bottom to provide straighter hole 
without loss of gage, the “Blue Demon” 
4-Blade Insert Bits are available in all 
popular shot hole sizes, from all author- 
ized Hawthorne dealers throughout the 
world. 
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NEW FROM SIE 


Hi Fidelity 


Low Distortion 
Extended Range 


Future methods demand Hi Fidelity Full Spectrum FM Magnetic Recordings ! 


TRANSISTORIZED With the recent miniaturization of FM modules, more 
MINIATURE MODULES Geophysicists now specify the full spectrum SIE FM 
Recording System  “RecorData” — the only low power- 
drain all-purpose system for refraction, conventional and 

high-resolution seismic recording. 
In addition to the latest advancement in portability, the 

new “RecorData” (PMR-20) gives you: 


© System performance not dependent on tape quality. 
@ No modulation noise. 

© Wide frequency response (1 to 500 cps). 

© Negligible phase shift from 5 to 500 cps. 

© A complete recorder in two packages. 


DIVISION 


DRESSER ELECTRONICS 
10201 Westheimer + P. 0. Box 22187 + Houston 27, Texas * SUnset 2-2000 
CABLE: SIECO HOUSTON TWX: HO-1185 
MEXICO CANADA EUROPE 
SIE Mexico South ics (Canada) Limited SIE Division of 
Edison 40-1 5513 Third Street S. E., Calgary, Alberta, Canada Zuricl 


Phone: 21-5787 Phone; Chestnut 3.0153, 3-0937, 3.0964, 3.6922 Telephone: 8a a7; /89 52 6 83 
Cobles: Dresserzur Zurich 


VISIT OUR GEOPHYSICAL 
SERVICE CENTER 
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you can relieve reviewing 
geophysicists from this expensive, 
time-consuming 


waste motion 


COPY 


the ideal method of reproducing seismic records for review and 


exchange work Seismo Copy is an exclusive, new method of reproducing 
seismic records directly from field documents without advanced assembly 
from you. Individual documents are photographed in their entirety on 16 mil- 
limeter film. Records are individually exposed into a negative phase and indi- 
vidually exposed back into a positive phase. The results are a continuous record 
section, properly tied together without splicing. Seismo Copy is a high fidelity 
reproduction process that positively captures all the quality of your field 
documents. 

Through centralized processing The Industrial Photographers Company 
is able to offer this fast, efficient reproduction service at an exceptionally 
low cost. All seismic work is under the direction of an experienced geophysi- 
cist; and is, of course, held in the strictest confidence. 

For immediate service or for descriptive details of Industrial Photog- 
raphers’ complete reproduction and microfilm services write or call: 


& 


photographers company 


Reproduction Specialists M & M BUILDING, HOUSTON 2, TEXAS 
for the oil industry Telephone: CApito! 2-1973 
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WET. 
Sensitive 
SURVEYING 
ALTIMETERS 


TYPE 
FA-185 


i 
3 
4 

4 


DURABLE 
SELF-BALANCING 
CUSTOM CALIBRATED 


ACCURACY 0.1% @ wo 
FA-181 


STANDARD RANGES 

Minus 1000 to 3000 feet 
Minus 1000 to 6000 feet 
Minus 1000 to 15000 feet 
Special Ranges Available 


Write for additional information 


Y¥ WALLACE & TIERNAN INCORPORATED 


25 MAIN STREET, BELLEVILLE 9,NEW JERSEY 
IN CANADA! WALLACE & TIERNAN LTD.,WARDEN AVE., TORONTO 13, ONT. 
A-118.42 
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SUBSIDIARY OF WESTINGHOUSE AIR BRAKE COMPANY 
ENID, OKLAHOMA, U. S. A. 


TROJAN PLUGS | 
SURE GRIP TAMPS 
Z DYNA-POINTS 


DYNA PUNCHES 


Uf WS & KNIVES 
— 


EXPLORATION COMPANY 


900 W. Main St. - Phone 6254 - Nite 2-5150 -- Box 445 -- Palestine, Texas 
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GEOPHYSICAL SURVEYS 
IN 
MINING, HYDROLOGICAL AND 
ENGINEERING PROJECTS 


270 pages, Cloth 
N. fis. 16.50 ($4.50) 


N. fls. 12.—($3.25) to 2.A.E.G. and S.E.G. Members 


CONTENTS 


Twenty-one (17 new) separate papers, of which thirteen are 
concerned with mining problems from many countries and 


cover a wide variety of geological conditions. Three further 


papers discuss the application of electrical techniques to the 


very important problem of water investigations, the remainder 


being devoted to civil engineering problems, mainly dam sites. 


PUBLISHED BY 
THE EUROPEAN ASSOCIATION OF EXPLORATION GEOPHYSICISTS 
30, CAREL VAN BYLANDTLAAN, THE HAGUE, HOLLAND 


SUPPORTED BY 
THE SOCIETY OF EXPLORATION GEOPHYSICISTS 


Orders to be sent to the Secretary-Treasurer of the E.A.E.G. 
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Western Conducts More 
Marine Seismic Surveys 


Than All Other ) 
Contractors Combined! d 


Every year since 1954, Western Geophysical Company 
boats have done well over half of ALL contract offshore 
seismic explorations in the world. 

This record reflects the considered judgment of many of 
the most respected people in the petroleum industry on the 7/ ’ 

competence of Western men, techniques and instrumenta- / Y, ify» 

tion. When you consider your next geophysical survey, VOLT 

let us tell you more about the Western services that have GEOPHYSICAL COMPANY 
earned and held such confidence. A SUBSIDIARY OF LITTON INDUSTRIES 


Principal Office: 933 NORTH LA BREA AVENUE, LOS ANGELES 38, CALIFORNIA 


Western's land seismic parties 
and gravity crews 
are at work on four continents. 


AFFILIATE ANO REGIONAL OFFICES THROUGHOUT THE WORLD 


Please mention GeopHysics when answering advertisers 


58 
ats 
id 
— 
7 
ie 


For quick-look appraisal of taped data, 
use Kodak’s direct-writing paper 


It responds well 
to this preservation 
treatment 


The new direct-writing types of oscillographs are great allies 
of magnetic tape. They make the data not only visible but instantly 
visible. They give you the assurance you need that all is captured and 
it is safe to move on. 

If your direct-writing oscillograph happens to be loaded with 
Kodak Linagraph Direct Print Paper, you have another assurance. 
This paper responds to a simple preservation treatment with Kodak 
Linagraph Permanizing Developer. Then you are sure that your visible 
record will stay visible, eminently visible—even if you put the original 
record through the bright light of a diazo or photocopying machine. 

We shall give you specific instructions if you will drop us a note 
on how you handle your records and what you want to be tia to do 
with them. Write 


EASTMAN KODAK COMPANY 
=e Photorecording Methods Division 
Rochester 4, N. Y. 
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EXPLORATIONS, INC. 


Wholly «wned subsidiary, 
conducting foreign seis- 
mic explorations. Pres. 
J. F. Rollins 


Scientific Service Laboratories, Inc. 
Wholly owned subsidiary, designing, 
developing, producing scientific 
systems. Pres. W. B. Huckabay 


6923 Snider Plaza 
Phone EMerson 3-1531 
Dallas 5, Texas 

Cable: RAYFLEX 


The parent company, 
conducting domestic 
seismic explorations. 
Pres. M. C. Kelsey 


Wholly owned subsidiary, conduct- 
ing marine explorations with new 
Omni-Search system. 

Pres, E. F. McMullin 


MARINE, INC. 


NOW 


4 companies 
forged into a 
stronger 
chain of 
service... 


Please mention GropHysics when answering advertisers 


: 
bet 
| 
60 
ae 
. 
EXPLORATION CO. 
2 
By 
9 
: 
GROUP 
‘ 
| 
ae 


UNSCRAMBLING seismic data is the foremost reason for 
any seismic facility. In accomplishing this,recording accuracy, 
ease of programming, playback speed, and dependability 


are paramount. It is equally important that the system be 


versatile and easily adaptable to any of the multitude 


of data handling methods in popular use throughout the indus- 
try today. Data systems by Electro-Tech lead the field. 
Write or call for details today on our portable direct recording 


unit, and our pulse width and FM systems. 


ELECTRO-TECHNICAL LABS. 
MANDREL INDUSTRIES, INC. 


5134 Glenmont Houston 36, Texas 
Cable Address: ELECTROTEX 
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ENGINEERED SEISMIC SURVEYS 


6111 MAPLE AVENUE DALLAS TEXAS 
R. D. Arnett (.G. McBurney J. H. Pernell 
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DO YOU HAVE YOUR COPY? 


$5.50 (84.25 to SEG members) 
(Add 2% sales tax in Oklahoma; 


50¢ additional per copy on foreign orders) 


Order now from 


SOCIETY OF EXPLORATION GEOPHYSICISTS 
810 Shell Building Tulsa 1, Okla. 


Please mention GropHysics when answering advertisers 
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WHEN YOUR CONTRACT IS WITH ff GENERAL 


General’s geophysical services — exploration techniques, 
sinterpretation — are expressly coordinated to give you a 
on your oil prospect investment. THE PROOF? You enjoy 
of successful exploration. 


GEOPHYSICAL COMPANY 
WOUSTON CLUB BUILDING + HOUSTON, TEXAS 


In Canada: 10509 81st Avenue, Edmonton, Alberta, Canada 
General Geophysical Company de France (SARL), 4 Square Rapp 
General Geophysical Company de Venezuela, C. A., Apartado 1871 
General Geophysical Company (Bahamas) Ltd., Bogota, Colombia 


GENERAL GegpHYSicAl CO 


! WHEN YOUR CONTRACT IS WITH GENERAL, THE PERCENTAGE FOR SUCCES 


equipment design, 


better return 


a higher percentage 


SFUL EXPLORATION 


IN YOUR FAVOR 
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SEG PUBLICATIONS 


EARLY GEOPHYSICAL PAPERS 
(Compiled in 1947) 
Members of SEG 
Foreign postage, Raf, COPY 


Edited by L. L. Nettleton. 848 pages. 614 x 
914. Fully illustrated. Cloth. 


$12.00 
$ 8.00 
$ 0.50 


Papers presented at meetings of the Society of Pe- 
troleum Geophysicists during the years 1929 through 
1935 were published for the society in the journals 
of the American Association of Petroleum Geologists 
and the American Physical Society. Along with three 
special papers which had been mimeographed and 
distributed to SPG members in 1931 and 1932, the 
57 papers were assembled in 1947 with permission 
from the AAPG and APS and republished by offset 
process. The second printing of this historical work 
was prepared for publication in January, 1957. 
These papers contain much of the early history and 
technical fundamentals of the present geophysical 
industry, making this volume necessary to a com- 
plete geophysical library. 


INDEX OF WELLS SHOT FOR VELOCITY 
(1953) $5.00 
Members of SEG $2.50 


Edited by B. G. Swan, Continental Oil Co. 
68 pages. 634 x 914. Paper. 


This publication contains a list of over 2,000 wells 
previously published in Gropnysics. The original 
index appeared in GEopuHysics, v. 9, n. 4 (October 
1944), and supplements were published in v. 11, n. 
4; v. 14, n. 1; and v. 16, n. 1. This is a composite 
of the information contained in the four separate 
lists. Wells are listed by State and County, giving 
the name of the company, lease, location, survey 
depth, by whom shot, date, and by whom sponsored. 


INDEX OF WELLS SHOT FOR VELOCITY 
FOURTH SUPPLEMENT (1956) $1.00 
Members of SEG $0.50 


Edited by V. U. Gaither, Continental Oil 
Co. 23 pages. 634 x 914. Paper. 
This is a reprint of the supplement published in 
GEopHysics, v. 21, n. 1 (January 1956), listing wells 
shot during the period from July 1952 to October 
1955. 


FIFTH SUPPLEMENT (1957) $1.00 
Members of SEG $0.50 


This is a reprint of the supplement published in 
GEOPHYSICS, v. 22, n. 1 (January 1957). 


GEOPHYSICAL CASE HISTORIES, 
VOLUME 1 (1948) $10.00 
Members of SEG, AAPG & AIME $ 6.00 
Foreign postage, per copy $ 0.50 


Edited by L. L. Nettleton. 680 pages. Fully 
illustrated. 7 x 10. Cloth. 


This volume, first of a series to be published at five 
to seven year intervals, is a collection of 60 papers 
by 61 authors on geophysical observations made 
under a wide variety of field circumstances. The 
purpose of the series is to provide material by which 
geophysical surveys can be judged from later devel 
opment and thus aid in the interpretation and eval- 
uation of other geophysical work. Contents: 3 general 
and historical papers; 21 salt dome case histories 
(Texas, Louisiana and Mississippi); 17 mid-continent 
case histories (Arkansas, Illinois, Oklahoma and 
Texas); 4 Rocky Mountain case histories: 1] Califor- 
nia case histories; 4 foreign case histories. Features 
a complete title index cross-referenced to a classified 
index of all maps and figures. 


GEOPHYSICAL CASE HISTORIES, 
VOLUME II (1956) $10.00 
Members of SEG, AAPG & AIME $ 6.00 
Foreign postage, per copy $ 0.50 


Edited by Paul L. Lyons, Sinclair Oil & Gas 
Co. 690 pages. Fully ilJustrated. 7 x 10. 
Cloth. 


The second volume of the series contains 53 papers 
by 75 authors. In this volume the case histories are 
grouped by type of trap rather than by geologic 
province. It contains 17 foreign papers—l3 more 
than Volume I, 3 more than half the material is 
original, being published for the first time. The 
second volume also contains a complete index cross- 
referenced to a Classified index of all maps and 
figures. Contents: #6 general and historical, 6 salt 
dome, 6 reefs, 23‘ anticlines, 3 stratigraphic traps, 
3 mining and 5 new uses. 


CUMULATIVE IMDEX, 1931-1953 (1955) $5.00 
Members of SEG $3.00 
Foreign postage, per copy : $0.50 


Edited by Kenneth L. Cook, Department of 
Geophysics, University of Utah. 322 pages. 
7 x 10. Cloth. 


All publications of the Society (GEOPHYSICS, EARLY 
GEOPHYSICAL PAPERS and GEOPHYSICAL. CASE HISTORIES, 
VOLUME 1) during the period are indexed by subject, 
title and authors. Geophysical patents abstracted in 
ceopHysics from 1939 through 1953 are indexed 
separately by subject, patent number and inventor. 
Coverage of all literature reviews substantially ex- 
tends the usefulness of this index, since most of the 
significant literature of exploration geophysics since 
1936 not published by the Society has been reviewed 
in GEopHYsics. Whether or not the reader has a 
complete file of seG publications, he will find this 
index most useful. 


Society of Exploration Geophysicists 
Box 1536, Tulsa 1, Oklahoma 


Please mention GeopHysics when answering advertisers 
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SEG PUBLICATIONS 


SIXTH SUPPLEMENT $1.00 
Members of SEG $0.50 


This is a reprint of the supplement published in the 
1957 GEOPHYSICS YEARBOOK (December 1957). 


SEVENTH SUPPLEMENT $1.00 
5, Members of SEG $0.50 


This is a reprint of the supplement published in 
Geopuysics, v. 23, n. 5 (December, 1958). 


LESSONS IN SEISMIC COMPUTING 
(1959) $5.50 
Members of SEG $4.25 
Foreign postage, per copy $0.50 


Written by M. M. Slotnick, edited by R. A. 
Geyer. 350 pages, 7 x 934. Fully illustrated. 
Cloth. 


An elementary text and problem book. Its 350 pages 
contain 44 lessons in seismology arranged for selec- 
tion or combination to cover the normal 36-week 
course, or for condensation into an 18-week course. 
Written by Dr. Slotnick for the Geophysics Depart- 
ment of Humble Oil & Refining Co., this series of 
lessons presents the basic physical laws relevant to 
seismic interpretation. The lessons begin without 
assuming more than secondary school mathematics. 
An elementary knowledge of the Calculus is desirable, 
though not required, for the last half of the book. In 
its scope, detail and clarity of style, this work pro- 
vides an authoritative reference for individual study, 
as well as a textbook for formal courses. 


YEARBOOK of the Society of Explora- 

tion Geophysicists (Available Jan. 
1961) $3.00 
Members of SEG $2.50 
Foreign postage $0.20 


Published annually in January 


World-wide professional directory of exploration geo- 
physicists. Alphabetical listing of SEG membership 
with the grade of membership indicated. Geographic 
listing by companies. Annual report of the SEG. New 
committees for the current year. Directory of SEG 
local sections. Listing of affiliated societies. Pictures 
of SEG officers. Report of the annual international 
meeting of SEG and abstracts of the papers. Consti- 
tution and by-laws, as amended. Calendar of future 
meetings. Code of ethics. 


GEOPHYSICS (Bimonthly) 
Subscription, per year $10.00 
Foreign postage, per year $ .50 
Members of SEG No charge 


The official journal of the SEG is published in Feb- 
ruary, April, June, August, October and December 


each year beginning 1960. (The first issue was pub- 
lished in January, 1936 as Volume I, No. 1. Publica- 
tion was quarterly through 1957. There were five 
issues in both 1958 and 1959.) Volumes run on a 
calendar year basis. An average of eight technical 
og appear in each issue, all contributed by mem- 

Ts and others engaged in exploration geophysics. 


GEOPHYSICS—Clothbound Volumes 


SEG Foreign 
Members Others Postage 


$ 6.00 
6.00 


Reprints 


s 
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Vol 
Original Printing— 
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Vol. 25 (1960) 


Available issues, each 
Members of SEG 
Foreign postage 

Issue 


y 
October 
January 
April 
July 
October 
December 


Yea 


December 
February 
April 

July 
October 
December 
February 
April 

June 
August 
October 
December 
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a Vol. 2 (1937) 10.00 
Vol. 3 (1938) 10.00 
a Vol. 4 (1939) 10.00 
Vol. 5 (1940) 10.00 
Ee Vol. 6 (1941) 10.00 
ee Vol. 7 (1942) 12.00 
Vol. 8 (1943) 12.00 
Vol. 9 (1944) 12.00 
Vol. 10 12.00 
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ENGINEER'S FIELO NOTES 
Praject No 10-863-S4 


Report from Rogers 


The true value of the final, the complete geophysical report stems from 
many mechanics, many techniques, many interpretations that follow the 
first survey. The final Rogers Geophysical Report is an expert and thus 
complete analysis of the known and unknown of your potential oil province. 
That is because the expertness of Rogers extends to all the vital links 
between the Report and the first planning. Completeness and expertness 
are Rogers’ bywords in geophysical prospecting. 


OGERS Geophysical Companies 
3676 WEST ALABAMA « HOUSTON, TEXAS 
Edificio Republica « Caracas, Venezuela 
Mogadiscio « Somalia 
34 Ave. des Champs Elysees « Paris, France 
FOREIGN OFFICES 1-3 Arlington St., St. James's « London 1, England 
Hote! Castellana Hilton « Madrid, Spain 


Tripoli, Libya « Algeria 
ROGER CREW Go EVERVYVWHERE 
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Hot ship plays it cool on desert exp/ora- 
tions. The 305 hp Hiller 12 E is as untroubled 
by wind-blown sand as a well-built truck; as 
powertul as the biggest bul/dozer. When you 
put a 12 E to work in the desert, it often car- 
ries the same seismic crew and gear on one 
flight that other copters could transport only 
in two or more. 

The Hiller has proved in commercial service 
that it operates effectively and economically 
under full load at high altitude and at any 
latitude. It can even carry the fue/ for six 


hours’ cruise with profitable payload re- 
maining. It's this margin of power that could 
be vital to your operation. 

All this and more reasons, too. Your Hiller 
72 E operator will demonstrate them to you— 
show you why more and more charter con- 
tracts specify Hiller 12 E. Write for booklet, 
“New Workhorse for Petroleum.” It gives full 
information on the 12 E and on Hiller charter 
operators’ trip or contract services every- 
where. Address Commercial Division 


HILLER® 


AIRCRAFT 
CORPORATION 


PALO ALTO, CALIFORNIA 
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COMPLETELY ENCLOSED, SELF LEVELING 


provides accurate 
gravity readings 
even in rough 
water and on 

soft bottoms 


The exclusive seismic compensation on the LaCoste & Romberg Underwater 
Gravity Meter is added assurance of the ultimate in accuracy in offshore 
gravity surveys. To compensate for seismic motion of the ocean bottom, a 
servo-operated elevator* automatically raises and lowers the meter. With 
this seismic compensation, the accuracy of the LaCoste & Romberg Under- 
water Meter is within 0.02 mgl under good conditions and is better than 
0.1 mgl under such adverse conditions as rough weather and soft bottoms. 
Drift is within 1 mgl per month. 


You will find, too, that LaCoste & Romberg Gravity Meters—either 
underwater or land—are exceptionally rugged and dependable. Experi- 
ence since 1946 with underwater meters has enabled LaCoste & Romberg 
to build gravity meters that meet the practical requirements of commercial 
offshore exploration work. *PATENT NOS. 2,589,709 AND 2,589,710 


For detailed information about LaCoste & Romberg Gravity Meters, Write— 


LaCoste & Romberg e 6606 North Lamar @ AUSTIN, TEXAS 
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Geograph means ADDED service, ADDED advantages from the Robert 
H. Ray Companies. It is the newest, proven method for gathering sub- 
surface data. 

The Geograph weight replaces dynamite as an energy source. This 
technique has proven superior in reflection efficiency in all types of terrain 
—from frozen muskeg to desert sands. 

Geograph, known as the “Thumper”, improves data quality, especially 
in the poor record areas. In other regions, results have been comparable 
to, or better, than those obtained with other seismic methods. 

Just as research and development, modern instrumentation and elec- 
tronic processes have made weight dropping a proven, practical service, 
they have also provided an economical method of seismic coverage. Costs 
per mile have been substantially reduced. 

A proven method, improved record quality, economy and versatility 
are just a few of the advantages ADDED to your geophysical program 
when you specify Geograph. 

Whether your prospect requires Gravity, Magnetic or Seismic, the 
Robert H. Ray Companies provide these services world wide. 


For Geograph in Europe and the French Zone contact GEOGRAFRANCE, 48 Bivd. de Latour-Maubourg, Paris 7e. 


RAY GEOPHYSICAL € 
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EXPLORATION BY GEOPHYSICS 


HANNOVER - HAARSTRASSE 5 - PHONE: 86661 - TELEX: 922847 - CABLE: PRAKLA 
GERMANY 
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FAULTS ® SALT OVERHANGS 
UNCONFORMITIES ® CROSS BEDDING 
REGIONAL AND STRUCTURAL DIPS 


BARS AND PINCHOUTS 


the Schlumberger 
Continuous Dipmeter 


marks every 


geological feature! 


Closely spaced dip calculations have 
clearly marked the fault on the graph 
above. Other structural or depositional 
features are likewise recognizable by their 
characteristic dip patterns. 


FO TME te DUSTRY® 
a * The Continuous Dipmeter gives you a 


SCHLUVURA BERGER wealth of geological information from a 
single bore hole. Call your nearest 
Schlumberger engineer for a discussion 
of the outstanding results obtained with 
the Continuous Dipmeter. 


| 
4 


GEOPHYSICS, DECEMBER, 1960 


AVAILABLE PUBLICATIONS OF 


The American Association 
of Petroleum Geologists 


Comprehensive Index of the Publications of the A.A.P.G., 1917-1945. Compiled by Daisy 

Winifred Heath. 603 pp. 6.75 x 9.50 inches. Cloth. To members, $3.00 

Comprehensive Index of Publications of the A.A.P.G., 1946-1955. Compiled by Daisy 

Winifred Heath and June McFarland. 302 pp., 6.75 x $.50 inches, Cloth. To members, 

Both Indexes at once, covering all Association publications, 1917-1955 at special price. 

To members, $5.00 

Appalachian Basin Ordovician Symposium. From August, 1948, Bulletin. 264 pp illus. 

6 x 9 inches. Cloth. To members, $1 

Possible Future Oil Provinces of the United States and Canada. 4th printing. From 

August, 1941, Bulletin. 154 pp., 83 figs. 6 x 9 inches. Paper. To members, 

Problems of Petroleum Geology. 2d printing. Originally published, 1934. 43 papers. 1,073 

pp. 200 illus. 5.75 x 8.5 inches. Cloth. To members, $4.00 

Possible Future Petroleum Provinces of North America. From February, 1951, Bulletin. 

360 pp., 153 figs. 6 x 9 inches. Cloth. To members, $2.50 

Geological Cross Section of Paleozoic Rocks: Central Mississippi to Northern Michigan. 

Prepared under auspices of Geologic Names and Correlations Committee. 5 cross sections, 

vertical scale 500 feet to the inch. 29 pp. of explanatory text, index. 8 x 10 inches, Press- 

board, sections folded in pocket. To members, 

Miocene Stratigraphy of California (1938). By Robert M. Kleinpell. 45 -» 14 figs., 22 

pls., 18 tables. Offset reprinted. 5.5 x 8.5 inches. Cloth. To members, $4.50 

Petroleum Geology of Southern Oklahoma. 24 articles. 402 pp., 110 figs., 6 pls., 27 tables. 

6.75 x 9.5 inches. Cloth. To members, $5.00 

Structure of Typical American Oil Fields. Vol. I (1929) 4th printing. 510 pp., 190 illus. 

5.5 x 8.5 inches. Cloth. To members $4.00 

Structure of Typical American Oil Fields. Vol. II (1929). 4th printing. 750 pp., 235 illus., 

5.5 x 8.5 inches. Cloth. To members $5.00 

Structure of Typical American Oil Fields. Vol. III. McCoy Memorial Volume. 24 papers. 

516 pp., 219 illus. Cloth. To members, $3 

Geology of California (1933). By R. D. Reed. 355 pp., 58 figs. , 26 tables. Structural Evolu- 

tion of Southern California (1936). By R. D. Reed and J. S. Hollister. 157 pp., 57 figs., 14 

photographs, 9-color tectonic map. Both offset reprinted. 2d printing. 5.5 x 8.5 inches. 

Clothbound together. To members, $6.00 

Western Canada ge ge Basin. Rutherford Memorial Volume. Edited by Leslie M. 

Clark (1954). 30 papers. 521 2d printing. Stratigraphy of Plains of Southern Alberta. 

Dowling Memorial Volume {1931). A symposium. 14 papers. 166 pp. 2d printing, by offset. 

6.75 x 9.5 inches. Clothbound together. To members, $6.00 

Jurassic and Carboniferous of Western Canada. Allan Memorial Volume. E 

Goodman, 24 papers. 514 pp. Cloth. To members, $6.00 

Habitat of Oil. Edited by Lewis G. Weeks. 1,364 PP.» cloth bound. 56 papers on oil occur- 

rence—20 on North America; 5 on South America; 7 on Europe; 5 on Middle East; 4 on 

Far East; 9 topical papers on oil occurrence and migration, hydrocarbons, and basin 

development; and an analysis by the Editor. To members, $9.00 

Lower Tertiary Biostratigraphy on the California Coast Ranges. By V. Standish Mallory. 

Companion volume to Miocene Stratigraphy of California, by Robert M. Kleinpell. 297 

pp. of text; 7 line amas: 42 plates of Foraminifera; 18 tables; index. 6 x 9 inches. 

Cloth. To members, $11.00 

tery Geology of Southern Oklahoma, Voi. II. 17 articles, 350 pp., 141 figs., 3 pls. 
6.75 x 9.5 inches. Cloth. To members, $5.00 

Stratigraphic Type Oil Fields (1941). 37 papers. 902 pp., 299 figs., 3 pls. Offset reprinted. 

5.5 x 8.5 inches. Cloth. To members, $4.50 

Stratigraphic Cross Section of Paleozoic Rocks, West Texas to Northern Montana. Pre- 

pared under the auspices of the Committee on Stratigraphic Correlations. 6 cross sections; 

vertical scale 400 feet to the inch. Explanatory text, index. 8 x 10 inches. Pressboard, 

sections folded in pocket. To members, 

Bulletin of The American Association of Petroleum Geologists. Official monthly publication. 
Each number, approximately 150 pages of articles, maps, discussions, reviews. Annual 
subscription, $18.% (outside United States, $19.00). Descriptive price list of back num- 
bers on request. 

Available back issues ef the Bulletin through Vol. 39 (1955) are distributed at $2.00 per issue, 
plus postage, by: Walter J. Johnson, Inc., 111 Fifth Avenue, New York 3, N.Y. Send orders 
and remittances to that address. 


(Prices, postpaid. Write for discount to colleges and public libraries.) 


THE AMERICAN ASSOCIATION OF PETROLEUM GEOLOGISTS 
BOX 979, TULSA |, OKLAHOMA, U.S.A. 


Please mention GreopHysics when answering advertisers 
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CENTURY serves the world 


Century Geophysical! Corporation, with its 44 geophysical units presently spread across the globe, is backed and 
strengthened by the research and development laboratories and the manufacturing facilities of Century Elec- 
tronics & Instruments, Inc. These two affiliates are working partners in providing many industries with technical 
services and electronic instrumentation and equipment. 


Seismic and gravity meter crew operations, review geophysical interpretations, velocity surveys, data center 
processing, uranium exploration logging and radiometric analysis are some of the general areas where Century 
has gained a large backlog of experience and wide acceptance. : 


Century’s facilities are available to assist you in gathering, pro- 
cessing and evaluating scientific data anywhere in the world. 


Century Geophysical Corporation Century Electronics & Instruments, Inc. 


TWX TU-1407 — TULSA, OKLAHOMA 


Please mention GropHysics when answering advertisers 


75 
che 
Ti 
O\ | 
4 
| 
| 


GEOPHYSICS, DECEMBER, 1960 


INPUT FOR SIMULATORS DISPLAY DEVICES 


SURVEYING GEOPHYSICS 


VEHICLE TRACKING AVIONICS 


INSTRUMENTATION 


e 
e 
e 
e 
e 
e 
e 
e 
e 
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Aero Systems Engineering has 
unique capabilities in these fields 


Aero Service Corporation has made many significant contributions in 
fields related to photogrammetry, geophysics, special simulators, train- 
ing devices, avionic testing, and data handling and analysis. The 
authoritative position we enjoy in many of these sciences is based on 
skills developed through 41 years of experience. Perhaps we may help 
you, too. For our new facilities report, simply use the coupon below. 


er® SYSTEMS ENGINEERING 


a new division of Gentlemen — Please send a facilities report to: 
AERO SERVICE Nome 
CCRPORATION Company 


210 East Courtiand St. Address 
Philadelphia 20, Pennsylvania 


Please mention GEeopHysics when answering advertisers 
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. . . you'll get an accurate picture of 
oil-producing possibilities. Tidelands’ 
experienced crews and modern equip- 
ment assure positive results and high 
production. 


A Complete Geophysical Service 


&THEIMER ROAD. HOUSTON 27, TEXAS 


In the Hands 
of Experience... 


A 


useful 


ing problem of this type. 


Data recorded as curves can often be made more 


if it is converted to digital form. This is espec- 


ially true where the data is to be used with electronic 
computers. 

Seismograph Service Corporation now provides a 
complete data conversion service, at a cost far below 
manual conversion methods. This service includes 


who can work with you on amy data handl- 


Data to be converted must be in rectangular coordinate 
form. The input curves are tracked by a curve follower 
to convert the information into electrical impulses. The 
final output, coded punched tape, can be fed directly into 
an electronic computer or translated to punched cards, or 
listed in tabular form. A 10-inch original curve width will 
convert to as many as 1000 data samples on the punched 


tape. The long dimension of the curve can be conformed 


to any one of the many scales available. 


Some examples of how this data conversion is used 
in geophysics: Digitized data from velocity logs are 
the basic raw materials used for constructing syn- 
thetic seismograms, which inciude multiples. Seismic 
record traces can be digitized for frequency analysis, 
auto-correlation, etc. 

This service is available NOW. Call or write 
SSC today and find out how data conversion can 
be effectively applied to your needs. 
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... trace for trace 

the HS-J consistently 
matches the performance 
of larger and heavier 
single coil geophones. 


What the seismic records cannot show is the economy, efficiency 
and improved field techniques which the HS-J makes possible. 

Why do the job the hard way when this new small geophone opens 
a whole new area of faster, lower cost field operations? 


ECONOMY Remarkably low initial cost, 
lighter less expensive cabling, savings in 
crew man-hours—these are some of the rea- 
sons why more and more of the industry is 
turning to the HS-J. 


PORTABILITY The HS-J’s small size and 
light weight make the use of small hand 
reels and mechanical handling equipment 
routine procedure, speed up field work and 
ease the going in rugged terrain. 


Watch for exciting new automatic vehicu- 
lar handling techniques, now being field 
tested. 
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From stern to bow along the 420-ton, 174-foot TEXIN: 
the electric-powered seistnic cable reel; shooter's platform; 
storage racks for charges; starboard firing reel; 

auxiliary cable reel; jet-engine launches; instrument room; 
radar room; and bridge. 


ATION 


Remote coastal areas. . . far distant from supply ports 
. +. are the scene of many seismic surveys. These areas 
are best worked with the single-ship method of explora- 
tion. GSI’s new TEXIN, similar to her sister ship the 
SONIC, travels to your prospect under her own power, 
arriving equipped and provisioned for an extensive off- 
shore program. Auxiliary, jet-engine launches proiect 
data-gathering facilities into water no deeper than 18 
inches ... giving you a complete marine survey... 
from deep water to shoreline. 


Write now for Bulletin $ 60-1, 
describing both the TEXIN and SONIC. 


Geopuysicat Service Inc. 


A TEXAS INSTRUMENTS COMPANY 
9CO EXCHANGE BANK BLOG. © DALLAS 3S. TEXAS 


Please mention GropHysics when answering advertisers 
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